Last active
March 14, 2018 23:10
-
-
Save jhamman/b804a64a4ecbf1f26fd17f6494c5cbe8 to your computer and use it in GitHub Desktop.
ESMF Unstructured Mesh Regridding
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 12, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"%matplotlib inline\n", | |
"\n", | |
"from geopandas import GeoDataFrame\n", | |
"import matplotlib.pyplot as plt" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 6, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"gdf = GeoDataFrame.from_file('/gcs/camels/basin_set_full_res/HCDN_nhru_final_671.shp')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 7, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/html": [ | |
"<div>\n", | |
"<style scoped>\n", | |
" .dataframe tbody tr th:only-of-type {\n", | |
" vertical-align: middle;\n", | |
" }\n", | |
"\n", | |
" .dataframe tbody tr th {\n", | |
" vertical-align: top;\n", | |
" }\n", | |
"\n", | |
" .dataframe thead th {\n", | |
" text-align: right;\n", | |
" }\n", | |
"</style>\n", | |
"<table border=\"1\" class=\"dataframe\">\n", | |
" <thead>\n", | |
" <tr style=\"text-align: right;\">\n", | |
" <th></th>\n", | |
" <th>hru_id</th>\n", | |
" <th>ann_P</th>\n", | |
" <th>lon_cen</th>\n", | |
" <th>lat_cen</th>\n", | |
" <th>AREA</th>\n", | |
" <th>elev_mean</th>\n", | |
" <th>ave_T</th>\n", | |
" <th>july_T</th>\n", | |
" <th>Perimeter</th>\n", | |
" <th>geometry</th>\n", | |
" </tr>\n", | |
" </thead>\n", | |
" <tbody>\n", | |
" <tr>\n", | |
" <th>0</th>\n", | |
" <td>1013500</td>\n", | |
" <td>0.0</td>\n", | |
" <td>-68.56551</td>\n", | |
" <td>47.01169</td>\n", | |
" <td>2.303988e+09</td>\n", | |
" <td>277.493490</td>\n", | |
" <td>None</td>\n", | |
" <td>None</td>\n", | |
" <td>647993</td>\n", | |
" <td>(POLYGON ((-68.35649871820158 46.9031143188376...</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>1</th>\n", | |
" <td>1022500</td>\n", | |
" <td>0.0</td>\n", | |
" <td>-68.07313</td>\n", | |
" <td>44.79691</td>\n", | |
" <td>6.203873e+08</td>\n", | |
" <td>103.604200</td>\n", | |
" <td>None</td>\n", | |
" <td>None</td>\n", | |
" <td>312624</td>\n", | |
" <td>POLYGON ((-67.97836303701675 44.61310195937608...</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>2</th>\n", | |
" <td>1030500</td>\n", | |
" <td>0.0</td>\n", | |
" <td>-68.14985</td>\n", | |
" <td>45.80527</td>\n", | |
" <td>3.676155e+09</td>\n", | |
" <td>174.433896</td>\n", | |
" <td>None</td>\n", | |
" <td>None</td>\n", | |
" <td>662248</td>\n", | |
" <td>(POLYGON ((-67.83991241436119 45.3661422724800...</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>3</th>\n", | |
" <td>1031500</td>\n", | |
" <td>0.0</td>\n", | |
" <td>-69.58119</td>\n", | |
" <td>45.23568</td>\n", | |
" <td>7.665447e+08</td>\n", | |
" <td>304.743490</td>\n", | |
" <td>None</td>\n", | |
" <td>None</td>\n", | |
" <td>309614</td>\n", | |
" <td>(POLYGON ((-69.33810424852713 45.1231651307486...</td>\n", | |
" </tr>\n", | |
" <tr>\n", | |
" <th>4</th>\n", | |
" <td>1047000</td>\n", | |
" <td>0.0</td>\n", | |
" <td>-70.16213</td>\n", | |
" <td>44.98744</td>\n", | |
" <td>9.049562e+08</td>\n", | |
" <td>379.779978</td>\n", | |
" <td>None</td>\n", | |
" <td>None</td>\n", | |
" <td>310157</td>\n", | |
" <td>POLYGON ((-70.10847473135688 45.21669387774366...</td>\n", | |
" </tr>\n", | |
" </tbody>\n", | |
"</table>\n", | |
"</div>" | |
], | |
"text/plain": [ | |
" hru_id ann_P lon_cen lat_cen AREA elev_mean ave_T july_T \\\n", | |
"0 1013500 0.0 -68.56551 47.01169 2.303988e+09 277.493490 None None \n", | |
"1 1022500 0.0 -68.07313 44.79691 6.203873e+08 103.604200 None None \n", | |
"2 1030500 0.0 -68.14985 45.80527 3.676155e+09 174.433896 None None \n", | |
"3 1031500 0.0 -69.58119 45.23568 7.665447e+08 304.743490 None None \n", | |
"4 1047000 0.0 -70.16213 44.98744 9.049562e+08 379.779978 None None \n", | |
"\n", | |
" Perimeter geometry \n", | |
"0 647993 (POLYGON ((-68.35649871820158 46.9031143188376... \n", | |
"1 312624 POLYGON ((-67.97836303701675 44.61310195937608... \n", | |
"2 662248 (POLYGON ((-67.83991241436119 45.3661422724800... \n", | |
"3 309614 (POLYGON ((-69.33810424852713 45.1231651307486... \n", | |
"4 310157 POLYGON ((-70.10847473135688 45.21669387774366... " | |
] | |
}, | |
"execution_count": 7, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"gdf.head()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 16, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<Figure size 720x432 with 0 Axes>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAACkCAYAAACHKtXIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvFvnyVgAAIABJREFUeJzsnWd4XNW1sN8zXaPeu6zi3nsDd4oNBhtsajB26BACXEJyST4gkFwgkAakkNA7BoxtsAEbN2zj3qtc1HvXSBpNP2d/P+Qm1EbWjJrnfR4ePOfssmZ0Zs3ea68iCSHw4cOHDx89H1VXC+DDhw8fPjyDT6H78OHDRy/Bp9B9+PDho5fgU+g+fPjw0UvwKXQfPnz46CX4FLoPHz589BI07jSSJCkHqANkwCWEGCtJUhjwGZAM5AA3CyGqvSOmDx8+fPhoi/as0GcIIUYKIcaeef0ksEEI0Q/YcOa1Dx8+fPjoIjpicpkHvH/m3+8D8zsujg8fPnz4uFgkdyJFJUnKBqoBAfxXCPGGJEkmIUTIBW2qhRChzfS9D7gPwN/ff8zAgQM9JrwPHz58XArs27evQggR2VY7t2zowGVCiCJJkqKAdZIknXBXECHEG8AbAGPHjhV79+51t6sPHz58+AAkScp1p51bJhchRNGZ/5cBK4DxQKkkSbFnJosFyi5OVB8+fPjw4QnaVOiSJPlLkhR49t/AVcBR4Gtg8Zlmi4GvvCWkDx8+fPhoG3dMLtHACkmSzrb/RAixRpKkPcDnkiTdDeQBN3lPTB8+fPjw0RZtKnQhRBYwopnrlcAsbwjlw8dPMZmtZBZWICExPC0WtVrFmUWGDx8+zuDuoWivpcaRh0Edil4d2KVyZNZtIsowiEBtTJfK0V05dLqIJ/7RYNULDfRjQJ8onr//WoL8DV0smQ8f3YdLJvRfCIVjpu/IMe9BEfK568G6pC5X5i7FTqUtg+Omr5AVR5fK0l3JKz0fhFxdZ2Xn0Vzyy0xdKJEPH92PXr1CNzkK0Uh6zK4KIg1pJBnHsaP8bxyteoeR4UtI8J/Q1SICYJNrqHbkkGPeRkbteuYn/Rt/bZsup5cUd1w9hhF941CrVVhsDr7bkU5IgF+HxrTYHBzPLmH7kRwGp8Rwxbj+7eovKwpq1SWzJvLRA+i1Ct3iMrG59F/kWw4Qoo2n2lGAURNCjD4agzqYcH37vrzeJEAbxdXxz7O19G8cN33FltK/ckXcs2hVPnPCWeqsdt7+eif+fjpunD4cgG2Hsrj5ilHtHquipp6SilocLpnaehuBRgPTR/dt9zivfb6V6joL982bBEB4sD8rNh9h4Yzh6LS99qvloxvjVqSop+hIYJEQgmpHEUHaSDQqXaPrAgWVpG62X7W9gBO1G9hT+TFqSceo0IVMjlpyUTK4i8VVBQiMmvB29XMpdraVvUZ+/S6MmnAGh1zPwOBrvSNkD8Fqd/LeN7sRCI6cLmbvifxz9zRqFddeNpjHb5uO0dDwTJgtdjILK4iLCCYyNMCrsq3fc4pn3vwOp0smLT4cnVbDidxS7po7gdEDEhk/OMmr8/u4dJAkad8FebRabtdTFHqWeS/HTZtICRiDgosIfTJGTTAbS/5NonEEY8NvbLGvorg4YvoGf204fQMvv1jx3UIRMl/nPUCl/TS3pHxGgDa62TYSLXtpmJ1lFFkOIFDoHzTbK94cmaWVJEeGdonJoN7uwF+va7sh8Mpnm/l47T4AZozuy+7jedRbHShn7quAUf3jKas28+r/3EBRRS2P/G05eq2alS/fTUSId5X68h8O8/JHG5AVQWx4EDHhgVTXWTHVWbhv/mRumjnSq/P7uDRwV6H3mH1hpD4Zs6uKg6ZvUaFiROgcwnWJ9A2cTKg2rsV+RZbjqCUtI8LmdYqcEipGhS+h3lVGvau8WYVudpZg1ISjkZo3qQRoo+gffLVX5UyLbt/uwZOsP5pBTEggE9IS22z7wPzJ9IkJ5f1v97D/VAE2uxMCVIgwNXaXC0mS2G8rx0+SuPOlT7DWNhwq250y9734Oe8+fRvBHbS1t0ZIgIHwYH/Kqs0UV9ZSVWthQJ8oNCqJlNgwr83rw0dz9BiFbnZVYlD741AsDAu5mkHB09hb+SVpgRMJOaPQhVCQpPMrzmJrOnqVP1tK3yJCn8SkqDvRqPRelVOSJPoEXHZebmcp64p+S4gumclR/4NeHUiQLr5RHyEEDsXc5d42ncW8MYPdbmvQa7lh2nC+3HSY4vJaZD+wBAlwuZAkSAgLZmJqIkcKSqix2ggP8MdcayMmPJAKUz27juVy1YS2E8LVmK1kF1WRGh+Ow+kiJMAPjaZ5M96FTBqWwlej+vLOqp2czCtnwpAkckuqWTRnHDFhl8bf00f3occo9Fi/AcxPfOrc64y6nVQ5ihipiTxnknAqNnRq47k2YbokTpg2UecqYUTodV5X5heyt+K/5Ju3I1BIDbyCoaG3Npn/rIuiWqXDqVh6vUJ3uFw89N5XGHVanpo/kwCDHrUkoXfjAPGjZ+/AVGfh0U9Wc7S4lHmjBvPUdTMAUJ0xG20+nMkTf2/wVddq1Hz03B1umZQyCip49G/LKTeZ0arVOFwy/7toFgtnNomna4KfXgvAffMnt9m2oxQWVxMdFYxG7fOs6e4IIagsriYirnN3aT1Gof+UakcxJ2q3MjJ0LlGGFIBGyrzYegKX4mBE+FyGhV5z7kvfGZidJZyqWYUQAn9tJEerl3Kw8l2mxj5NauD54No6VzFGdQRqdNTYC6i0n6ZPgHdt/F3Jyn3HKayupazWzIwX3kSjajhHmDE4lceuvow+EU2yLzciJNDI+/ff3OL9acPT+L/7r+FYdgnJMWFo1Gq3Ao9yi6sYkhrDwVOFCAGvPnQtw/q2bMZrDxaLneoaC/Gxrb83d9Bq1D5l3s3ZsmwH3729gYqCKkJjQljyx1sZPLHzPOp6zKHoTzG7qpCQ8Nc0/0WRFSd7q75kTNiNSJJETt16Yo1jMWo6x7/bLtdhdVUSok9GVmzsKX0USdhJ9J9MbNA9jUxDZ1GE3KK3Tm9g47FMaqw25o0ezLbTOTy9bB11Njs2p4tBcVG8dc8CQoxd56ppsTk4nFHExKHJHhtz++4MJCQmjU/z2Jg+uic2i53bEu7HbKoncWA8L69/hrCYEI8sJnvdoehPCdA03crUOsqxyCZi/PphcpSQV38YndqfQUFTiDQM7TRlDqBXB6JXByKEoMD0AmrnWgAKa35AJakJNc5Gr2lsS+/Nyhxg5pDzSm3KgBQ2/vZebC4Xh/OKGdUnzi3TizuYLFYyy6oYkxzfduMLMBp0HlXmAJPHt9+/3UfPJC+9AJ1By5irRvCLV3/e6eYW6KErdKdiZ13xv4g3DiZEG02C/zCOmtZxzLSJKaEjiFLWozLeQbEcQaJxWLOrYW9jc+Zgc2VjcRynsObP566rpUA0qjDscgFaVQRJoU8T5j+30+XrCh7/ZDUbj2eSGBZCiNGAVqPm7qnjmNg3sYmtO7eimsyyKqYPTEWlap/bZr3dwT3vLseg1fD2zxe0u39HKCo28dmKPYwanoRKkph6WcvbbSGE2y6p36w9zK592Tz+iyupqbWi0ahaNeOYaizkFVTRNzUSnU6LRq2isLgavU5DRHjvPqvpSqz1Nvy8kF+oV6/Qqx2FxPkNIkgbgclZQh9pFCNC5zAidA4AQtwKqEjqomx8DrmMI8UzAQVQY9D0RwgrLqUOWZiQ5ToAnEo5wX4zu0TGrmBPViFOWSGrvOrctV2Z+YQHGLl+1CDmjBiAXqPhuZXrSS8sI8jPwMC4SOJCgto1j79ex0MzJ/KXNVvo7EegpKyGk6dLWLfpGL96uHnX05XfHOD4ySJ27c0mMiKAP/6/+cRGhzTb9iwzpgygb2oUOp2GkGAjGk3zi5T9h/I4mVnCm+9vQVYEWo2aAX2jMZtt5BVUERbqz9yrhzNlcn/6pkR1+P36aIw3lHl76JEKPcqQSpQhtcX7kqTGWfcKIKMN/FXnCQZUWdagKJYLrsjYXKeabRtsmIZaZWz2Xm9Er1GjliTkn+wKK80W3t26j3e37mt0/aphSRi12ouaa1xKAvdOG48QdKpSHz2iD//5+yKcThmttnkTmhBQXlFHclI4Uyb3JzoyuM1xjUY9A/rFYLM52Xcol6SEMBLjQtHrz38++YVVHDtZiNXmJCzUn/JKM06XTFFJDXfcNIH0k8XsPZhDdn4lq577khmXD+AXd83o1B3MpcKadzZSW1nHNfdeQUCIf6fN2yNNLgDCdRql/kMkTT8k461IUuMvvlDMCMWESpPgkfnckkkISuveoaj2NWSlps328cG/Ji74F50gWfeioMrE6xt2sSe7gMLq2mbbaNQq3r1nIaPbaQfv7VRV15NfVE1WTjkHjuQxa8pApl02AIB6i52fP/weS9+6j6ISE0uX7yYxPowhA+MYOsj3OXYWpvIabk98AKfDxc1PXM+9Ly/q8Ji9LvT/p7gqrgXXSQBUAb9BFXAfAIqrANm6ApV+Empdm++/XQjhbPLD0RyKolBQ8yKldW82uaeS/FCEFYCogCX0CXvWozL2NKrNFk6WVnC8sIx1R08zbUAqQxKiGN0nHn+De+kBejMWq4P8wiqCAv2IjW68kt+xJ5Nv1x0hMMDAbx6ZDUBZeS1Rke0zUfnwPM/e+DLbVu5BpZL4xWt3c/1DHYv87tUKXbhyEPb1CKUKlX4WaEciXeAhItt3odKNRJI8F0ikyGXUVMwnMOw9VJpkQELVinJXFIVDReNwKdVoVMEowkGo32xigx/iVNliHHIBQ2K+xahzP2qyN9GeA8HeiktWOJZTQkVNPWP6JzRJB2yut7N1xymMfjrKK83MnjWUAP+mz3S1qZ7QTtzW+2ib7V/t4fc3vAzAzNsv58kPH+nQ8+5xhS41aMy9QKEQYq4kSe8B04CztoUlQoiDrY3hyRV6WzhthxBKIRr9dFRqz9ipzyohRdgorX2H2OCHLmocqzMDk3UjMYHN+6O7y4HsQlbsPsaVw/sxsX8SWnX3dHustztQq1S4ZJlF//0cf72OnIpqhsRHE6DXMSA2kgGxEQxNiCGrrIphiTGU1pipt9sJ9NOTGNb6gWF7sDlcbNx/mn4JEfRL6Lqc80II/uffX7PlcBbRoQHcMn0kS2aPOy+nzcmTf1zO6cxS/ufBK7li2qAuk9VH81QUVhIUHoiumZ2kEII/3/UvNny4hYiEcF5a9wwJ/WIvei5veLk8CqQDF+7nfi2EWNZe4TyBEA5A28qvnhVbzaOotWPwj1jqkTnPzqWSDNTZ96K3fEOYsf3pbf20ffHTdtw/udZqZ8XuY+RX1jC8TyzBxgaFbrbZMeq0nRod2xw55dU8u2I9h/OLkYXAJSuN7v94KgeANUcaHxpr1WqcckNVqQdnTuDhKzseVi+EoMJcz6Mfryb3WDlCFtx73UTuvHpsp2ecNFvt7ErP47KhyWjUKq6ZMIjBfRoncdNq1dwyfyyr1x7mdFapxxV6aXktew7kcM0Vw3yHohfJ64+/T2L/OH729AK0usa7dUmS+M27D/OLV+9CpZLw82KCuAtxS6FLkpQAXAs8DzzuVYncRJJat69qDRNRR65HoMVe9280hitQaz0Xgts38nUqzJ8jK1bUqub/WELIjUxBniY2NAiNWkWQn54qs4X3tuzjaH4J+7IKQYKbJw7nN3Onuq3Y92YVEB5oJCXSMwERoQF+BBsN2F1y241p8EYRgnPKHOBUSYVHZDlaWMrqQyc4VVrBk7dNIzk8lNH9ErpEmQX46Zk1uh8A8y8bSnpeGVE/SfOrVquYNC6NCWNSsVg9X5awqMTE0pV7mDa5P4EBvkIqzZFzLJ+g8ADCYhr8/YUQLPvrKooyS3jk3/dyy2/msfSllWxbsZvpt1zW7Bj+QZ3rxeaWyUWSpGXAi0Ag8MQFJpdJgB3YADwphLA30/c+4D6ApKSkMbm5uR0WWrGtA+FA5dd0dSyEHcV1CpU6FUnVYFeUnVm47NvRGq9Fpep4Tg13Kap9D7NtK1GBd6KWjAQaxrXdqT3jV9Wy4WgG+7MK0WrVfHvoZKP7Kkli53MPnSv+0BplNXU88uFqFo4fysLxwzwmo83pYsvJbCx2B+EB/gQadEQE+hOg11Frs1Nvd5AU3mBSqaizsC+ngIKqWkYmxTI4PoogPwM6N7IeusNZ04+hG1UTKqqsZd5T7/DSfXOZOapzokpPZ5Wy+vsjPHLvTNS+3DDU1lkJCmy8KHM6nDw+7ffc9/Iihk0ZhBACRVF4bsFfOLDhCB9k/ovQqNbdTX9csYuqYlOHD0TBgzZ0SZLmAtcIIR6SJGk65xV6LFAC6IA3gEwhxB9aG8sTNnQhBIrpQST9DFTGW35yz4ls/xEh7Kh1o1Gpo5Cdx1FpUpFayD3ubTIrfolDLiXEbyaxQQ94ZQ5FEVTXW/jVh99gc7k4WlAKgAS8de8CJvRzr3LOhmMZ9IuJOKdgfXQO+eUm9BoNUV6usAQNB62/fvYLdDotrz5/S9sdejkuWeFXzy3jleduQpIkbBY7BmPDwbOiKKhUKlxOF39a9BoPvfJzwmJCyTyUQ/LQRNRtnFlVFFbyy4m/4x87XyAivmP1Bzyp0F8EFgEuwECDDX25EOKOC9pM54yib20sj3m5CBnkYqRmfMxl5zFc1q/QGu9Apen6EmCFNa/gkquJDrobg5fkccoys59/h2qzleSYUE6VVBAfGsQ9M8Zxw9ihLUYVtkR2eZXHzC4Ol3xRK2whBDVWG3U2O0adjvCASycAy9Ocyixl7aZjqNUqjp0oIiu3nGeeuI5JY1sOzruUyM6vICUxArOpnneeWsoj/7y7SRubxY5Wp0HdjmfZYXNgKqvBL9CPwA7+WHvsUFQI8Vvgt2cGnU6D4r5DkqRYIUSx1HBSOB842iGJ3UQoNcjVD4JzN6rQD1Dpzx+YKXIJLstXyLavkdTJqFwJaAxTO0OsFokLehRJknDI5ThcZeg0ng+31qrVPHz1JPpEhfLVvuNYnU7qrHaGJES3W5kDxAQHoiiiw/blradyOJhXxC+vaN+hZm6lif9u2sXK/ceBhiCjOyeP4p5p47s0G2NPpU9iOBnZ5cREBTF1Uj9ee+FWn6nlAsL9DbicLgJC/Hnkn3fjsDkwmyzoDFp+XL6Lq38+g+/e2sDACf0YNKGf2+PqDDqikjrXk6pdfug/MblsBCJp2NkfBB4QQphb6+8Rk4tSj3BsRUjBKK4MtP6No7AUxYaEgtSNQurLzF+Sa3oZCTUDI98gQD/Uq/MdyCni+ZVrOFFUww1jhxDgp2f+2MEMiHX/4XLJCiaLlYjAi/dvXncsg1qrjQVjW3+/GaWV/HXNVnIqqrE4nJTX1Tfb7j+L5zN1QEqrY3194DhjkhOID/UF1/honbO6b8/3hzi5O4M7nlqAJEns/u4AtZV1jJszEgmJoPDAc227KnbCK8m5hBA/AD+c+XeXZJWSVP5IhtkNof3N/H6oVN1vBeenTSUt7AXUqkBMtq1eVehO2YQw/IH75wSyNyOUnNIK9uTY+M3c9u1UBB0P/BmbHM/bW/YAcLyojIgAI1FBjbeeVoeTRW98To3V1uI4KgkemDGBkUmt+/E6ZZnM8ipsThc3jx/eIdl7OqXltWzdeRpzvZ16i51b5o0jItz7NvqewicvruDr19fitLt4YfWTLHp64bl74+eMatK+pwTBdZ/j/nYiqQLQGK4891oo9SjOA6j13a/iT6D+/AMSbBjv1bnscjUx/rNxyh8zZ6yVkrKxjEic1O4HUqtWd9huHervx+NXTwHgj19tYHBcNE9dP+OcLPtyChrMKpJAr1WjVauJCPAnzN8PJCipMaMIhVqrne+PZXC4oBS9Rs3guCiOF5Uxc3AaN4wecm4+q8PJsj1HsTqd9I+JYGSSZ6oOeYOaehsPvvolN00dzg2Xe86r6CxREYHs2JvFngM5AKT2iWTOLO/uDLuK9NPFnMgoYczwPiTFu3f2s2fNQaqKTQDNBgb1VHqsQv8pLutSnOY30Id9ilp76R72BOhSCNClEBswr0FxRrfdx5uctcN/+uBtja7vyMjj3veWo1xg8rO7ZMx2BzmV1U3GySirJKOsEoAN6ZlIEswZPqBRmwC9nukDU6mqt9A/pvOjQEur68gurmLi4D5ttlVJEBMaSFpchFdk2bkvG4RgcP9YXLJMeEj3MUF6mkH9YhnUzijM51b8mpLsMqL7RGA84yu+7sMt9B+bSp9BnZfQz9P0GoWu1gxDEvVIjg1wiSp0p1KHCjVqlbHRilwIgRB0i4hAi8PBPe8s52B+cYfGGRQbxeyhjQPFVCqJ5xdc5ZED3Ythx/Fcckur8dNryS6pYtqINEJbiBAMNBr424PXe02WcSP7kJFVxpCBcfRJDCP8Ai+LOrMNf6O+WzwPXUVQWABBYY1NUOPnjMRu8XwQV2fSI5NztYRQ6kDy75IKRV2NEDICpdmEYS99/QPXjxlMcmQofrqLyy/uCcrrzCz458dUmC1tN76AuJBAhIC40CDEmRQCUUEBzBqcRoifgakDUrqFjdPudPHLf67kmTuuYMnLn2G22nl60ZXMGTewWynPg0fzGDEksVt8Zt0dh83RLUwy7h6K9irNJ6kCL0llDpBR/Qo/5l+J3VXe5N4VQ/vy8Htfsf1Ux6N0O0J1vbVFZR4RYMRPq22UVyVAr+OZ62cyLCEGSZLYl1PI/twiDheUMDgukqp6K//csIM73/qCv63dSnpRWWe9lWbRazU8MHcitRY7bzy+EINOw9PvruHW5z/iHyt/5Hhu6bm2sqI0yW3TWYwcmuRT5m5iq28S/N6t6VErdOE8irBvQSBQGa5C0vRrMCc4tiPpJnv8IT0bKSYLBRVSt/4SyIqdXUU3MjD8acL8Jja5X2OxEdxFPtyH8oq5573l1NubbmeTwoJZOHYot04YScCZlZDZZienwsTQhIYDgLMmo2NFpfxwIgudRsO9U8ehUkkIIbA6XezIyGPVoXQKq2oY1SeOmYPSmJDatavQnJIqfvbiJ1jtTgASIoL5+v/u4rkPvue73SdIjQ3niZunMbpf97TZnjhdwsB+MV6fJyu3nJSkiHb9rYQQlOZWEJN8cWclddVmco7lM+zyhqRnF0aIdkd6ZU1RodSAKhjJMA9o+JJIkoRAdea1Z7dG3xbv5IOcNVwdO54J4YMZGtx9bfNqlR6DJg7/Fs4PukqZr9h/jOdWbsBxJuGWSpKICQ5kVFIsiy8bQ3JEMAGGxrIFGPTnlDk0/I0lCYYlxDAsobGCkSQJo07LrMFpzBqchhACu0umzmrr8h/g5Jgw3v/NreSUVJGeX8bqHccprKjhZ7NGM3v8QEL8DQxI7L51Pfuldo5suw5kk5LUvsPhd59dRlxaNLPbqdD3rjvMqBlDKMkuZ893B88p9OwjuQya4LnkfV1Fj1qhdzZ7q06woXQf96TOJVzfdt3HrsbsOI1eHY1W3X2Cag7lFfP9sdOYbXauGtqfcSnx6DQ9ah3hMSpq6gn2N6D1ULKx9nIxK+HOZNeBbMor6ph7ZdsxBDnHC+gzKL7Je1EUhdzjhUQmhuMf5NfkvtlU36k1Pj1Fr6xYJOybEdZVSAEPImnSPChZ1+Kp6j3/Or6FcL0/t6aN8YBUvRurw8l72/ezYPSQJsFOvZX8wioS4kK7rUJ3yQo2m7PZqkxHThTSJz6sSVbEs5zYm8mKf33PkW0nGXvFMGqr6ggKDeDhv9/Z5FDz9V99wOwl00kZ1vW5ntyldx6KKnVUOc3UOgrPXRJKPUJuehDYU1CEgyzT31GEs8NjORWZtYXpVNhazcDgAzBoNRSaatvtceNtvtp+zGuHpYnxYV2izE1V9bz56vfU1lhbbffOp9t49a0Nzd4bNjC+RWUOUJRZyqxbJ/PzZxeQm17I7MXTSBwQhyI3XbAuenoBhl6aA75nrdCF4NXTv8OoDmBh4r2E6bqv/fEs9Y5TDelzDZd3ypfJ5LASpDWg6qarMB8t43C6eOK/q3n+7jkE+nXOAd0Xq/bhZ9C6ZeZwh5paK8VlNZRV1DJpTBparZqczDJeeupLHE6ZN5Y+hLqZhHFL393K3uxiXPU2pqZFc2MHcoh/+Y/vUKvVzH/oqhbbOB0utLrzpr/MQ7mkjWg7IKyr6JUrdIHAIps5ZT7c1aJgdllZVbgdh9z6ytopV5JZ+UdqbDs7Ra4QnV+vVeZ1NjtldT1z91FdZ6Gg3NRqG51Ww2sPzz+nzGstLee38RRajYoV3xzAavNMQE1QoAGXS+aH7afQahvOCpLTolj80CxmzxuFdIE//ub1x6itseKwObj+5nEsmTuK9He28Pnf1/Dw9D9wcl8WFUVVHNyczpoPtmAqr3VLhhsfns31D1zRapsLlTmAy+GkMKOkne+2+9GjTqdUkoqFCfdy2LSTYG0YtU4Tp+uOMjxkPFpV5zr/ayU16bW5WGQbtyS1nKcsxG8SYcYZmGzbCPGb5FEZSm0mDlXncFXsSI+O2xKnSivoH+2dUPWWKK01s+FkJpnlVaSXlLE/v4jbx47grsljqHc4iQ0OIMjQ/bfPoYFGQgPdC78vqqzlz1/8wMyRfZk7YZBXd3YnM0rpmxKF1kNVnCRJYujAeAakNfZGmjilPxOnnPciqTfbMNdaeeeZz1i/dAcjLh/IUx88SFK/GKbOH8uaj37kyI5T7Pn+CCq1CqfdRXV5LVf97HLCY1ovwNLgFdW+z6zv6FQQosnKvafRo0wuilA4WXeYQUENCuztrL9wtHYv/QKG8lDfp7DJFiwuM2H67mWKEcJFneMoQfrOUbzeYH16BkIRTBuQ0qleKrKisOb4aarqLfyYmcvm09lAg/ujIgTRgQH8+sopzB02sNNkao72HmxbHU5cLpnSajMqlURq7PmKNtuO5TA8JYbATnA19dSB/MVgs9i5f9IzlOZVMnvxFB752yLyThShN+qJjA+lqqSGnWsOsfGLnQyZ0Jfs44WMnTWUaxZPxfCTg1NLnRWdQYtao27X+7HUWTFeYJt3OV1oulGJwrP0Si8XRSh8UfAWC+J/jkalZVPZar4t/gxZuHhq8GvY1REVAAAgAElEQVS8n/MKtyU9SIyhewZq+OgYDllmw4kMvjxwjGqLlWPFDZGhIX4GXr9tHqMS3c+uWFZnJkCvx9iBVAiKIrjzpU9IjglDCMEv508hJiywzX4/Hs3mQEYhkiRhtTvJLzcRFx7Ek7d2SUbqFhFCYLE68PdSwM3Wr/by9ZsbuWbxNC6fN6bVlbEQgj/c8S9qKur4y3f/26jw+Y9f7eXL19bwP/++m8Nb05l7j/ufY1f+oLWHXqnQZSHzQ9lqXMLF1TELAKiwl3Cy7jBDg8firw5Eo+q6XCXeJMtcQmqA96P2fsqG9EwGxESQENp9/PAVIXDKMu/t2M+qI+lklFdh0Gj4z+3zmJjinivaX9dvJczfyM8ntd/F0+F08fnmQ9Q47OzNKyTeP4A1O06ycOrwbqeULxYhBKezynj/8x2MGZ7EjdeOvuixykpqCAg0YGzGHfFCnA4X+zYcZdT0wej9mppQy/IrCY8NaVIGTpYVKouqCQrzx+Df/c1vF0OvVOhnWVX0CbOjF6JRaXvEr2tHkYXCbdv/xJ0pVzA3zrv51H9Kdb0VJAg1tuwy1pU4ZJld2flo1WompiS63a+0zkyY0Q9tG4V+W6LEVMfegkIqzBbmjxiMkAUul8ypgnJ0Wg0Dk6IIMOh65POZk1fBus3p3HNHg2dWlamesA4E43y7fB/VVWZ+ds+0Zu+bymsRQhAa1b5Fw9nUHJcCvTL0/ywzoq7D5Kzk84L/oJV0jAmbypjQKV0tVqvsqkhnfPjAi/qCH6rOotpex7L8rVwVPQqduvN2IaH+3VORn0WnVjOlb3K7+0UHdiyYKCYkkLkhTe32VXUWEiJDWLf/FCfzyxmZFsfUYan46Rv+ZiazjV+/vQq704Veo+HJW2aSFtu+ivCbDmVSbjJTZjKTnlfKrxZOa2SD7yj7juRxOrvsnLmlI8pcdilcPW8UDrurxTZ2qxO/gPabdb78x1oWPjK7R/5oegu3FbokSWpgL1B4pqZoCrAUCAP2A4uEEF5PJuxUnGwqW0WZrQhFyJyqP0ydy9SqQi+2FhDr17V29YtV5gB1LgtDQ5JxKjJZ9SUMDHJ/JeoN6u0O/PVdn1K0O3I2N8v8yUORFYVvdqXzxrc7WXLVOIwGLQv/730q684HM9mdLSu6lqius/DNzuNoCm0MSYjyqDIHWHDtaBa0YWI5daKI/gPbPrNwuWT0Bi1+xpafl+ikxvLX11rIOlpA8uB4Alv5Mbnp0Tltzn+p0Z79yqNA+gWvXwL+LoToB1QDd3tSsOZQhEKxLY/r4m7nqpgFTAibRR/jAOL9mi8cLISgzFbCprI1bCr7ztvitUpHVhHB2gCssoNcSym7K095UKqLozcp81WH0zlZWuGVsdUqFddPGsKjN0wh2N/A5z8cQiOpOPskXDt+EAMT2u+RdePlw3jzsZuYMbovN9/mWVfY1qivt/PH331BbY2FlDT3SmHpDe3bTcqywptPf8HBLemc3Jt9MWJe0ri1QpckKQG4FngeeFxq0E4zgdvPNHkfeBZ43QsynkMlqUgyNuRwSTSmkGhMASR+rPgeRSioLsiFLoTgVN0x9lVvZ2jwaIYE91yXwZGhqbw25kHsspNAbVMTSLXdgp9Gi6ETTTHdlbMl7dwNrgo2+vGbFd/x1QOLPCpHdZ0FP70OwwWeG6P6xrNy21FG9Ytn/cHTpMSEXXThC51ewx1Lmu5KhRA4HC70es8/C/7+eq6eOxK9XnsuaMjTqNUqHnt1sVfGvhRw1+TyCvAb4KxPVjhgEkKc3S8WAPHNdZQk6T7gPoCkJM8mw7G6LJicVRRac7HIZgI057MMSpJEuC6SGxMWYVB3bztwW1hlB98V7SFaF8Jl0UOa3N9Rns2AoCjSgtxPJWqTnRytLmJsRPvCnZ2KjCIEenXnH7/U2mzIimjxgNYpy7ywdjM7svN4YMp45g8f3Op4Dlnmtc3b6RcVgawojYprdBSjQccb3zREB18zfiBhQf4M7hPNF8/cCcAzdid6D/s7f/bxDqoqzSy+eyp4KXPA+Mn9LqpfdaWZ0PDmzy2EENRU1BES6ZksoS6nq93+6L2FNp9gSZLmAmVCiH0XXm6mabPuMkKIN4QQY4UQYyMjPVe41yZb0KsNbK1Yg4LMd8Wf81OPnQhDdI9X5gB+ah03Jl7WrDIHuCZhSLuUOYBBrW2XMrfLLjJMFWTVVlFsqWnXXB3lZGkFj3yxilmvvXMusKg5tGo1v79mJmt+sYR5wwa1OW5GWSVmu4NDRcVkVzQtTN0RiiprGZgYxYDESB7991dU1NQ3uu+n13q0LJ0QgoW3TuC2RZPbdA9sifKyWlYt34eieNbz7eTxQn6+8F+kHyngu5X7ObTtFPIFCcgkSWpWmRdll/H1mxuRZYWNn+9k/dLtbs2n0WouSWUObrgtSpL0IrAIcAEGIAhYAVwNxAghXJIkTQKeFUK0mlHHU26LilAosuaiUekQQmFn5UbGhF1Ool9qu/+QTsWFVtUjnX06nby6KpICw6i2WwjVe7+KfHFNHauOpJNZUcXlaclcltaHMA+7TxbX1BFo0PHMNxt49pqZBBkMHCgo4oNdB3j+uqs6FHiUW1pNfEQwapXExgMZaNQqpo1oX9rnL7YeYnhyLIFGPRa7k75xERw+mEdUdBAxsa2HwF+IEIL0Y4UMHtr5zgFOp0xhXiXGAD0lhdWEBhpIdKMSUnFOOdnHCqgoqmbK/DEEBPu3GXwEHTuv6q54xQ9dkqTpwBNnvFy+AL4UQiyVJOk/wGEhxL9b6+8phe5Q7KTXHibRmEKYriG3iCIUpIsoE1djNxOsD6DWWU+QtvMS38tCQX2J1j91h7XppzHbHcwZ3L/dSrWy3sKrP2xnQHQk1w0dQJDBwMrDx4kK8GdCcmKzppW16af5w3cbWfuLJRh1Oha+/Qnh/kb+c8u8DptiXv5sEycLygkN8mPasDRmjx3QZpGLY7klDOkTw4n8MvrHRzZazX/64TYCAg1cN9+9oKjqKjNZmWVERgWR1KfjuXhKik3t+jE5izth9bJLxm51Ygw04LA5Kc2vdEv5A+SdKCI4IpDgiLajdXsanZFt8X9pOCDNoMGm/nYHxmoXOpWeESHjzilzAJOzkgOm9mc0lFEQQlDrsCCLzivaW+foeB5uxUvyZtRW8P/2fkNB/fnsgEeLS5uYtLzJ1YP6sWDkkItaIYcZ/UgICea/P+7m99805Nee0CeRj/ce4heff43Ddd5V0Gx3kFlRxdWD+pEYGowkSagkiddvnke/yHCOl3S88HRqXBjH8krZcCCDZz5Yy7vf72m1vRACh7OhZN/AxKgmppmsjDLqze4XLw4NC2Dk6GR0Wk2r5pQtG9P56J2tfLfqAPtb8DCpr7fz4jPLqa5yP+tljcnC6i/38uAdb1Jjavzcn9ibRfqeTACsZhsv3vMGGYcaipnrDFoS+jbvTeOwO/ni1e8aPZNJA+N6pTJvFw0FeDvnvzFjxojuRr3TKpyySzhk50X1t8tOUeOo97BUbaMoivgsd7PXxt9ZmiOOVRWLEkuteODHz8VJU5nX5vIWiqI0uXa6rELYnQ1/64zySjHyT/8QN775kXC4XF6VZePBDDHqob+Jx17/SuxMz3W73870XFFRY250LT+votn31hqbNx0X/+/XS0X60fwW2zz3u2Xiry+uEldM+qPYtuVki+3qzTZRUlTt9tzLl+4UV457Tiy54TXx98c+OHfdbnOIzct3CyGEqCw2iQenPScWpD0qMg7ntTmm2VQvNi/fLcy1nf/d6wqAvcINHdsjQ/+7E6+eXMHqwp0sn/J7ArXetytfiOikxEJ22YVKktCquqYWprewu1zkVJpYduAIGRVV/M/Myxge1/58OQ6XC5VKhbqNtK1VdRa0GnWnFa84i+xSePDut7n7/hmMm5jW4mGs0ykjSbBt80n6pEaS1Ceiwwe3367Yx9GDeRzYk03VznRueXQ2wRGBSJLE/AdmoVKpOL4nkz/d9xZWs423dv2R4DD3oniLssrYu/4I1983q0My9gR6dS6X7oQiFMpsJmL8wrpaFB8XSXZlNd8cPcEvpk5s1w+krCgs/nAZ+/KL0GvULJkwhsdmTPaipBePzebE0I4gH6dTxuFw4X+RHjMAhXmVPP/bZRSWmAhVCUr2ZfLyqidI35NFQt9oqkpruOzaUTx/939J35tFUFgAb2x7Dv9gI1azjYDgzl0gdWd6dS6X7oRKUnWKMs+rLyPJv3vlee8tpISH8vC09kdcqlUqXrz+KmRFsHT/YRaOaupWujMnnxHxMfhp3Veme07lk11SiawIbps+qt1yNUd7lDmAVqu+6OChirJaDu7LIS4hjPiUCKprLWhrzKSOTsM/2Mgtj80h53ghA0ansPaTbRzblUmfgXG8suZJ/M7U+nRXmZvKaz3mv94b8LlZ9AAq7DVsqzjW1WL0GmptNtafysQhyx0eKzE0hOTwUJ68choJIU2zBX558CjbMhsO+YQQ5FaZzkWztsS4/oncNGUEKdE9c9en1WkoLTaxZ1cGB/fnMDQ1AqdLYcFDs0gd0uA2mTw4noi4UK6/Zwajpw9i5k0T0LXzR8dSZyX3RFGnHtZ3dy5Zk0uFvZwIvecCnTxNndNKvcvqM+V4gUqLhX35RUxLS0bfQvWljPJKHl22mgC9niUTRzNncP9m27WE2W5n3YkMnv5mA8/MmYFTVth4KpOtWbm8e/uNXJbadQWJq6vrCQ3tHBfd7Mwylv3je2orzfzhg/vPmbSO7MxgyPhUVCrVJZUG92K55E0uTsXRap3Rs8rcodixyTaCtMGU2kr4umgZo0PGMSZsQmeJ2oSt5UdZXbiTsWH9uSlpqkfGLLPWEeXX+1y6Pj1wmMNFJVidLn4943Lig9vefocbjVw1oC/Q8sFycngor98yj5jgQLQXKBuXovDGtj0smTD6nEulEAKXoqBVq1GEYMWh4/xt449U1De46C07eIyDhcXnxnh1844uUehCCH73q09JP1bIl9/8CrXGu0r0/Ze/YfZtExk+sS96Q0PtgqqyWt54bjn9hicxYFQyOr2qVWW+b+MxVBoVo6a2HfnroxebXBQhqHfVt9lOp9ITpG3YKlfbqzhec5gdlVu9LV6LfJG3hReOfcquyhNsKD3Av0+t4uvCnR3eVkYY3PMc+DrvcKf6418sQghWHDnOM2s28G36KaalJRNsMOBshxmlst7CkqXLyTc1TWWgUalICguhrM7cSOG/u3Mfr/6wnWPFpQAcKizm5nc+5WRZBYoQfLj7AL9b9f05ZQ6N82Skhofx2PT2H5za7c529/kpkiQxeGgCkiSRl1vBnl2ZOJ3uf15CCN5+fSMbvz+Ky9V6v70/HOfUkTzWfraLq26ZwNTrG84CAkOM3PTQFdx43wx0+rbXk6OmD6IoqwxTea3bcl7K9HqTS42zhm+LV2OT67gy+hoSjO4nCPtpBsfOIMtcjCwUgrRGAjV+SEJiTdFxroofgr/G+2lrT9aU0i8oyu1shV2BEIL39x5Aq1ZzpLiEvflF5FabCDf6MTwuhv0FRQTq9QQa9MwfOgh/rZbJKX1IbMbG3RpWp5OfffwFy5fcfu7a+hMZ7M4t4NEZk/HX6SiqqeWXX6zmz/Nnsy+/kKdWr280hl6tQlY11B8NNOj531lTWDBiaLvf884dp5kwsW+H3VTPuiYu+3QnWZllFORV8cp/FqM7E1K/dVM6kdFBDBzcNNeeEIL/vLaO5Z/tpv/AWJKSI0hJi+LKOcMJDTtvwjl1KJc3nltJQVYZD/5hAX2HJaBWq4n5Sd7zuup6AjvJ9NPTuaTdFgssBQgUEs8o70xzBh/kvEGafz/uSL67TSUthEAWMhnm0wRrg4n1c7/4sDc4XVtGvyCfh0tz5FRVM+eN9xkRH0tWZRXVVhvLl9xGangY35/M4G+bt2F3ydw4bBD3TRrfrlwwmZVV7Mkv4NaRwxtd359XxENffM2Xd99OfEgQ2ZXVbDiZiVat4uX1W3Ep53c48SGBFNTWnXs9NjGeT+68ueNv3EOUlpiIjmkI41+z+iB/fXE1dz8wg1sXXdZseyEEp04Uc2BPNp9+uJ20vlE899LNaHUadm8/jSwr9O8bTVVZDTGJ4UTGhXJ8bzYfv7KGP7x3H2qNmvdeWo3NYmfGDWMZMLLrzhJ6Epe0QoemttEqRyWh2jC3VjiKUCixFRPn12xGYB/diDd37uXlTQ0msikpfVg4YggBOj1T05KBhucgs7KK+OCgdrkOtsXJ0gp+zMrh7kmNv2Pv7NhHmdlMcW0dxaYaDpWUNUlDOig6kq/uucNjsniKv7ywirXfHCK1bxQPPXY1er2m2ZX6WfbvyQKVCr1ew0dvbSYk1J87752GXq9h5QfbKMuvJCYqiJseuoLivEpSB8ehUqk4tP00gSFG1Bo1JbnlFOVUcsO90yktqOLzf61nwX0ziEvpvg4LXcElr9B9XBq4FIU9eQUsPXiEHzKyWbb4VvpFdjwBVUd5Z9c+/rR+C4OjIympM1PvcBBuNFJZb8Ehy3x7/2LSIrqXB9Oh/TmUltYiyzJz5o5ClhXU6uZ3s1arg7/95Vvi4kK5ccE4gkMa/MaFELz7n02sWr6XXz5xDSNHJhIWHdzi4fM/f/c5a5fu5IqF47n/2RswGDs3irancMl7ufi4NNCoVExKTmJin0RsLhd1dveTVnmTkfGxAMwe1I8wo5FduQXo1GoKamrIrqym3uH18rvtZsTo5EavW1LmABUVdZhMFvJyKigrNGHQa7jpZ5OISwjjrgdnMmp4AiW5FdSaLASE+PPnRz/kyX/eifonWSbv+/2NzFowjrVLd/LXxz/ht/9e7HNh7ACX7ArdqTiRkND4cqFfNB8eP0D/kAiGRzaNhFx68DDHS8u5e/wY+oS2P9VqT8dst1NlsRIbFIhWfV6JnV2pOmW50fWeRG2tlQ/e20p+VjkHdmWdu37DTWN58PGGws211fXs23yCqPhQhoxLpbSgioiY4CYK/UJkl4xao25wA3XKreY+v9TojPS5PQKzy3zO5U8RCh/nfsyRmiNsLFvH53mf+KLMOsC1qQPoExzCmhOnqbM1XhlnV1bz6YHD/O83a0kvLe8iCT2HrChUW61utw/Q60kKDWmitM+aHXqqMgcICvLjllsnknWqBEmSuHxyX/oE6vCTBXZbw84jKNSfGfPHMGRcKgDRCWGoNWoObmsocr5x+V4+fmXNuTH3bTnB4ze8Sm11PZIk8adfvM+jc/9KxtH8zn+DPZhe/xNY7ajGqXESqgtFJam4LOIydlXtIlgTyFUxc3pldRNPYnE6MGqbd5cMMzTYTW8Y1rR259mgm30FRejaKObQ3XHKMhqVCmMzh6qKEOzKzSclLIyYoPO+/ma7nQB977UH79pxmoiIAG5ZOJ60/tGoFo6jKKccp11Gb2i5X3xqg7fW9PmjWbt0F7Iso1arGTN1IIHBRoLOuDE+9cZdCCF85pd20usVeqIxEUUorC1Zy+jQ0ST7J5Psn9zVYvUYnt+9mdzaaq5I6suSIaPPXS+1mNlZnMe8tOYLMYf4+RGo11Nnt3PtWx/w2NTJPDBpfGeJ7TGqLFZu+uBTbh81gqsG9G3ky64Iwb1LV7I1KwdocEn8aNFNqCSJ79JPc9PI9vub9xSc5XWEyDI33DkZtUZNnakep91JYU45A0a0HOvhZ2xYHKhUKubc3jghWv8L+kltpCL20TyXxM+fSlIxPmw8B6oP+Ews7cRks2LUaEkJDm10/YXdP7ApP4sau63ZfovHjuLj228CQBYCq9PVbLvuTpjRj88X3cq+gkIKaxpHK9ZYbczsl8rz115JQkgQw+NizgVk9URlvu77o/yw6XirbTauO4qiKFy76HIe//NtnP06BYb4M3b6oFaVObifRdHHxXHJHor6cI/SejMhBgN6dePNnCIETkUm01TF4PCWg57WnjzNwytWs2jMSJ65cka75y+qreXr4ydILysn2GAgyKDnhiGDSQoJ7lZ26HqHAz+ttlMjbIsKq4mOCW7VG6U9bP4hnW0/nmLCxL7MumII69YeYdaVQ88VubBaHTx395tcuXACqYPjWPfFbkZPHcCYqQNbXE3X1VjY+0M6M+a5V//UR/N4zA9dkiQDsAXQ02CiWSaE+L0kSe8B04CziTCWCCEOtjaWT6F3L07UlNA/KNrrSmh7Th6hfn4Mim47WGR3fgGFtbVYHU4+O3yUY6XN1/TsExLC3EEDeOzySd1ua251OvnZp19w97gxXDtogFfmWLp0J+++s5nIyCAeefQqxo9P88o8eXmVZGaUMmNmg2kt52QxOr2GwFB/VCoV/oHNG8zXLdvNlQvHI8sKm1cdYKabBa19NI8nFboE+AshzJIkaYEfgUeBB4DVQohl7grVnRS62WVGFjLB2vbl9+hNbCvLZHBwLKH67rENPlpSykcHDlFjs1Frs2Oy2XApCiNiY5g3eCDfnTxNqJ8fq9NPUGquJzE4mG/vWoS6mx2c1ZyRO9zovc/15Mli6ups9OkTTkiI/0UXo2iOnJxyYmKCMRh0fPftQYaPSCI+vn1BUNZ6O37NVDuy2xzoDY0P2YtyKsg7XUJ1eR0zbxyDRqvx2K6jt+CVSFFJkow0KPQHz/zXYxW6j56LcqYgrt0ln/Om8dGY8vI6Vq7Yy+Ill6Nrx2dUX2/n7y9/g9XuJEat4q4n56LRadDrPfM5O+yuJlkWhRD866llfPPhNiZeOZTqijpe+PhBjAGtuMtcYnjUD12SJLUkSQeBMmCdEGLXmVvPS5J0WJKkv0uS1KyPliRJ90mStFeSpL3l5T3fH7k74FJk3jy1lTdObWlXutjegkqSUKtU3UqZHywuZu3J021WI+oINpuTcjfTyK5csZecnIpWlXldjaXJNX9/PVdOH8TAhDAGjU7GP9DPY8ocaDZlriRJXHfn5UydOxKdQcsLHz/Y7OreR9u0d4UeAqwAfglUAiWADngDyBRC/KG1/t15hX6k5ggSEkODu793gsXlYHd5NivyDjIoJIYHBkzrapF6NUIINmVmM7NvaottZEXxuunH5ZJRFHEu1W1rfPH5bqZOG0B0dOsmRV+1oJ6BVyJFhRAm4AdgthCiWDRgB94Fep6T8QUMCx7mtjL/ofQQOypad+/yJkaNjsSAMIaExpIWeOlkpRNCsPzoMdadznR7Z7LudCYZlZUdmleSJGakpVBmNrPsyDFKas28vWcfK4+dfwY6w45vs9ipqWpctEV2ybiaKVJx083j21TmAC5X9y9m4sN92nwKJUmKPLMyR5IkP+AK4IQkSbFnrknAfOCoNwXtDpTbyxFCkGiMZGBgYpfKkhwQzj39LufKuPOBPRaXA7Ozeb/w3oAkSZSZ6/nLlh9Zc/K0W32+PXGSue9+1KK3THvmDtIbGBsfR2SAkZ+NHM6UlOQOjdle1Go1kTHnlbTLKfOP369AUdqvlA/vzgRwa7Xvo+fgzrIiFtgkSdJhYA8NNvTVwMeSJB0BjgARwP95T8zuQYQuAkmSSAuMI1TfNfU5zU4bOeYKym115JjPrzyfOrCSEmsNt299m6PVhQBYXQ4+ztrF1/mHukRWb/DAxPGsvXsx1w0e6Fb7huRYKk6Udfz8xqDVkBwWilqlwqDVetWLpTl+alfWaNU89vxCdBdh4x5+kW6On/17PYd3ZvgC9Lopbf48CyEOA6OauT7TKxJ1Y9zxd5aFgtpLZetkobCtLJOTNaU8Mngmn+fsJSWg4UdGr9ISojMyMSKFY6YiUgIieO7war4pOEK/wCjmJgzr9HJ63YHfTJvCo5dNQlZ8CqijKIrCd5/u5L0/f8u9/+96brxneleL5OMnXHrfcC9SYqnmo5wNOOWOF/RtDrWk4ur4ITwyeCZCCFYXHMZ2Zq5fD72KML0/vx02h1tSxrGnModgrR9PDbuGz6fd16IyF0JQYTOzOv8IFTazV+RuCZvLhXwR5gJ3qLZa+dOmLQDoNZpzHjGKEFRY2i4e3t1RFIUnF7/ZrP3ckxzZlYksN/yNVCoV9/zuOn7zys/oOzTBq/P6uDh8BjQPEmMMxajW8/ujH/LCiLu8OpckSXxw+fk5DGrtuesA02MGMD2m5ShFRQjWF53gw8xd7K3MBWDZ9HuJMAS02MeT2FwuFn+zDI1KxQtTr6JPsGdzpm/PzWu2fqgE5NfVEmH0XHHikro68mtqGBYdjcGDZe5aw1Jvp77OiqnKTIQbh5/uUl5kIjLu/N9i2ITGppnL54zw2Fw+PI8vl4sXsMtO9Oru4SNtl53oVJom5qItJae5f8cn516rkNh2zROEdHLU6Ob8bFZnnKTaZmV+v0GMjYknJqDt8wlZUfgu6xTXpg3o8tD/Vekn+N3a70kJDePrxZ1TK1SWFbLSi0hIiWzRZ9tuc2Kz2AkMMbrtmvj0z99i9q0TuOzqYZ4U10cH8RW46EK6izK3uBzcvf1Dfrd/JYpobNrQqTVcEXt+BR+oNXSJ3NMSU3hp+tWUWsw8vH4105e+zZuH9lDvbL1EW6nFjEat7hRlvi0nlw/2t5yp87pBA9l6/708PWu612U5i1qtot/QhFYDcHasP8ZT97zDqo92uD1uUt9IohO6V61TH+7jU+i9GKNGx38n/Qx/jY7an7gzToxMYXHfSSzsM4oBQdEsn3k/fpqu+SFSSRJfzLuVhQOGEOFnZOXpdN461PpOLi4giNkp/TpFPocs84cNm8iqqm6xTYifH+MSOteuXFFaQ1HeeU8nRVHIzSgFIDejlLdf/pbHX7yJK250PzHWXU9eR98h8a22URSF9AM5FyWzD+/iM7n46FYoQrQr++Oe4gLC/IykhXh3Vbk5KxuNSsVlyX28Ok9HeOK217l+0WSmXjMCWVaoqTQTFhXU1WL58ADumlx8h6I+uhXtUeaKENy7ZiVRRtVZ+foAAA/CSURBVH+u6zuQh0dP9JoJZlpqilfG9RRCCGbOG8XUaxoOLdVqlU+ZX4L4TC69BCEE+yvzulqMTqXaZqXe6SBAp2NDbhb1Tu+4i3Y2F7NrliSJa26d2K4+dTWWRlGmToeLytKaVnqcx2F3cWRXpi/AqJvhU+i9BEmSGBgc09VidCrhfkb+MmMO8/sNZsUNtxOga76YdXdFCEFlvaXJtZd+tZTf/fwtryvLFx75iCO7s8693rf1FP/34Hsc3N52WgW71c4f7nuHbz9x/8DVh/fxmVx6EYpQ+P2Br/nt8Dnn/NJ7O/P6Der0OU9WlpNlqmZOWv+L6i+E4OODh/n2xClMViv/v707j46qyhM4/v2lEiAJCZWQBRNISICwpQElQmRRBBqlaVAEwRWnmXHpnjONTffM6KCeaYeZblvneHRkRh2dVjmtLO0CLdDiAgiKCAjIIhBCWELIAiQhe0K480deSIWkKgmQVNXL73POO6m6b8nvl1dcbt13331r58+7tG7bFz+wac0e0icOblP3UdbB0yT0j8ERWPegi4LTRXy3JYOS4nLOF5Zx7HAucYlRPLpo+qXjhjtDePHJlbz52T/jcASQPmkI6ZMa5gbKP1VIfk4hqTc2nWWye48Qlu18ltpabaH7Em2h20TNxVqWZm5j5fHvOFic6+1wbC0uLJzbk698hI2IMDYxgV05pxkc0zBbZn5OEWesLo97fzHJ7f7r39/Bl2u/b1S27NUN/O5X71JWWknthVr+dvIfWP/+dtYt38b2jQfJ3H+KXn0iG/0n8cvFsxg2qh8f/t+Xzf6emPgIhqY1f+1ARHAEOpqd31x5j1boNhEgQnBgEI+l3MxQZ5y3w/Gaygs1vLTja97e9x1FlRXt8jvCunR123o2xlBSVeVx/01Hs3hi3Xp6BgcTF95w4fJcwXnysgt5+pUHSB50ndv9EwfE4IxqfEfv6ImDmff4bYSEdqWspJK4xJ7UXqjl5Q9+SdfgIMIiQnlt8WqKzjZM7xAaFszC5+YwY95Yt7/L2zdtqbbRCt0mHBLAnMSRJIdFEdgJJ+GqFxTg4OHhaXR1BPKbDX+loB3nbSmrqeb5bZuBujtXjTGICN+cOOl2H2MML23Zys5TOaRER7Fw/JhL61J+1JtDe0+Sd6qIwMueEXrx4kXWrfgWgIHDEhg2KrnRMZe/+gUORwAiQpgzhH9/62F+/dxcQsO68ez/zqe0uIKoXj3oHt50OoQrma1R+SYdh24jRdXl7CvMYUxMcqecWfFyi7/eQHRIKI+OaJ9nr2QVFbLvTB6BEsCL279iStIAfjN6nNvtK2su8GVWFk+t/5xfjx/L5P796BnaeKoFYwzF58pw9mzcAq+tvciDE37Hf3+0oMk6gOrqCx7nNv/8o52Ul1Yx/YExbrfxpKXjq/al49A7IWeXEMbF9vd2GD7jqTG3tstx66clSHJGkOSMYEv2MZxdg1l1+AfG90lkdFzzDz/Zl5fHkq3fsuL+uSQ6nXyacYQpKU374gMcdd0ceacKKTxTwqDhCZSXVlFeVkXh2dJmK/QuXQIxxrDy9Y3cPmcU4RGNJx+bdGfr7xZ1tebdrby35DNKiiuIiArjttk3MvcXE/WxdT5Kz4qypYqaGvYWXN3FYWMM7x7Y02T44OfHjvLBof2X3kd2C+Hb09lEh4TSt0eE2+Ol9Y5n1UP30zcigoKyMvbl5VNdW9ukz/35f1rBEz97g11fZxAbXzfz4Wcf7mTYjcnE941ye3wRYfJdaYQ5r90EayJwNu881ZU15GWfY+lLn1yzY6trTyt0ZUtBDgfHiouovHChzfv+af9u1mUexgBfHD9KUVXjeXBmDBjEiNiGi5YDI6P4x9HjmD1oKLGhrZt+OLeklNAuQQQGBLBk6zeXykWEhP6xhIV345sNB1n++iZys89x6/QR/P0zd7TY7REZHXZNL2TKZV13jy6aoRdKfZj2oSu/98Hh/QyMjGZoVEyj8gNn8skvL2Nwz2hiQ7tjjGF3fi7Xx7ofQQLwN2veZ9OJLP488z6CAwNZm3nYY994vQsXLxLYyq4IYwwZZ86SEh1FcUUFPYIbLlaWFFew/LUNpE8aQurIvq06Xns5V3Ce+RN/T1VFDdHXOSnIKeKj/f9B1256IbUjXbPpc0Wkm4h8KyJ7RGS/iPzWKk8SkW0ikiEiy0XEv27TU7axJvMQZyqajmZJdkaQHtfnUqv5q1MneHrzZ/zpgOdnrL41bRZvTJ3J5pNZDO4ZzYOpI1oVR2src6hriadE13WfuFbmAGE9gnno8SntXpnn5xRSe8HzE48io8N5ZNEMAJyRoWAMuSfPetxHeU+LLXSp+34VaowpFZEgYAuwAFgIfGCMWSYirwJ7jDH/4+lY2kJXV6r+c+r6dT+3rITo4FAcLhXpqZLzFFVVcvJ8EcsP7uXpMbeSfIUzMa7JPMS0fu6f+uTvNn68m27BXRrdHQpwPCOPE0fyGD91GABlJRW88+In5GefI+vgacZNHU6ffjHcNqd9Rg+ppq7ZKBdT9y+p/m6EIGsxwETgPqv8beBfAY8VulJX6o/7vmNvQS5nK8sJCnAwoU8SqzIO8kT6eG6Ijb80S+Ou/Bw+PHyAE+eLGRIVfVVzvE9NvrJb+/3FLdOGU1XZdEKzmDgn1yX0vPQ+NCyYnz9zJ9VVFxCBvy7bRlzfqEvj7pXvaNWwRRFxADuB/sASIBMoMsbUX3HKBpqdFV9EHgEeAUhISLjaeFUntf5YBidKiikoKyUmJJTPT2QyNWkAs1e/x4NDRvDs2MmICD/tN4if9hvktrK5aAzrjx3h5t6JhAR57iVsy1S+/khE6Bbc9G/g7ilI9bf5D7upH1vWfs/QtL5aofuYVnX6GWNqjTEjgN7AKKC5GZGa7bsxxrxujEkzxqRFR0c3t4lSLZo5YAg5pedJj+tDfPdwEsOcfJ+fy9j4BHp370FlbePRLO4qmj/u3cnCL9byQcaBjgjblhIH9OL+BVN0LLoPatONRcaYIhHZCKQDThEJtFrpvYGcdohPKQB+kjwQR0AA/Z2RbDyRxche8ZRWV7E5+zgPDB3htmtl1ZEfOFp0jpSIKCb0SeL57VsYGBFFas/YDs5AqfbXYoUuItFAjVWZBwOTgeeADcBsYBnwELCqPQNVnVtYl67MTkkFYERM3eRj+eWl/LjvABwevvZX1NTw0s6txHcPZ2rSAAIlgD0FuX43d7pSrdGaFvp1wNtWP3oAsMIY87GIHACWichiYBfwZjvGqVQTMSEt38Qze2Aq+eVl9ArtztHiQn6WegPlNdWsyzrMP0Tc1AFR2tfpE2ebTMmrvKvFTjBjzPfGmOuNMcOMManGmGet8qPGmFHGmP7GmLuNMZ7nDFXKCwIDArhn0I/4JCuD+z5ewSPD04gMCWV2ylBvh+b3jmfk8umft3s7DOVC7xRVtvbC9s3cPTCV2ODuBDkcjcasqytTW3sRh0P/jh3pmt0pqpQ/O36+iLl/Wcb1S5ew9MBub4fj9/Kyz/H0/DfYsemgt0NRzdAKXdnaf02azorp9zIteSCl1doreLVie0cy9+cTeW3xaqqr2j7xmWpfOh+6sr2EcCcvTJjq7TBsY9joftw2dxSBQdoe9DV6RpTfO1ZcyNIDu7wdRqchIsz+uwl6Y5EP0jOi/FpO6XnmrV3JjtxT3g7F9vJzCinIKfJ2GMoDrdCVX+sWGMjtSSmkRLh/ko+6NmLiIijILeKlRSupra3lDwvfJefYmUvrN6/dQ0lxuRcjVNqHrvxaZLcQ/iV9grfD6DQO7MjiZGY++TnFbFqzhzFTUomzHos3LL0/QUEOL0fYuWmFrpRqtVkPTyAg0EFEVHcmzxzJ6IkN8/T1iAz1sKfqCNrlopRqNRHhrvk3U1leTVlJJUVnmz4pSnmPttCVUm3m7Nmdp5bM83YY6jLaQldKKZvQCl0ppWxCK3SllLIJrdCVUsomtEJXSimb6ND50EWkADjeYb+wbaKAMy1u5V/smBNoXv7EjjlBx+eVaIyJbmmjDq3QfZmI7GjNBPL+xI45geblT+yYE/huXtrlopRSNqEVulJK2YRW6A1e93YA7cCOOYHm5U/smBP4aF7ah66UUjahLXSllLIJrdCVUsomOl2FLiJ3i8h+EbkoImku5T8WkZ0istf6OdFl3Uir/IiIvCwi4p3o3fOQV08R2SAipSLyymX7+G1e1ronrdgPichtLuW3W2VHROSJjo+69URkuIhstc7DX0Qk3GVds/n5AxEZISLfiMhuEdkhIqOscrE+a0dE5HsRucHbsbaFiCy3ctotIsdEZLfLOu+fL2NMp1qAwcBAYCOQ5lJ+PRBnvU4FTrms+xa4CRBgHTDV23m0Ia9QYBzwGPDKZfv4c15DgD1AVyAJyAQc1pIJJANdrG2GeDsPD/ltB26xXs8H/s1Tft6Otw15ra//PAE/ATa6vF5nfebSgW3ejvUqcvxP4BlfOl+droVujPnBGHOomfJdxpgc6+1+oJuIdBWR64BwY8xWU3fm3gHu7MCQW8VDXmXGmC1ApWu5v+cF3AEsM8ZUGWOygCPAKGs5Yow5aoypBpZZ2/qqgcCX1utPgVnWa3f5+QsD1H/b6AHU/9u6A3jH1PkGcFqfRb9ifZudA7xnFfnE+ep0FXorzQJ2GWOqgHgg22VdtlXm7/w9r3jgpMv7+vjdlfuqfcAM6/XdQB/rtb/lcbnHgedF5CTwAvCkVe7vedUbD+QZYzKs9z6Rly2fWCQinwG9mlm1yBizqoV9hwLPAVPqi5rZzCtjPa8mr+YO10yZP+XlLv7mGileHZvrKT/qulleFpFngNVAdf1uzWzvU2OMW8hrEvArY8z7IjIHeBOYjJ/n5fJ5vJeG1jn4SF62rNCNMZOvZD8R6Q18CMwzxmRaxdlAb5fNetPw9bFDXWlebvh7Xtk0tGahcfzuyr2iFflNARCRFGCaVeYpP5/gKS8ReQdYYL1dCbxhvfbrvABEJBC4CxjpUuwTeWmXi0VEnMAa4EljzFf15caY00CJiKRb/WbzgLa2hn2ODfJaDdxjXedIAgZQd5F3OzBARJJEpAtwj7WtTxKRGOtnAPAU8Kq1yl1+/iIHuMV6PRGo75pYDcyzRrukA8XWZ9GfTAYOGmNcuyx943x5+0pxRy/ATOr+N60C8oBPrPKngDJgt8sSY61Lo66vMxN4BesOW19a3OVlrTsGnANKrW2G2CSvRVbsh3AZoUPdSIrD1rpF3s6hhfwWWLEeBn7veg7c5ecPC3Ujq3ZSN/JjGzDSKhdgiZXXXlxGLvnLArwFPNZMudfPl976r5RSNqFdLkopZRNaoSullE1oha6UUjahFbpSStmEVuhKKWUTWqErpZRNaIWulFI28f9ZVJf9bc3EbAAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<Figure size 432x288 with 1 Axes>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"plt.figure(figsize=(10, 6))\n", | |
"ax = gdf.plot(column='hru_id')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": {}, | |
"outputs": [], | |
"source": [] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": {}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python (camels)", | |
"language": "python", | |
"name": "camels" | |
}, | |
"language_info": { | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 3 | |
}, | |
"file_extension": ".py", | |
"mimetype": "text/x-python", | |
"name": "python", | |
"nbconvert_exporter": "python", | |
"pygments_lexer": "ipython3", | |
"version": "3.6.4" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 2 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment