{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Dfsu - Export to shapefile\n", "\n", "1. Read selected item and timestep from dfsu\n", "2. Extract geometry\n", "3. Create GeoPandas dataframe\n", "4. Save to ESRI shapefile" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import geopandas as gpd\n", "import mikeio" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Step 1. read the selected data" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiIAAAFMCAYAAAD2otedAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAA//dJREFUeJzsnXe8JVWVtp9dVSfd3IkO0HQ33WRJIrREQRjCCIggiCKggxFMIIKoCI4EQUQUHVA/x+wgjoJpRMFBBFSgB1BQMg3ddI43nVRVe31/VDiVzrn3No2NUC+/Q99btWvvXVXn1n5rrXetpUREyJEjR44cOXLk2AIwtvQEcuTIkSNHjhyvXOREJEeOHDly5MixxZATkRw5cuTIkSPHFkNORHLkyJEjR44cWww5EcmRI0eOHDlybDHkRCRHjhw5cuTIscWQE5EcOXLkyJEjxxZDTkRy5MiRI0eOHFsMORHJkSNHjhw5cmwx5EQkR44cOXLkyLHFkBORf0Lcf//9fOADH2DXXXelu7ubbbfdlpNPPpknnngis/2jjz7KUUcdRU9PD5MnT+a0005jzZo1qXaXXXYZxx13HNOnT0cpxSWXXJLZ39y5c1FKZX623377cZ3DeMcCWLZsGSeffDIDAwP09fXxxje+kWeeeWZc4zz22GOcf/757LnnnvT29jJz5kze8IY3sGjRohflvP74xz9y4IEH0tXVxYwZM/jQhz7EyMhIql2j0eCCCy5g1qxZVCoVFi5cyG233TauMQJ885vfZOedd6ZcLrP99ttz3XXXZbZ7IddvonMd7/nnyJEjRwjJ8U+HE088UWbMmCEf/OAH5Rvf+IZ89rOflenTp0t3d7c8/PDDsbZLly6VqVOnyvz58+VLX/qSXHbZZTJp0iTZY489pNFoxNoCMmPGDDnyyCMFkIsvvjhz/Jtvvlm+973vxT6XXnqpAHLWWWeN6xzGO9bw8LBsv/32stVWW8mVV14p11xzjcyePVu22WYbWbt27ZjjfPSjH5WBgQE588wz5Wtf+5pcddVVMn/+fDFNU2677bbNel4PPviglMtl2WuvveT666+XT37yk1IqleSoo45KtT3llFPEsiw577zz5Gtf+5rst99+YlmW3HXXXWOOIyJyww03CCAnnniifP3rX5fTTjtNAPnc5z4Xa/dCr99E5jqR88+RI0eOADkR+SfEPffckyIRTzzxhJRKJTn11FNj29///vdLpVKR5557Ltx22223CSBf+9rXYm0XL14sIiJr1qzpSA6y8NnPflYAueeee8bVfrxjXXnllQLIfffdF2579NFHxTRNufDCC8ccZ9GiRTI8PBzbtnbtWpk2bZoccMABYx4/kfM6+uijZebMmTI4OBhu+8Y3viGA/OY3vwm33XvvvQLI5z//+XBbrVaT+fPny3777TfmONVqVaZMmSJveMMbYttPPfVU6e7ulvXr14fbXuj1m8hcx3v+OXLkyBFFTkReRnj1q18tr371q2PbttpqKznppJNSbXfYYQc57LDDMvvZFCKy8847y7x58yY03/GMtc8++8g+++yT2n7EEUfI/PnzJzxegBNOOEEmT548Zrvxntfg4KBYliUf+9jHYtsbjYb09PTImWeeGW772Mc+JqZpxhZsEZHLL79cAFmyZEm4bc2aNfLoo4/K6OhouO1Xv/qVAPKrX/0qdvwf//hHAeR73/teuG0i1++5556TRx99NLZtvHOdyPnnyJEjRxS5RuRlAhFh1apVTJ06Ndy2bNkyVq9ezWte85pU+3333ZcHH3xws4z94IMP8uijj/K2t71ts/QXQGvNX//617bzf/rppxkeHt6kvleuXBm7VlmYyHk9/PDDOI6TmmuxWGTPPfeMXesHH3yQHXbYgb6+vljbfffdF4CHHnoo3PaVr3yFnXfemfvuuy92PJAaa++998YwjHD/RK/f6aefzs4775y6BuOZ60TOP0eOHDmiyInIywQ/+MEPWLZsGW95y1vCbStWrABg5syZqfYzZ85k/fr1NBqNzTI2wKmnnvqC+4oimF+7+QMsX758wv3edddd/OlPf4pdqyxM5LzGutbRea5YseIFndOKFSswTZOtttoqtr1YLDJlypTw+M1x/cY714mcf44cOXJEkRORlwEee+wxzj77bPbbbz/OOOOMcHutVgOgVCqljimXy7E2mwqtNTfeeCN77bVX6m36heLFmP/q1at529vexrx58zj//PPbtpvoeY011+g8a7XauM/pkksuQUQ45JBDYscXi8XMeUTHmuj1+/3vf4+IpM5rPMdP5Pxz5MiRIwprS08gxwvDypUrecMb3kB/fz///d//jWma4b5KpQKQafWo1+uxNpuKO++8k2XLlnHOOedkzi2K/v7+CY03kfmPZ6zR0VGOOeYYhoeHufvuu+np6Wk7dqfz2pS5RudSqVRe0D2pVCo0m83MfdGxNsf9H+9cJ3L+OXLkyBFFbhH5J8bg4CBHH300Gzdu5NZbb2XWrFmx/YGZPDCbR7FixQomT56c+QY7EfzgBz/AMAze+ta3pvbNnDkz9vnRj340ob6D+bWbPxCe81hjNZtNTjjhBP7617/ys5/9jFe96lWbfF5ZGOtaR+/NzJkzx3VOncZyXZfVq1fHtjebTdatWxceP5Hr12ms8V7/6PZk27HGyZEjxysXuUXknxT1ep1jjz2WJ554gttvv51ddtkl1Wbrrbdm2rRpmcm77rvvPvbcc88XNIdGo8FPfvITDjnkkMyFJpn0atddd51Q/4ZhsNtuu2XO/95772W77bajt7d3zLG01px++un87ne/46abbuJ1r3tdx3HHOq8svOpVr8KyLBYtWsTJJ58cbm82mzz00EOxbXvuuSd33HEHQ0NDMRHovffeG+7vhGD/okWL+Nd//ddw+6JFi9Bah/sncv06jTWeuU7k/HPkyJEjhi0ctZNjE+A4jhx33HFiWVYqhDOJ973vfVKpVGIhobfffrsAcv3112ceM97w3Z/+9KcCyDe/+c0Jn8N4x/rc5z4ngNx///3htscee0xM05QLLrhgXGOcddZZmXlT2mFTz+uoo46SmTNnytDQULjt//2//yeA/PrXvw63/fnPf07l5qjX67JgwQJZuHBhrM+s8N1qtSqTJ0+WY445Jtb27W9/u3R1dcm6devCbRO5flnhuxOZ63jPP0eOHDmiUCIJdVqOlzw+8pGP8KUvfYljjz02803z7W9/e/jz0qVL2WuvvRgYGODDH/4wIyMjfP7zn2ebbbbh/vvvj7lmvve97/Hcc89RrVa54oorOPTQQ3n9618PwGmnncacOXNi47z5zW/ml7/8JatWraK/v39C5zDesYaHh9lrr70YHh7mvPPOo1AocM011+C6Lg899BDTpk3rOM61117LOeecw3777cdZZ52V2v+mN72J7u7uzXJeDzzwAPvvvz+77LIL73nPe3j++ef5whe+wMEHH8xvfvObWNuTTz6Zm2++mXPOOYcFCxbwne98h/vuu4/f/e53HHzwwWG7Sy65hM985jPccccdMcHqf/zHf3D22Wfz5je/mSOPPJK77rqL7373u1x22WV84hOfCNtN5Podcsgh3HnnnSnB6njnOpHzz5EjR44QW5gI5dgEvO51rxOg7SeJRx55RI444gjp6uqSgYEBOfXUU2XlypUT6veOO+6ItR0cHJRyuSwnnHDCZj+H5FhLly6VN7/5zdLX1yc9PT1yzDHHyJNPPjmucc4444yO1yrI8Lq5zuuuu+6S/fffX8rlskybNk3OPvvsmIUgQK1Wk/POO09mzJghpVJJ9tlnH7n11ltT7S6++OLMayIi8vWvf1123HFHKRaLMn/+fPniF78oWutUu/Fev+CebOpcJ3L+OXLkyBEgt4jkyJEjR44cObYY8qiZHDly5MiRI8cWQ05EcuTIkSNHjhxbDDkRyZEjR44cOXJsMeREJEeOHDly5NiM+MMf/sCxxx7LrFmzUEpxyy23xPZfcskl7LTTTnR3dzNp0iQOP/zwMDdPO7iuy0UXXcS8efOoVCrMnz+fz372s6kot0cffZTjjjuO/v5+uru72WeffViyZEm4v16vc/bZZzNlyhR6eno48cQTWbVq1WY7901BTkRy5MiRI0eOzYjR0VH22GMPvvrVr2bu32GHHfjKV77Cww8/zN13383cuXM54ogjWLNmTds+r7zySq6//nq+8pWv8Oijj3LllVdy1VVXcd1114Vtnn76aQ488EB22mknfv/73/PXv/6Viy66KKwNBXDOOefwi1/8gh//+MfceeedLF++nBNOOGHznfwmII+ayZEjR44cOV4kKKW4+eabOf7449u2GRoaor+/n9tvv53DDjsss80xxxzD9OnT+eY3vxluO/HEE6lUKnz/+98H4JRTTqFQKPC9730vs4/BwUGmTZvGD3/4Q9785jcDXtHUnXfemT/96U+89rWv3cSzfGF42ad411qzfPlyent7UUpt6enkyJHjZQIRYXh4mFmzZmEYuXF5S6Ner7ctBrk5ICKpNaRUKr3gel3NZpOvf/3r9Pf3s8cee7Rtt//++/P1r3+dJ554gh122IG//OUv3H333VxzzTWAt9b96le/4vzzz+fII4/kwQcfZN68eVx44YUhCfq///s/bNvm8MMPD/vdaaed2HbbbbcoEXnZJzRbunRpx4RW+Sf/5J/880I+S5cu3dKPuVc8arWaTNvKeFHvc09PT2rbWGUwREQAufnmm1Pbf/GLX0h3d7copWTWrFly3333dezHdV254IILRCkllmWJUkouv/zycP+KFSsEkK6uLrnmmmvkwQcflCuuuEKUUvL73/9eRER+8IMfSLFYTPW9zz77yPnnnz/mubxYeNlbRIKiXldffQ2zZ8+lu7s3FPc0mjZWuYAC72sF2I5DQ1y6KqWwnQDD9Tq9XeWYMKjh2JSsQmy80WYDyzAomma4zXZd6q5LT7EYbhMRao5NpVCIbdtQr9NfKmIohfa3jzQbWIaibBXCiYpA3bWpWPFbWHWadCXmtKExQl+xjKVab22jTh3L0HSZLd+hFo2t65TNCMMXqEuVilEOrwXAqD1Il1nCMgrhNsetosSmaEVSpougZRRLdcXm5OhhTKMLResNo+msxzQqmKo1J1fXEKlSMCeH8xEELUOYqs/f5M1Au0MooweFEc7UcdejKGAYPZHzrOPqISyzld5cBFw9hGG0CrsBuO4Qyuj2+wzmvhHEwAzbClpsbL2BorlV2KEAth7CMvpQ/hgADb2GgupH0bpPjowgOFhqIHLpHRrOGsrWjNicbHej12fk7azhrMYyejFVKbhMiDTRYmMZXd41EhBxcdxVFK2ZsT6bznKKxlagWt8nrUdQykKpUmRbFa2HsYLzDM7f3YhhDkQ34TrLMM2tUKrgbxK0ux6lrNZ1FhCpIXojpjkzvJciAu7zGOZsiHxHkCFQXSCtvy/Rg2hxwYhfO0ePoFTrfooIVWc5mNuG91MQ1jWHUPRjGaWgYezajIrNiYccx9zps2Pbh4aGmD179piFA3O8+Gg2m6xZrfnDfVvR07P5Ld8jI8LB+65m6dKlseKPL8Qacuihh/LQQw+xdu1avvGNb3DyySdz7733stVWW2W2v+mmm/jBD37AD3/4Q3bddVceeughPvKRjzBr1izOOOMMtPZWjDe+8Y2cc845gFeU8o9//CM33HDDmMU+tyRe9kQkeFg/8shfeeKJx7bwbHJsOha/CH2mS9a/cKx+ge2WZ2zbkLEta+7rxzk2wNMZ256bwPFZ92NjxrZlbY7POqes8ZdkbMs6FmBkHNsUsDSj3cY2fXr4b/OX/Pu/XZC5L3f5vnTQ06Po6X0x3GTeIt/X1xcjIi8E3d3dLFiwgAULFvDa176W7bffnm9+85tceOGFme0/9rGP8fGPf5xTTjkFgN12243nnnuOK664gjPOOIOpU6diWVaqEvvOO+/M3XffDcCMGTNoNpts3LiRgYGBsM2qVauYMSP+wvOPxMueiCSxauUUAMRUNPvip68NcMvglIm8hCnvZwWu5b+IKUAJYrSaBJDg9+izSeO1TT2vxItbMmJbPJiRtzIlqIIm63lnGBqV+LuzTBdU8GLnTcYyHMoFp9XG0KDAVJreQiN2uiXlYCrtt3MpGk64XwDHP5mCcui3at5hvllJASYuA2YVpUAhlIzWuAAGGpPI+Qg0sBAUJppuI+7nDfouoaO3I3Y9FfG3WAA30sCI7Hdbp5q6JwLoxEbdao0OR0q2EUD5xyq8d+3EF0FF+8vGuNY05fUcve0GLoZKXwOFhO2U0phKYsdFr5uJYEb6UMF8IttaX+PE2fn3Jdq3qYhZvIKxTBRm8KVVgutvL2JgBB2Fc1IUVPoxpZR3hQMLCiI0aH3PRrQb9htAgLW6wjKnP9zmYvDQ8Gw22l3+GREb36hC32+FEmVy5HgxobWm0Wi03V+tVlNaJNM0Q0tIsVhkn3324fHHH4+1eeKJJ8IionvvvTeFQoHf/e53nHjiiQA8/vjjLFmyhP32229zns6E8IohIsViH67bZHioB1DY3SZN2wyf/Nr0iYhAoxJfDQISoU1wK942bUlm8LMgYMY2oLTyHpkWaTKiBCyJtUcBBYl0ABgao5gmI4ahMSOkxTA0BctFBERaEzSVZqCrFvZZMHX4s20qBoq1sG9XCQXTW667TM+FAx6xscWk4C+lBpr+0mjYT7hwIZjKYZo5jFJQUc3YIlnACRc8EahLAR1eTGGyUc0mXQhlXJ/gtF+0RcBJXGjDn58Atu+6yeKGWrzFKeyLFvkI5ihCa9Hy4Yo3KS2ES79uEx2vSR+/KTCJkgYXK4OIGOFIgoVGqezzBvz90WPAiDWU2Fcb0pzbjC7igBGjIh6KGBgqmJXgRIiIlWTVKMoUUlaHGAnxUdc2GkFEKIqTaO/dm5XNMt3SetjXdIG+co111W6yrooYwRi51SPH+DEyMsJTTz0V/r548WIeeughJk+ezJQpU7jssss47rjjmDlzJmvXruWrX/0qy5Yt46STTgqPOeyww3jTm97EBz7wAQCOPfZYLrvsMrbddlt23XVXHnzwQa655hr+7d/+LTzmYx/7GG95y1s4+OCDOfTQQ7n11lv5xS9+we9//3sA+vv7OfPMMzn33HOZPHkyfX19fPCDH2S//fbbckJVXkFEZNttt2PJkr9TKjdp1Eu4ZdUiIQZo/0oYNt5KEz4oW1Cut0GUZD+X0utAuE2hwJWWRSW6X8axTRuIDaoYf58WCZZXf/6G9klIeqEMTiu+R1FzC5i2pq/oPaBtMSnjoJDQMiICjpixBVRj4IoKSUVrrVA0pMAGt5tJ5iiOMiiE9gkJrRMi0JBCwgLhWRXMjIupUTQwKIkOjA4pMiISt4S0h0JnkJHWXY/f/6SVJDYm+FYDFZKdrK9CdIxO+zc/JCQZ7UZXSHgtNapF3DKucbzn4Hj/90j75Fc42Bb9eke/zbrNH5BnW1KJLQlrh4CpDLS4IbFJDtxE0W3YjLgtk2dTTPqsWsZMM04sR45xYtGiRRx66KHh7+eeey4AZ5xxBjfccAOPPfYY3/nOd1i7di1Tpkxhn3324a677mLXXXcNj3n66adZu3Zt+Pt1113HRRddxFlnncXq1auZNWsW733ve/n0pz8dtnnTm97EDTfcwBVXXMGHPvQhdtxxR37yk59w4IEHhm2++MUvYhgGJ554Io1GgyOPPJL/+I//eDEvx5h42ecRCeKzf/e7P/CjH32fdWsHWL9+gOqsIiiFKHCLxJ5D9SkqJCYCMQuHUwZdECSuB40928WMEBUNKkIKRGWQEVMTMQh4Owv+YovE2irLRVkSWxwsyw0tEQXLIf2e6qG/XMMyNYZvok82GSjW6LJsALrNBmXDpttq+iTEiLgdWphWGKLisbdMK0W/McqAWfXbeFaUovKsLU1tYpO8GNBj1Ckpl2wIBYSi0om3dQ+uZJOGwEWhASdmqUiTEUdUSLhcAotIxEqSsIgk3Tmuf3zKNdOm/aYiarlQSXeXDzOyvdVWUrYaI+WWaV2XoM+s41rtW5/ofSkkzlMAC4WlDBRgR6iIQlFRSZsLWJiheyaLhAT9ahEaYlMTl5rECXtdPIK70S3zrDM13L7G7qEpFg8MzsGR+Ngi0BiymHVLk9n77Myn3vWR2P7g2TI4OLjZdAM5Ng3BvXjg79NfFI3IyLDm1busyu/1i4RXTPD7HnvsgohBpauOLnnWEAF0gdRaYTSj+gxiJEM5IPF1LP2CGZN3JOwPolBJkYBkLEpthATiGOCqOPHxj2+98WabaxztTTxslsDGZpmGa3ptxaBgeGTARfluhvT7ra3NsM+sN+dB3cWILoUvlIa/jNhiYJPlqwJXOn0tFTYGtqRfUnUbEhKfcbq/5LbQYhN+EvcwaYWJ/izBV6bdfXhxjPxZdg7lW5+MlBUs6/ikniJ9XTr1EbV0JLcn4YZusuSYkkpXHbQP97exJSlaDzM70ofnTmx9L7p8QhzsczBRCvqsWur8lIJRu+i3fVm/r+XIsUXxiiEihUKB3t5pVCp1dNl7COlCNgcwHOKrXNTv7TL2U10S/yagtPJetUO7v4o8xQO7dvgamjwabRve9PzjPa2SYBi6w+QUbkBEMrQEQZv1jS5sbeJoA0u5uKI6EANFc8wIcMU6t5dRXcLTGGhc6Xxc8s00C03MuFl/DJdMQCayiIqgfN1GMGNvq2e5SFtsgqOyEFimxvqSZIlrJ4rWfH33Wcrp1nKDZR0Xb5uGRkWu8diqlpCMSHxbcr8msCpljZnVr+9uHOOaKaVQQuiaCcayIzMv4IbfHCdiU+q16qn+m65JU79ivNc5cmwxvGKICMBOO+2MUlDubiAWbSJZfJ2IyiYCSmhZNNo9F9u9Hkb70Z33e/vaEwbdMMPdgUUkK2oiCs8i0rmNoFhX76IpJlpaETLt4D2ox1pUFcudSTjiPfgbkmGGis5zXF9LRVNa4tF2jpwosqwb0X1hzxGtRPvRo8cmCUCyx87Hbyok8v9wm7T2WhmumvixY6NThM94+o3aMGLuL3SmJiRrm0LBGN9tCPQp8RC0ZsIypRR0q5YWKkCfVSd5V4abpdZ7QW4RyZHjRcMriogceuhCALrLtUyXTADlAG3e2ATBrPtt2pCJdJxAdhvl0nrSO4kGosZYLRRu0yMjImCZnawh/hDa8CQnY0xPoxhqlHAyw3wSfWKgx9WnwbBTpjYGCQHfQpFlqkqNrXB1cOnGbp+1cLegWm/q4q1749FxZFoX1HgtHptjcYsHyLZ+kg73JE2cOrUN3E3jRZS8RG1D0bN1Jdu+kUVETDEYyxoC3kBmNMkbLa1PFF1GExDsiFWuYjYxIjN3taLqtBIIDtvVscfPkSPHJuEVRUTmz5+HKxa95mjHdgrfPRP8EtunUK7CbCisusKsgVVTmHUwmqBsPIIRJSoacD3yomxQTTAaCqtuYI4G+w3PRTMRaIU7alIy3TGJQDD7YLHtBBGFUool1ckM2yX0GITIHmeC3paGZWyMZRXRAiO65Ifijq9P7zQ6tVVoBVWxqLd1yYxjnHEs3GrMuYw5CgW8L5qJxvAltcGtMjKX4OjR0bl0/kKMRXGjnsikV7Jdz0FwWvK7mNaNgCjtEZeE9iNLw2Lhi2A7EMmKsgFFU7csIkp57pmg7xG7iGnoMFS9r9Kd1VWOHDk2A15RDlClFK7ZT0Wvo0SThhTbPmGNpuAWVfqVziUWMaPw2gQi1KC5oZWXH2QMhmBgIA2NWKAaCl2UVpROu0M1YCufuCgawwpUHauo4/PNgKtNlNXe4C7iWU6qTpGy5bLe7mGj3UWPVae3UM/IVeG9WVZSJp14nwDLnQGmMcwkw3u7bH9pPLFsMSNyxhFFVUpUfcuKdoV+o0nR95e9kCSXAmitwjwjnW6fUp0JnQ6XZUX7m9Jpuz9Om5lauJgKDPHEqJ6lSxA0GjMz/LkdxiIiYdr1yJySRCP0VorCQajgGfWyzl7Eaz8oQkVByZ+B921utRS/8QZxaEgDE0W/MjGVl53EQAUp0MI8I0opTAzsDufUZXiRYHYiK0pfocag45UiqDmWl1PF9frpLnWl+smRI8fmwRa1iFxyySWewCzy2WmnncL9hxxySGr/+973vhc05rtP97LJzahtwBpxUW4otIi1i1lE/DXFrHtiVbN98jsCraLn3hnfqqj8/0BhNFW24MEnQTQMaPjZ1/Ce9No2qQ2VqdUKOG5cyJrsJIicSe3x305rtkXNLjLSLOHqQOBpMORUWFabxNpGd+xN0nuzbM9ntUBdFxh0u6jqMkvtqTzZmMGwrqD93CZZVyQqkBXx8j1scCus1T1UpRVv3aDAkC6xQZepieUlFGuzBrXVh/jHOGJQI3AdGaHwsp2wMks+KRLk4fBycQSWChVKP9vbCpTf3sR7E/fSvEmYfDf4WEiYOMxUcdeKAoqpdG5ZM4+OO3bb4BoEX0MN4ZkFEUIiyjf+GdSE8F5EiYxHACDInVsTGPWywCF4CdCCdlqEDdqh4YfiuggNvEysDpqmuNTFpiY2NbfJsNtgvdtkRCS811n3raACB1D8zPsiFhHxtSaG47Usl/LMqjlyvFjY4haRXXfdldtvvz383UoUcXv3u9/Nv//7v4e/d3W9sDeT1y7ci+9+p8xAeZgNzw54S0RR4VQUTpeB64f2Gg1BuSCGZ+kw6j4nMXwXiwOZHglpkRixQY+nJlL05RmFaoIUBEr+dleBozoTGwHtGjS0gVJCoeBiWR6jiSYac9w0EdGiqDsWdcfyzNuG9yAftov0Fxsxl8qoW2LULVMybPoKNSqGTSPjQriiqOuCHx0Tn3eNIuvcHoZ0hS6jQZ9Rw1K6dQnwIxoE6mIxKiWcVE7PyDlhYuBSkwI1sSgrl7KfvbWTgSggLLYYYQRFfOFS4RKqImQkuuhH2wckJEuDEe8xmFew9CbbxNvHF02JaRmSUKgxRcuRGY+znWdhGM/DoqXjVjREKKiWgc/Vgq1UimfbwLAI3agw1LiBZlinU5NVxaWISmVabaJpoGlIK3MuxC0ysesa6ThIwtZlNJhUGGXELlIwXGxtei8UQFclJyI5crxY2OJExLKsjsV2urq6NmsxHsMw2Hnn3fn73++jq6tGtdqF2RTMplAa1IgBTgmsqoMullCulw4+XHT8B6VZBycrK7RuaQSUlz1r7KucWC0VCmUr9DCR9O8Zy2n06RrRl4gomk2LZtOkUHApFFqExIlYM1ytqDkFmm6gh4gntaraBfqKzcTy7FsidIE1jQKWcukvVJlhbcQ0wdYmdSmMGW1T1UX6zDqjusyoLlFWNr1GnbJhewu6Mliju/3w2c7QeOJWw89uWheLupgUlUtFuTGS42rBNLz2jihfi6LCfrLRIiTRf6JoR0LaQeEt0GmdQ9xuE7zXB+2SScs23ROlwrMajztL0SJu7dq7kRmJP4ItXnh5yYBRVNu7qfHISBGhIS6jkh0HJUBNNF2RxGeuCA2fAnm2srRFJHr/FNBr1ljndoeZfbUoRtyS55pRMLU8yrp6V2gRqZQrHa5Ojhw5Xgi2uFj1ySefZNasWWy33XaceuqpLFkSr7b5gx/8gKlTp/KqV72KCy+8kGq1s3q90WgwNDQU+yRx0knHIAJ9fcOpfcoRyuttClXNVv9Xo7JGwodRrJ144tTYSqJBufHFwU862hHt1gEjSF42ruVGIYlEZ6CwbYtqtUizafnmbkXTNhhulBhslGm6LYtF1MQfHF9zrI5aCEcM1jV7WN7sZ9CpMKLLfh6QznN2MLEjVQPrUmSN28dKu5+atjy9xri/nipD3OrlKhnURQbdYjhWVYo0tEldTN/KEl2wxrIQ+C4I8EObA5fExEhIHDGbSmYPwbYsm82mIehHiApcxz6qw7627iuoYbBemzTFHFMovVZrhtuQkAB1PPGqN65Qj2ZnVVAITYxJBPdPMbewnt0qy+g1a2GyvqiL0VAwtVylpL25dJdyIpIjx4uFLUpEFi5cyLe//W1uvfVWrr/+ehYvXsxBBx3E8LBHEN72trfx/e9/nzvuuIMLL7yQ733ve7z97W/v2OcVV1xBf39/+Jk9e3aqzaxZ0+nqnkZPdxXTiDz0XI1Zt1F+mIi2FFgGpaHsp6fRjJh4fZdMciFRoV6kw6Tb7dNgNpLkIuiY1GDitE885jgm1WqBes1i+YaBiBWkNQkzw6Tf0Gab5dE/ShS2NnmmupWnaJjAOlnTxdS2uhR4ojmLdW5Pql5OJwRvtWl4rpshXWKVU2FQKj4pUanjx0MktEBDLKoUGRWP4DSk8zVqBxX5v/dT+4ungsiYTTeBEBAPT2sinntDqQ6Onta8Wo6hbA1O1BoSjBQgsDR5+pLOJ+AK40kZQtUnKw2vzF1sXyvXb3soBcO6xIBZI/gDS+bEUQp6/Zwj3ZVcrJojx4uFLeqaOfroo8Ofd999dxYuXMicOXO46aabOPPMM3nPe94T7t9tt92YOXMmhx12GE8//TTz58/P7PPCCy8MCwyBV4Mgi4wcdOAB/Pa3t9DbO8LGwX6U42I0Io9JBW7Ze6RZVVB9ghTSvn+j7lXkVVlBI/4zzbDBneiV9nWNAigbJL1mZ8CziqhorZvodMQIs6s2mxalkhM5sp3JXVF1LHoKdtZOBLC1gU2JZ0anMKM8RJdpj1ksDVpWkUIk572XgVWx1umlp1ifQPSHZxUptkttJlCniGBQlSI9qhndRfsIFg9avCiLwJUTWEFcTFD4sSovxFrR+disDKkT7b+9syzqzMjam+wJP0rH+93TmqaPlvB/LdhiYJIdbh5kx22f+TfSD0JduzgZbZUCS3THEHCPQFqYStNj1BnRFWwxIw4rv6/g67TFbcc5crx88ZL68xoYGGCHHXaIlU+OYuFCLyFZu/0ApVKJvr6+2CcLxxzzekQM+vqGMZo2ZiPxribglFvvVm2tIq5HFDKfnapluDCaGfsjY6V+d4NjPb1I7LW1w8t7e6sIROuA1epFPzW8N2CWNSRAU5temfsMNCKWlUG7wqpmP8sb/TR0ITyVTqjpYviG3dQmTT82WmOwxukd4+g42llFRKCBRVCyrUEh4hbqbA3RAnUxqVGIuXKSx+hNdppE3sA7tnshppB0sboMadOERw4iijo5UlJllZC2MT22KFzaRVKlUUW3zXjaSVkk4llDAhvUJLMKCA3XSv8p+hrVoeGR8U0qR44cE8ZLioiMjIzw9NNPM3PmzMz9Dz30EEDb/RNBqVRi5swFlEo2JSs7HjewiIBvFXEzHOECRsu62xbKJV3sLgsBCYn0p1AY7Vw06ZEytCLp+Yl4ZMQ7YizrhacVScLRKkZQGtrC0YqGLrCi2c+KRt+YtTocPK2GCIzoeGTCsK5Q04UJuHuytCLe23o9EdkzKqWIriFLfdoiIG4iuZlnDUkv7e4L+HMaO5dH+h6Nh/rEquhuwgyytgpepAv4JZM6uNDSlhKFLWkXXqtWkBpfRtux/t4UbS1USrW+awooGQ5lZWOnSmrDSNHbtnLd6jHnlCNHjk3DFiUi5513HnfeeSfPPvssf/zjH3nTm96EaZq89a1v5emnn+azn/0s//d//8ezzz7Lz3/+c04//XQOPvhgdt99980y/pFHvh6AvslpAawuKMJkDfgPrMG0y8NoQnnQIyqd1oXAKmKOgjUCRoPwVTJIiha2zXjFDAvlZa0oNl6CsyA1fcIqIgI6JXr1ImscR7WxhsS3NbUZ5hUJ+my46dDcqtPyIdV1keWNflY2esPcI1kLSFUXqepixgKkWO2kLVoiHunZ4HSxuDGNxY0pDLqe9iPLKpKVVt7FoOFVY0nt86qyqhQBCeAdk7VAj623yMZ4mNb4LQWt2bRzx2R31E6KFN3v4oXHagwaKJp+5FHW3DzBcRapkFT6u6huJ1oPMnOe/v2xJZteBcdabUSrrig/IV6rvz6jmsiP42VXHSx48ffrNqzvMKMcOXK8EGxRjcjzzz/PW9/6VtatW8e0adM48MAD+fOf/8y0adOo1+vcfvvtXHvttYyOjjJ79mxOPPFEPvWpT2228ffbb0++/90uuidVWb+8P/Zm51TSxt3CKNRtQSyvnTUKpv9qaFXBGUPPprTvxsHTlFiAGKAtcP0cXcrB91O3ebKrxM/RHFm28kiJodAiqJIvbhT8xsk+FbZtUinGmY9IK54i1HqIouoW6PWzUjZDEWvMo07NLdAjzUguC0VNl1jWKDLJGqXHqlMxbAxF2LeLQcM1MIz0otKUAhvdLvqMKjUp+aSl5FskvIug0GzUXWzU3RRw6TZqDJg1isrBESOVQTNATSwKPuuL6hkcDOq6QMFI04q4NSRjkcPIfKdvF/oahtBmzjDRR3jM2C3HesPIsoAF9zxLchokKoufceTvBQNL4mLaVjRRanQcURT8PC9a4rlFXFFe5EumzgmaGAyJhYWmK4u1R87RkETtG/FKA8REwgpKygnPVgRGnQKDzRLKD5YZ3DjYdpwcOXK8MGxRi8iNN97I8uXLaTQaPP/889x4442hCHX27NnceeedrFu3jnq9zpNPPslVV13VVvOxKVBKsecer8E0ha7+Wmxf1C0TRWlQMBwoDrZICHgkwxrLRSMtl0vwGFQazKZPQJq0d9YLPnPBc900Fcao/3HC3ryPVlC3kMECetRE2tSwEYGyn2NEBGzHoNooMFQrMVovMtKssLbazWC9TM2xaDhm6I5pZiRGC+AVC0tCsd7uYUV9gGfrU1ne8MJ9bf8ttEF7F8w6t4c1Tj+rnX5GdCVCBJT//9absY3JRt3N082teKIxg6puZWFNYtQtsd7p9jO3WoxImXW6hw26p62bZazoGEcbMeuPK4qmmDSwaGCFVpvsd/XOEMYupyh46fH1GEUT22efbSGgerbv8hprvknLSOdrpcL/Jy0bLu1dhXWfhHgUTjHaNjuvh6RVJIiWSaIaieCqORaDzbI3hqFwSorqUK4RyTF+/OEPf+DYY49l1qxZKKW45ZZbwn22bXPBBRew22670d3dzaxZszj99NNZvnx5xz5d1+Wiiy5i3rx5VCoV5s+fz2c/+9lEDSbh05/+NDNnzqRSqXD44Yfz5JNPpvr61a9+xcKFC6lUKkyaNInjjz9+c536JuElpRHZEnjzW96ACPRG3DMCuKUsIuLdcGskQ5yqfPdMB5iN9m+zpY1Q3ACGnd3CsKGw0cAcVJhVA6OpUNpPDd+mCjAosDsU03MV1dESGzdW2DhcptosYvvi01b6dUXTtRhplllX62bJUD/PDg3EXDDJMatOMbU4BIuyYXht6rrIeqeHpY3JLKlPZnWzl0E3O3ulYPjRNNlIyWHE02s0KTAs2fkfmtqkKmU26C7Wur0M6S7q4kXVKPy3/6xz6PAn44liCzQohOTDxgrzVASLp4uXyVWjcP11snNhwfgYY7bBswI1MTP77aQJimYq0XhumImEJgfVmL2+xp6sK94nti3pKvP/HdUmIxLNeyPURLX/+uPlA4neMS2ePijs2/9ejvjfvaaj2OCTkAB2RdEcjr+o5MjRCaOjo+yxxx589atfTe2rVqs88MADXHTRRTzwwAP89Kc/5fHHH+e4447r2OeVV17J9ddfz1e+8hUeffRRrrzySq666iquu+66sM1VV13Fl7/8ZW644Qbuvfdeuru7OfLII6nX62Gbn/zkJ5x22mm8853v5C9/+Qv33HMPb3vb2zbfyW8Ctnhm1S2NyZMHGOibxUaWYxYcXNtClwzvCRZAAPESm4lptn0sK/GsIk6FuCdEPCJhNKHdC7rS3hPTqoGdlF64/n5AGRkOFlEoV8jOI6a82jTleKY18bUkdcfPI2+5FCsOptlyAHhVeOOuF9uPhhmsm5QtJ1NfIkDNLdBlxUN+G66V2gbgiLdgbrC76TGbobskOteqFCjhZUkdC25EGbHR7aKimrF5ihAhNgobM6yyGqApJl2JaznWstrwF8nAjZGlXmh51zzLksYrEqdQFMUZUzTcKczW4zRGxHGkaGJiiRdaPK4cJOL3o8Z28LQjRS4GiG57fHCYK7DRLVGlQEE0ReVQQGMqjzAYgVFHYEgsGjEXm4TXcUQb9JvtvxcWQtPX2IxKyXcFKbRvnQJFv1Vj565lPDw8k4IhiLho//7YFYPyyiZaawzjFf/ulmMcOProo2PpKaLo7+/ntttui237yle+wr777suSJUvYdtttM4/74x//yBvf+Ebe8IY3ADB37lz+67/+i/vuuw/wrCHXXnstn/rUp3jjG98IwHe/+12mT5/OLbfcwimnnILjOHz4wx/m85//PGeeeWbY9y677PKCz/mFIP+rAo448vVe8qJJVc8aknDLKEdTGNFY9Q6vrv5DPmUVEU98atY7+/ZDeYokiuoJmE6HI/2VTbmq/QBa+RlaI4fZURmj9zC265YvavXg6vZmbwGG6u2sFIpRO24V8VwS2W/ogeVFgHXNdLl18XUmnawi0b6iugRBMejGxTs1P613OLeESDHQgdiJuItOmV5diUfNtCvml75JiiBn6riiRaLahtj2lmEsqeNwMGgShGB3JiReP2rcFpp2XzqXuHA35OTifRequsByt5uqn5TdxmRUimyUCht0meEg2kpgUAoJEuKNG4zcRNGQ6D2Pw1DgaljvdvFccxINKXjkN0JYC4bL1OII08utytCmEgqGxqkolJCZpTnHKwvJrN2NRocKqBPA4OAgSikGBgbattl///353e9+xxNPPAHAX/7yF+6+++6Q8CxevJiVK1dy+OGHh8f09/ezcOFC/vSnPwHwwAMPsGzZMgzDYK+99mLmzJkcffTRPPLII5vlPDYVr3iLCMChr9+fH//4R/RMrrJxdW+YPwQRzLpg2K3Hv9UQnEr7J3loFemK/F4dnxgR/GWpgVfwzvAjaMbQnYT/atpTy6YBlrc0iCbtrtEKTGjWLYoVG8MArRWq7TfEc8H06gbFjLdRjUHdtahYXnyE7Wd0c8XEiFg84qnBFaO6RNUtUDRacRWBzsMWi6Y2KRjt83NGa8eEfUqJHqlRUQ6uKD+ja1xoGU3AFlSTrekCRdOba7ukXcE5xIv7eYt/IUu0GjkmSWwcTIwxrSIk5upB++PGt7augaf1MBERilHrkH9UrLquKKwOvHY8cMXws7f65AMDW8wwRDeYV3Ke3rl4xMkRlyFt4ar0lzpZ+G9UK4pGmmQF59egQMWwsZTGkYyxfStZX7FO0XC8sHO/yVClxGRsvv+bX/GBt5w68YuR4x+G347uRLn9Q2uTUR91gFWp5JgXX3wxl1xyyQvru17nggsu4K1vfWtHDeTHP/5xhoaG2GmnnTBNE9d1ueyyyzj1VO87uXLlSgCmT58eO2769OnhvmeeeQaASy65hGuuuYa5c+fyhS98gUMOOYQnnniCyZMnv6Bz2VTkFhHANE22n78LhZJLqbuJLiqUFgqjGjNCQkThWUUSr7vJNArWqGc5MRuCNZp4c82yCCRelBWegBVNOmNrYqzAU6FQnlWkHWkRoOktuGJn+HDC4xR2veAnPxsrZFQx3GxvpQi0Io42wkXX0Rm+pQTW2L3+wu9nbZXW2/CoLrU9vj1ZUKx3ej3TvBuPmAj2R/OPBN07mLh+4rNO1hAHI2O/yogyIRw76xTGZxXxXT+iwmsUJyFjHe2P5Rv3gpo5jj9fV7xzCerCjNswEum3KSY2JjYGNW0wLEXqUoi5zMbCsJRY4vZQJUv4nO7FRVH1dU1aB5YXRUMMRnWrvGCvUU/1BR6xG9Vebplp5bgwdem0Pqpli4fuvJtV69aOa/45Xp5YunQpg4OD4efCCy98Qf3Zts3JJ5+MiHD99dd3bHvTTTfxgx/8gB/+8Ic88MADfOc73+Hqq6/mO9/5zrjH034Wy09+8pOceOKJ7L333nzrW99CKcWPf/zjF3QuLwQ5EfFxwknHADAwaxjL1VgjOp2ATIHZiLyOZkGE8kZN33Oa0gZJiVrH/ZapPV1JxzgJie/vEMlIIFwVW6WZE0Ak+kBE0ax7USyu7kBuIEwZnwVbmzRcM5afwZX2i3MARwzW2T0gUE2EWrqYNCS7EF+7vB8ANhYb3S7sNkZAOxLpErXQ1KQQLvhZaFlDMvahMmqriN9f9nULkrul+/IIh4OB7YtghZb1Jtm2HWwxsCPEQ4diVF9xoXzXjK8vkTGiUsIxfULT8OoDA4F7K51rZqw5tqDaFspLJlpzUWwUizXaZASLQSkwKgXqYtGMENleo5E5uisGTbFQCrYqj4RtRKDuFHls3hRMV3PNt785rpnneHkimbW7VBrbXdwOAQl57rnnuO2228aMCP3Yxz7Gxz/+cU455RR22203TjvtNM455xyuuOIKgLBK/apVq2LHrVq1KtwXJAONakJKpRLbbbddquDsPxI5EfGx3XbbYhVm0NXVYLuBpRStjJzsCgwNht2GXYhQHNIYTShUwaqN/bjNzMmk/dTxnYq+CSmipFBjFAxTUG9ntkwsZuJpRlzXaMueRKDetHhqzRQajgpTfgeouxbPjw4wEtGLOGLEc01k9OtogydGpvNsdXJmZtZqhlakk+skwKDuaruoRt0z0X5q2mKFPcB6pyfzuOgiF0XLY2Z4acvDPZ3nmEyKJniLrIOJixmSBhFPtjqRJGeei8TyrTzt3UzBPLUomv78O8EVj7Qk3WKB9icLE6lWnJUHJnqkiwqtLXZohcp2+1hKU1FxwbQAIxFLWdFw6S94lhPbMdFisGZyN2smVdj4xDP86eGHxjXvHDnaISAhTz75JLfffjtTpkwZ85hqtZoSS5umGVo55s2bx4wZM/jd734X7h8aGuLee+9lv/32A2DvvfemVCrx+OOPx+by7LPPMmfOnM1xapuEnIhE8JXrLmGffQ+haNnMmbGM7vJoZjvPPZPYKEJhRGIhuoYDVjWxOndaOLQnarVqXvK0rFVGMn6JPXI7uWc6ji8pS4nWimbditSkiQyjFY2GRa3phbyuGell1C4yYhcZapQYapRoup4gcGOziw2Nip+DJBIWm/XmL17WVFss1tj9rG72sbbRQ9Vt5RnRGNQkLoZ1O1QXkTY/x6Fw/LmJP4+mNhnWXTQoMiolmtqMhDX7VoC2VpjWtmTekLEIU0BGnDDMNyttvVe0rZ0rJ+Pr6ecCUTFXV9bYSY2Ng0kzcm2C/jxrUCDqTROLTvFN4w8JzraKBEd7EUpRAW+6Z0mQo5R7RogJoaPumXqjQHDWj82bilbw7e9/P3z458iRhZGRER566KGwLMnixYt56KGHWLJkCbZt8+Y3v5lFixbxgx/8ANd1WblyJStXrqTZbL0AH3bYYXzlK18Jfz/22GO57LLL+NWvfsWzzz7LzTffzDXXXMOb3vQmAJRSfOQjH+HSSy/l5z//OQ8//DCnn346s2bNCvOE9PX18b73vY+LL76Y3/72tzz++OO8//3vB+Ckk076x1ycDORi1QiUUrzrzLeyw/bb8f3vf4etp61k7eBk1g8NIH7JdACzoaE3vjhYVaGQYQGxRgVtKHQp7dIJrSGulxzNiLhWlKswGqCzJA1Bm4wdyoU2ngIPZrulWPkxk/Ft2jWojRbp6mmGwsOmbeK6fkFA3wRTs4tYDYmE/4JpuGHYbENbrKl301No0leohlk1Y24XbTDiFKm5RZTyLCpdZpMmFk3bYhChYjbpNpvUKFI2m57QkrGEpF4xPUcUrqHoMpoUDTcm+gSwxcRSgitGJHtr61ps0N2UIuphA/FK4RljmSW8N3YL7Y/ZfhHWAo4yO1bybREC04uUEe39FBGYJnVJDq18IK5v7TAyzGftiE2Q0r0gGlOJZ6URhc4Qkrbm2J4kTUR7IrTEvxCoZFquKmJEJKrCSVuMlPKIyGq3N2wZ6EOibSYXq1jKpd4oElTkrVYKPLt1P9s9P8hXbrmJdxz+rxM4ixyvJCxatIhDDz00/D2oCH/GGWdwySWX8POf/xyAPffcM3bcHXfcwSGHHALA008/zdq1LU3Sddddx0UXXcRZZ53F6tWrmTVrFu9973v59Kc/HbY5//zzGR0d5T3veQ8bN27kwAMP5NZbb6VcbuVo+vznP49lWZx22mnUajUWLlzI//7v/zJp0qTNfRnGDSXtyle+TDA0NER/fz+Dg4MTysq6dOlyrvjcNbjOMMPVblYPTcb2H1YCjMwww1wjZkNTHMnwOxcB5T3Emn0KXVS4hQhRED+/SBtth1uAZj+xVcUtCJge4QhcN1KQmDFDFzSxl95oRnIlUJQ25EbA0pkSmGLJxiq62Amhq1KtSAXLcOkuN8PfTcPNWKQFS2m2Kg/TX6wBioa2qDpFGr4bRtEiCP2FGsXYBfIWGEu5zChuZHJhNNY+IBe2Nn0CEtcoVFTTL4jmUjEalHw1sHe8t8I3SNemCcbuUY3U9bFwY1E+AcyI7ywIKy7ghCWMQstKZOlsYlFSXhq0dnylpTERLzLHn7uJ9smFhPahJAmJzrlkpPO1xBd2IZ11X0ALdYp4Adntw4FtSY/raV1UR/FvFhRCd3jvPBtTOkIKgswg2reNtCNEy+yBMLGZIwaP12fG+hKB50Yncf8z83ADjZMSSobNIX94HjVnOld94LxNerbk2PwInvOfu/91lHtehKiZEYeP73Nnfq9fJOSumTaYPXsW13zhUgYmb0tv1yjzZyxlztSlTOreiGU4WA3/rbihKWSQECBcHRVQHPISokVCMrDq7UkI4OlA/EJ2ye3Rbam3PjdjZQiO6WhRbv+m7mqVIiHJIxxt4kRSv2dXZfVcIMuqAwzbRdY0etjQ7KahW4t/dGGrpdLFezsdMXm+MQktKmyvRVHXBQadCiO6TDNV7K61aLqYjOgu1rs91KSIq+H55mSGdZlO1yFLr+Bg0tDZAloI0rxb2OIJbRtiURfLT/te8Cw+/kcwqEu7rLXhWYTXIqo9cTFxxMD173E7EuLN2UjlCslyy6TGFkXdj2TRHfQj7SJuFNnEYCwEVpHgtywSEp+xhMdlzS1wzwjZkVRKQZ/R8EmIYFguVskJXyJ6S9lZgHPkyDFx5ESkA8rlMp+7/BO8/bR3sdWM+ZQLNtP717Jg+rPMt55nijNIaaSNMVulfy2tczFrLmjxdCZjQSWSmwWQhFsmOYFkzusYcVEdypumdSKpCSXaJ7upNRM1YzLHUbjaYF2jDzehV0iO0BQzVU231bXBertC1S2wtDHAymYfNV1s+7adbeMwqOoyq50BbCn4hKgdPELRLmInvk9wfZ2JLVZkkVexn1utAyuHN06jXSROJnloQYtBDc8N1Tk9u8KJaC/auUuS47mJeWsMP1Fauv9MItNmnPEg0Iq0DwOO0rJMiZW3T0GPUfdoXEIfEsVQo4xSGqvkYPiWQuV3apjtdTY5cuSYGHKNyBhQSnHQgftw0IH7UK/X+c1td3HXXffAxhX0qSq1/iLPj8yg4cYfZrFnoAiFEZfioENhyECKJigYnm2ldCOxw5TnSVEORKqWpyNjPG9Faz9+yvd2d9dRYLV5X03pRDy4joGRkbgsOX0tioZjhsX0staLIBy25hQoW3GXRtblqLdJDQ/wfHMqPW4jdE/UdZEes55ZOTcZ8gmtxaruX+DxCEmjeoUoPDICRRyUUn5F4XR+kaylOEUofMuRGXHRtLQwKtYu0JOI4IXcoqiLooDbMdLcxqDg170N3Bkdz108N0b6vL2U7tE0+smrI7HtE7eIeMcaftXe7P1pW07iDyMCQ0G3ajBCOYOICJZyeXrDVMxS/BoGXyvTzN/hcuTYXMj/miaAcrnMG4/9F66+6hKuvOpK+mbsSNlsMr9/KVPLG4guMMFLorI1ldVNShsdEDBdhVUXzJpQWdPeLyMQPkPNBvHXu/EYU1zf8pFVEUz8fZkDZ7t1xB1vqKii3iy0TRHu5SbxQy214eUpiQ6U1Z+brswbXZRbURWeKHTQ7WLYKaXmkEVEHDFiZeHtNnkromhKdqVgLR6h8eqZQJOM4n8Z/bVzYSTdKlkUJlq6vhlGr4Azrj/tlrujUyr/AHaHsF9PAEvYV/Z5esdvWu1hD80O707R+9v6fmRDgF6zjuPnDwGhoBwqRoM+s0a32aRUsFNELrCImEZuEcmRY3MhJyKbiEkDA3z+M+dy1gc+jFmoMKN7LfP6llEw/Dd3BYUhh66VTcym/8ZqqJhupLxBUxjKZgQq8uKrBHC8onm4CbdMFgSUrbCqJqrZpm1Qvya1HqiUa8X73UA7Y2dFDVBvFtJNBd/N0ppT3bVSbZJwMTLfxAPYGblGGmKxwemmHtFuhBJMT29J1S34eozWfDw3SmdoWnVlRDwNSENMmhR8q0iBjbob29d/dELLJZPt3rAjCc7a5eXQCRISHOtkZnaNo6HNsP5LuzkEY3QWmKpQd5LMxdKCV9yvNdLEyIhCQgvOmJBQrpr5MRBK2Aw6Fbp88tFletFUAfmYMTDCtJ7h+Bz8KeeumRw5Nh9yIvICsefuu/Claz7Httu9iu5CjQUDSxgoDVIYtCkOOnGppJl+OHevdDHrGbVaEoaC4ggUh738IlYNlE22iNX1XDlBRI3hKrDxLCPBxyc2uP6/yU6iBcQii5924snN2pv9FU3HRCfq2Qgk9B6KulOIE4U2fdbclm8q+bYrqDAVe7RvwSvvPuhUsHWrkrAtBqO6lJl7RGil+u6Ehlje23SY5TRe8C4Q6jalkMrbMZHlNyrSrGmLmi7QFE+UGlgg2mlBbMy2gwVulhpFhnSZESll5osBj4R0yj3SgjdXV9LaIRX5P3gPHp/KjdGnRxqCj9lZbe2PIvQbLlONJj2GnfnpVjaL6vOoUaJgpCOUgjwyU3qqTO5u5RMy/Itu5kQkR47NhpyIbAaUSiU+ecEHeevb/w2lTLbpWc2k0nB6acggIkqgZ5mDclN+BA/iJUYz/EXdEIVywWqANeolQPOzOqEchdJpkaDhKpREPlqhxEDVTLBVWswaibqIhzQqJEIuOttlPBdNFG5GGllXDK/+zBgLUkNbHSvCtstwCl5Uy6DbxYhTpOoWaCSsIKl5S2HM8FKNkRmeGiA61boUQouORHaO5T4I5y8GNV3w+vFT3NekyKiUGdRlv4hf9rlkWUU8N47lW2sUDSxcTOpYmW6VwLYwXow3NFeBnw4u7apRIfFoyV4DvUw2SfRmaCD0KXfMon1KwbxC+7oxTuTeTusZYaBS9cbwB7esXF6XI8fmQk5ENiMOOWghV1x+KVqbTJ26ERURTAq+ayYDyoHu5U7sCSv+2mw4CXGqiv+ogv2ZL4oSqBDbrvNKG9BQYSSNaF8P0uYY1wlcN+m33iRst0UOolaCJAJrx1hVZxvaart4u2K0ja7xoBjWFeodo2JabUfaRFIE6JywK7lwKy9E2J+fCrnf2EXgvPwjJiO6TOvPteXYEBRNrDaRK37m18iNckVRoxC6lqJoUqDqk6ZoNd5OBDALBaXbJGTL7ijqqmlZP9JXxvTFp0KajAREpU85YaK2sb5P2xfXMKcNGWlE3H1KwfS+YfrK9ZZGJLeI5Mix2ZATkc2MqVMGOOSQf8GyXAamDIXbs9wyARReZtbyupb/W8ggIQkInhdF2WTWpQkFs6h0Fd9kQ8eAZiBiVYhtIK7hEZPIHLTT+sqMWa4+IqR0dPtFt+GM5+3SF6122N/UYwlNPYLhdCAswZybEi/W5+3zzmVUFxnW5WwVhLQ4YXwuipqU0B20Lum5ePPxLB5pqa2/3BIQnexieX4VYPHOqdGmCF0A7RMbLwTYCLnseGGhMZRX02UiTqiWdaSdUkUiVyB9LYq49ClP3xE9vtK5+BLzC6uZbm1MzTVZ40gpmNk3yIwhz01jp32aOXLk2ETkRORFwCmnHIdpdjMweRjLD0/tRETAe3hW1mkKI9ojITpNQlJaRc8W7UXItNcZehjPuiAqUqtGeUIV1wDHQNyWq0ZnJUzL6s4XutZti7pt0XQMHFeFZdpb5zW+/hwxfUKTjfHqGEbcclsLApHFLqoV8USpBkO6HOYFaZd/IrCIpHQHKKp+RE3ncxbfBaGo6mKkbVQXE1hUPGja5x5p+AQkKwlY0ikSkBYJCYlFs0Mdn/jZSSxtfKc09cljjTEq00TDmD20LHMVHHqVCxl/BgWEtjYw5YXxzimsY6o5TPRKeBaRyLmMCFN+7zD3byO4RYPXv3a/cZ5bjhw5xkJORF4EmKbJKW85GcMQJm+10XvYj0FEAlRW2J7iP+MZnqxaG1o8xkEyvNwibXZKvF1mX9onJI7CbZip41JdCriuwbrhboaqniXA0SZN16LuFKj55KThmGiXjFTiadQdi5W1XobtEqNOgZpj0XBNbG3gaK9+ij0Oi4MnYk2H9yZhi+npQARqusCILoWLvwKaYqUyqgZWinZda4JImLFdMjVdjJGNzki7aMQXmTYotq24m55nfO4tR8l4ZuBFxYD3XU1aRdpdlcCt0ukLlSVSVQhd2PQYaUtItO8u1f4MlE/6tiuuYbI5Gs6h6Vpef1ro/rvLtF/YlJcLM3ebw9WXfY79d9iz7Vxz5MgxMeSKqxcJB7/utfzi578GVrJxQ4Oq6h77IBEMV+hdYjM8L+Nt23/SenqT1u8KUsXuRGXY1F2i4QrtoZVnjsn0PSjEMdC2gVHSwbRBvIgY0cqPllEoQ2M73pulSDK5lp93QhQFK13zJFzQxWunBUadAgZC2XTDBGbJSY4oi+nlUawxCtFpDAadLvrNKvGcYyqsV6MxWO30YaEpRAr4RdtWdYmC6YaLbFSomSyqF6BTfZYAo24JSRSUixOE7GW3KkV6aPiExwq3RxOfxfsMqENrn8ZIWShadqJsGG3OyUTapoEPejZ8a4clbmY143jfHjUqqCCBmkKk8/VUCroRhjPIn4q0mV9chduYwaDuoqEtrHVC/70uxfWC0Vfg1FNP5cA9c0tIjhybGzkReRHx3ve9k6s+fwVTp29kyaoKJCuVBiEUSoGjMeoOyjJxSyprjQ3N6Fkvt8rFK3TX8YGskKzMqVnig44LpcJteiJIL+Aly7MfLZjmkROzTeXfml1g8foB+ioNv5tOGg5PcGmp7Ek2dYENzS4qpk3RDLKGRt/KWz1pvKiSYsbbtpYgi6jXxhDBjJiUVNiLZ13x9COKLjPI8hpckzR5GXQrkTf4pKxVcMSLGamQTKgVF6lmQQuMShEr9VVLz0USP0fPKUmiPMLSTmyqY0QhnK3yNCNZ0VLxs2qJSw3R6IQrKLCWmHgFE+MZXBUuYz/ITAVdCNUIGVHguZL87gwFO5RW8VhjJsZfTKY+4KAE5u65gHP/7UOUSp0FzDly5Ng05ETkRcSC7ecyZ9udWbL0UXq7Rhmu9oRPd2W7GNUmRs3GqNkg4MzsB3yRqi1IMeHPV9kkBHz3TFI/kmrkGTokuYZnEZF2ZMRvq20TLYZXg6MNwYge77oGppn2DQXRNK6fkGs8DixHG1htqgUKBsNOhapbZFKxlvmmLBGSMcmqphddwbcmtDY6YlIkTWxEYKNbwcbCRNAYdJsNDCVIgi54dVJUKMzMcuO4ElhVFA5uZir5dtACTb9gXh/xIkXZ5AICspQ2nhlYsbGDecbbGintRhqW0jgd9DuxFOoI8SunKSiXotIx/Un8zMbBnYGiAkeEJnEyksTswnr+8uyC0D367ENP8ZELz2Vg68lsPXsbXrfP69ht3k5jjJYjR47xIteIvMg4+wPvxHUV0yatxxqqUlw+SPnpdZQXr6e4agRrqIGyNe7UHlq+FkVxGIwGMZLQSVagxCcj4YbgIGKhvW01ICTaZK1/GceJzp6USrhG3ERESyBklYSmYWwIdpvomGSUSrtsrK0FGDbaFVxfBxIgyz2gQ3Fka48rimFdxokkNLOxGHS7vAge4ufm1Y/xCFy7+dsRs1Zd4qnyg/6yasIEFX5pY7nQZItnW9cjuY/U+bbISFAxZmwSolQgWm1vTUm3dwGXsmowYNQpG24bEuKH/qq06LYdKir+0AuyrIh4SepGpQRKUX7TMBuOUYy8VlHfXtEsOqx7cjUP/+4Bvnz9l8cxUo4cOcaLnIi8yBgY6Ofg1x5CoeAytWcD5kgT5cYfvrqvjBS8Qnie60VhuoJp+wnLgifsWGkn/KRoEjy0AxKSeF0MNYQCYTaQ5FO8DcEIXQ7K98tLFoGQlPhURKUEmhIx2SuVrafIgqA65AxpbY/mMUmPG8Bgg9PlL+I+YWhTTdbxk7uIeKnRPfGqahED1ZrfSCRnSRCFEpCMLFIAcRISdFiP1bVJ524RAVsb2JK04CT1JYlrn3GGyfNNEh6Fpoj2aZf2060nHTzZsEKWHG+TlYTeUEK3simMEXoLEE3IPx4iohR0R/pVCK4oRqUUEjkRGCjVcKcY1Hc0GNnfYONxJuveZtCcAWrYRbdLQ5sjR44JIyci/wCc+s6TKNBF39wmViVBQgombn+l5bIBLyTXwRevglX1eMFY0alhrhAhLHinUKm77GVY9ceKiwNakIxPjM1ETyLLV5TcINiOZwnQWiEpH5OMkZAs3leWtSNeIE4yC9gFlonkZIecMsNOkeWNAUadEg1thYm9Ajhignh5PYIaNZ2sOFVdYK3dQ1UXaOgWyZDEcYHuJZhXdJ+LGStgF812G1hQksnJFB4RiV+j9JhjOTMCK5BX48WhpDzLhKGgoFwsJb7rKK4LSSKwiqiUAyXtEzTRFMeoGhxFUoA7HsuI6ZOR4NyqEkRDtSxOFdPGSoaZWQq3X6E0DA4Ojm+COXLkGBM5EfkHwDRNTjvj7SgDJu1YD7cLeC6ZxFNXlKcjCF4ilUDJd9XEPs3ExwGzKl6CM/8/IMOS0hLopVQZ/pNcifIK4/kJzlKLlmr9kHTPBJaSJFzXjFlBknB0tssiDZWydgQ5JZLt0paB2ORjbYecLrzaLia2WNR1kVFdouYWaGoz/DgRMWWn6TZ0AY2iKmVGpUxDTN/Vksi8GnPJpOfWctGo0AXk6UGsUE8S7SygtEkyEiVg483dohBKOCnXiMIXlyrxxb7jsF6oOAlPRsMU0BSUjoXiuj4ZyupdkQ77TlObyM8iaBEcaalQXBQVZceOccVAKRgoVFNuLtcPfntu1fKO55rjlY0//OEPHHvsscyaNQulFLfcckts/09/+lOOOOIIpkyZglKKhx56aMw+DznkEJRSqc8b3vCGsM073vGO1P6jjjoq3P/ss89y5plnMm/ePCqVCvPnz+fiiy+m2WxurlPfJORE5B+EhfvvzVb921KZ5lKe6pku3P4KkgxvwBOkCh6xCLcBYrWWKoWvCwleAQ0wUBhaYcQWHMlY24KaNJHOs6DJeMUNBo1uir9tG0ZyTH/5V5KRP6O1GHVKVpaemsrIA5JeXO1EdtR2VhctLXKQjFTxKupaVHWZaiJFfLs07Vo8a0ZQjE/wNByjUgoL1wXTT7tk0ufluWi86xroQdqF8EYzkEbJSCdNSDsEbqS0DalF6ZQSz4rRoecsrUjQ3kBTop0OROEEZCSx28oUMsWtIoInULVFsPFqPAZHBYRtK2OYEi0yEoQbDxSqqWtcLXsuvMWrlrU91xw5RkdH2WOPPfjqV7/adv+BBx7IlVdeOe4+f/rTn7JixYrw88gjj2CaJieddFKs3VFHHRVr91//9V/hvsceewytNV/72tf429/+xhe/+EVuuOEGPvGJT2zaiW4m5FEz/0B86Lz38MlPXMSkHessH+5D95Uzbdli+BYRt2XK7vj+GuUUQf244NBOa7vfRqHCRS7Z7bgzkmsFZmBOSQwj+Lk6Oofxgqf9SOfryJ68rU1KfiROuzd8T9hqh8XSstq1omhav2eRNwBbLCw/vXeWNiaA45MLR8zIgun14WAyrCsMa6FPjSLj8EMELhoR6RiB0hql9VtwbkFOvfG4ZTyKJXSR/aakkj8rKIimY/IzBYaIb5nx1DUW2g/F7gQvRNckGbEjfpba4KNa//rfZ1OEUpspRbU9W5nDrHD7cSLksdeqo2hF8IhAtVwAbJavWTnGnHO8HDE0NBT7vVQqZYZ1H3300Rx99NFt+znttNMAz0IxXkyePDn2+4033khXV1eKiJRKJWbMmJHZx1FHHRWzkGy33XY8/vjjXH/99Vx99dXjnsvmRk5E/oGYttU09t3rAO7/69307Chs2NhmMfAfnIavE0Epr6aMTutEJCAbwSJjtLQfoqQzEdEQehkyFl9psz1rDRNXgSEe4Ygck1ysOxMRzypiqPFoBBSONiiZ7pj1ZRwxKCjtLygZ1ouEq8QVhdkmSZZHAnwiklryW8cHic004S1MQUQxKF2UDZdCRBeRdclFhBEpUVH2GFlo21gVxMBCYaksMhZdbj1YvhhV2vjZkuejEIrKpYBnUXBRfo6PQPuifF4snqENoTSu+9wawcW7lqZ/dZ1YXp7EVVOeFaWKoijZ0T2B5cNz8QgzzCFWuP3hvTOV0GvVGXIq4TkPWUWgyrMrc9fMSxF3rtuBQr1dNepNhz3aBO5k9uzZse0XX3wxl1xyyWYfbzz45je/ySmnnEJ3dzxZ5u9//3u22morJk2axOtf/3ouvfRSpkyZ0rafwcHBFMn5RyMnIv9gvOO9b+XP776XSdOG2TA4QDImVwLnO77ZOyAfCgzbK/1CZMEi+Xvwkhc5rh2UeAtQeHDqWZ5hKWnbn/KMIZGMrCIgGqKFSrVE81mkF3tHG2EisrHgSpZ7Jg1bmxQM3dYtk9SROJKd5My7RF6V3/TS7bcREhYLA1eMjFopQX8mDTFwMCgTkIz44q9FURev5kwBTVG5bcmNIitM19tTlwKWaFTm6t8ar4gTWikCF5iRiDSJHmf5+o6AIJlIzHoRCIRdFKZyqUsBI9J+/FBo2nHrqDvSG7Pp/3HUBLoS+hTPOha3OJpopptDrHW6Cf4YBgpVhpwyoHC1YrTgxenUhkcmOvkcLwMsXbqUvr6+8PctleTuvvvu45FHHuGb3/xmbPtRRx3FCSecwLx583j66af5xCc+wdFHH82f/vSnzIrRTz31FNddd90WtYZATkT+4bAsi+1mz+e5VY9RsGxsO8HeE9KNaMZUwwG33NqXTHehIi/8hu4QgRuFbxVp556ZiKBAEknVQj1IzPQ+lvul8/7k9BytUllEk3C0idZ2xqm0Sz8eWBBa5KElmBRsMSi2KdzjZgSkuhgpoSZE08ELLiZV8chIIUI0HDH8YnbeRJpiobWXgt0IFvzItXIFCirurgjOqahcNEZmqvcSLhoo4fh5OVoqEAcoRCxEgQPERCgimULW2O++PsTAK0DXpVyG2xToGw8aYtHVqZy0QB0zvLp1DIoiMRdQlntKKSiIGxOv9lu1sN2oXQBlYhcNivW8+u4rEX19fTEisqXwzW9+k91224199903tv2UU04Jf95tt93YfffdmT9/Pr///e857LDDYm2XLVvGUUcdxUknncS73/3uf8i822GLilUvueSSlMJ3p53SGQtFhKOPPjpTffzPiAXbzwOgSD3lu4gWx1OA4UTeRKNpGzLuXFTrKYzhlgn6jK6PcS1h68VcIg06WVhU6y1dawWiUonNwn1kv9HDxKJn3HEIXAXlh7imB0wXg5OQAGSN57mEzMykYkgrY2vyuKQ1JiA2wdt7MM+aeBE6AE1t0pBCqj+NgYNFUwrUpEhVF6lrP7JHTJpi+O6mYC7e8Z5bJOu8PNphKe9j4tVyKaAJzjRwMYFnOSijKbfNdpqNwNhjt3GRjQc1bdEUq+13R3wSkkycP5II5W6XBUQpGHbL4Xe+ZDqUDI+YVP2XBqeksKv1Nj3kyPHiYnR0lBtvvJEzzzxzzLbbbbcdU6dO5amnnoptX758OYceeij7778/X//611+sqY4bWzxqZtddd40pfO++++5Um2uvvbaNOfmfE9vvvB0A5eYwxmgjRkaS3gMjcIyrgJjQKu2SNMFHrRFjBWJEEYQJx8wrkW0S3+YNkNZ/KEPQrsJ1FOIqRBN7Ww+H0wrX9f7VWuFq0JowMsJ2Ddx2K0UC4zXvt692my2AEVRYyTYZPuq0ubA6c4H1jkxmb221TcaaKBpiUpUiNmlTavbIHtFqUEDwkpvpDCJl+O6G5LXwKtC25hJaPpSnkyioIEhY062adHfIdDoWFNDwhbd1bU0gd4xHQhqemqZtm2aoIIlD41lGWr+376Oqi9R0K6ttybeQ1GwvYsouGuj6lg13zPHKxY9//GMajQZvf/vbx2z7/PPPs27dOmbOnBluW7ZsGYcccgh777033/rWtzCMLU4DtrxrxrKstgpfgIceeogvfOELLFq0KHYx/5kxf4e5ABS6BHPxEEzuQQ90ASCJlVUFpWeDMFQb8FzWMYR6ER+JKNO2SLlkArd5VC+StJTg74s57L3QXHFVopNsuNqM6TuUErS0zN1a/IxsCdi+mNU7BkxDt9VLhBijOmtrTO/EAgLSEJPA2RJdLy3/XLNTpme7KZILXyz9e+S4IOzXTFUrbg+hRWwMf9SmmKGbJRwlsFahMCLXJOpMyqq8C971LSkdZjsdi9Zl3XnB0/WIH2DuYFAThSm6Qy0ZD3XxSAiQSuAGvjYHI5b8LTmzmhgU0FgqO4w7EGePuiWWNyb7bi5N0zFYXyvj+H+DdsnAWGfjOLl7Jkc2RkZGYlaIxYsX89BDDzF58mS23XZb1q9fz5IlS1i+3BM9P/744wDMmDEjXA9PP/10tt56a6644opY39/85jc5/vjjUwLUkZERPvOZz3DiiScyY8YMnn76ac4//3wWLFjAkUceCbRIyJw5c7j66qtZs2ZNeHyndfjFxhYnIk8++SSzZs2iXC6z3377ccUVV7DtttsCUK1Wedvb3sZXv/rVcV+kRqNBo9Eq+JUMt3opoLe3F3EVVo+3BFnrR3BdF3dKb8w1gwhKg9lQuP6dMjS4Lp6FJDDpG7TCdoO3WZsxM7FC3IUTyzkS+TeoPSNhpbZg5U+4lbTyCwy3DnaaJsrQGKaErhtFmjhYCYFqM0x+1oJS3kJtGf6yLlB3LJTSGKHAVygYGkeMcGFTCK5WoVUpsHkYfnbNADW3gKAoRkhGXUNJOVSslm7A8JNuRY9taMvLuRHpXyEUlINleGNZuJj+mHG5J2FBPG+RFppCRobRLJVL3BITzSHSFIuiT0biZEPFlmsjIttMF7trjRREvGzq+1NgDQnGDL6wLiY1MbDEK3AXfD+09sK+62JSj2hKgusVjTJyUb44tfMMRsWiDyd1hoGg9q+NbRiWsj9H775usHvY2KhgGBrRCu3/jY6M5ILVHNlYtGgRhx56aPj7ueeeC8AZZ5zBt7/9bX7+85/zzne+M9wfaDuiUThLlixJWSsef/xx7r77bn7729+mxjRNk7/+9a985zvfYePGjcyaNYsjjjiCz372s6Go9rbbbuOpp57iqaeeYptttokdL+Pzh78oGBcR+fKXJ17k6Z3vfCe9vb0d2yxcuJBvf/vb7LjjjqxYsYLPfOYzHHTQQTzyyCP09vZyzjnnsP/++/PGN75x3ONeccUVfOYzn5nwfP+RUEpRkB4YGMaogLYtjKpNfRsDUQqz6okVpeSt3IYWwmVageVbhQXP8qE0mBpcBX7xVQwXz70xxrM5JWkY8y28jU6krQFE+ZMIHv0alZqTpFw4Ip7+I7oQBz+6YoQumbrrLdym0SIdTe39WwrJjYCrMDPm3ari22JxUWMQgCOKqlsImzWUQcVsOWgU4r+JZ7xlC9jaa2UrkwGj6ulIIifmZXNNunUMHPHcIlFEF19IuoMS1xCvvk3BTxbWWry9cOVWQrFIhAvpKr1AKNx1UBTHUDAnSVoAV1QrZNYn0hEFFA6mRyBFWOX2sdSZQlnZ9Bh1eowGXUbDFzGrkDAJnnWjEVFuu6KouyYokz4zruVw8Vw0sSyz4omAH2hsy4guh6RNgV9FuZXUrTLUZOqyOk5PYcznW45XLg455JCOC/s73vEO3vGOd3Ts4/e//31q24477ti230qlwm9+85uOfY5n3C2BcRGRj3zkI2yzzTaZ4T9ZWLp0Kcccc8yYf6jRhC+77747CxcuZM6cOdx0001MmzaN//3f/+XBBx8c15gBLrzwwpB9gmcRScZ+vxTwxhPfwE9+diO982BwSRGxDJQyQvGoW1ThamDYhK+AgZVDAJ24e4ZLi7CgUK6MnZAsXo4kjeSqnInxLEwKpbJdJFk6ElDp7QIo8aJGTB3bHPzUWqjjC6qjTcyOYcGtiWkxMPwbIQKijKCeYFjfxSDumjGzcmKIhOn6AZpiMuyUAZeiEY02SWo2/Hlg4IokUqG3zjmZFyXL0qR9QlMynNg+zyoSEJF4zy4qVkwOJBL10y5w2d8bWORIfyuaEbIm4hGgpN6jqkuscAawxUQh1KVI3S2w1u0DhC7VpMes01uoYykvJHup3c+wVKhpT7hrY1HAYa3Tx36VJ5lcGI1dpZq0IodEYESXeaAxm6Z4ZDOM2FJeArbAc2nYmu0XeRbWD73/7HE/D3PkyNEZ43bNLFq0iK222mpcbTf1TWFgYIAddtiBp556iocffpinn36agYGBWJsTTzyRgw46KJMtQvtMdy81HHbkwfz0lp/SM7PO4NIKuhAR0pkgVnyBMWxwS7TcKAYpa4fShLqN0HPSgUhIJBV7dir41BEk1r4EVCiszTo0a5fy3RxaFEok9caf7B4UWhueCYjg9LKWPa+qahCyGVTsjefDyIYrcfdEdAEPCEB6sU5DoyjEBvEK5lm4OL6Ly0BTNtslKfMys3ruirhrJUup0S7FuouJrRWFyLrZImqSOC5IxtYibSbRasrS6SvV6iVhFdECduQL61kYWiYRVxRrnD4GpSuYBrZjYhqBkNbbWJUSVadIQWt2KK9mUW1bhqXsfwNasxrV3jPgwfoc9lHP0GvWMYKvp4JBt0KvUWet28vDja1jAlcjmXMEQIS5Dw5Tqmr2P+5f2GPBzi9Jt2+OHP+MGBcRufjii+np6Rl3p5/4xCc2KVPbyMgITz/9NKeddhonn3wy73rXu2L7d9ttN774xS9y7LHHTrjvlxpM02TvPfdl0V/upmtqk8FG2XuDNhS6mH7MG7bgdClvQTey9R8CXt4Rw8sholDgSmYEjZCwloxJQmhPQiLbU2b9qFUjxUQiAkqtsBIZV6N9RQ/rlK20tUoqP3KktaA6WlGMjZG9cEer43o6lZjNBfCSpFlmZ7GiSLabowWFxqTuKipms5XkLdHGFpMijkemlEq1CJDOYBLMGb+In5f9NEiaplGxXCkBwqq7/vZ4DpTWkj+WRznkwuJFsyS/PMHXb0SXWOkMxKwjIrDR7cJ1DLYqDqfq//xhaHtu27gL21Q20FeopzQ3QeXkOkUebsym16zRYzToM2oUafLjjfvSZTTotRp0Gc1Y/9HrGGyftrjO5BVNJu04m9P/9cQxzjxHjhwTwbiJyERw4YUXjqvdeeedx7HHHsucOXNYvnw5F198MaZp8ta3vpVp06ZlClS33XZb5s2bN6H5vFTxlrcfz/0P3kPv1jU2LO1HlE9CMlYusylowzP161Z+qxgUPhGJpp6IpnFPNfZ+lLGWFMFLF99up4r/mg2VWL6SupDkot2JGfnWDZ91xBwcMfdDnLC0soSOMVXfmmIqiYkao2/4bsLCkrQSBHVtgvTynaAxQi1HZip4FHXtVTwsdoimSVpEIlfaT8Bm0BQDJZ4lxlReFtX0pfZIiknQLq09UQmhL2TfMdf/eqx2u1GApVws5dWYcbViuTPAsFRSxxkKBp0uNtoVphRGUwnhNjg9jLglDEPTV4jrQBQw7JQp+fqfQbeLitH06/tUsLXnIqrqMtVmmaJymFIYoWLa/tjxEyusE2Y+Mgp9JT71vnNeVqkEcuQYD0444YQJH3PDDTeM24sy4aiZWq2GiNDV5ZlQn3vuOW6++WZ22WUXjjjiiAn19fzzz/PWt76VdevWMW3aNA488ED+/Oc/M23atIlO658SfX29zJ6xgOdXP0mp18U2iUfNRKDEIxntSEjYLkh6FloSsq0iKWvIWM9Wczyum6Dz6PiRYSLm/Sy9SLIOTWxxj85dwHENClZLaCoime6NqHsGlJ9CXvu/taMinjXFVF54UrZVQ3kZVkPykGU5Cawy8axxWetYU1tYpo69jQcLvYsKFz/bD0HNTu8ev3btLCeea8nEFRMLl1IGuQmIiKV0mmSF47Unc4J3zda7XdSkwEbdRfSOisCwLndQm8CgU8HFZEl9Mtt1rY3ta/gCqdWNPhZ0rcEyWpqeUbeELSYlgoKInitmslX1xnXLsbk0xWJFc4Cy0WRyYRQzuF8C6lGDmfc2QMEH33sWPV3xuh45crwScMstt3DyySdTqaRfGrLwwx/+kJGRkRePiLzxjW/khBNO4H3vex8bN25k4cKFFAoF1q5dyzXXXMP73//+cfd14403TmjsLRle9GLhlNNO4OovXEn/5GGGap3dX2ZTsLNEFJGfDQG3CRRoPWuTVhFFS2syjmva3hpC5tt0bEJJl5CvSwh+jkL7dWmCsM0sB0BwTDwPRAbr8TcERCAgE9HaKZ3yVmgUjlaxt99k606p3oP5Jd0c7dmcou4WKBm2bznQ3lz99sp3azXFwlR2LA9ItOdAFjQe1qgFhqVMj9HwFt3IIS4GBfHCjZM9Bfcl0CGF+WAit6sqFsNukSftreg1G5EGLRgIjhiZpKrqFsOaPWvtHqY5w/SYjVBPZPv7NIpVzT62Lm8Mj13Z7KPLtGMsdtgpUxIbDMWzjamZ16euiyxvFKmoBgyCcVcBtcLA7DU56uRjedWCdNbnHDleKfjyl788bmLx3//93xPqe8IpAR544AEOOuigcLDp06fz3HPP8d3vfneTwnxf6dh+h+2oFCbT1zVKAQflCsqJfwzb+5Q3aJTrZVc13MhHtz5KSKVtV6Ji6eG1EtCgXOXtQ7XPeQ3+AqOyX3+zXAmpdq1Jua5CxMj8uNojCVr8PiSrrzEGT0AHw0eOcbQX8tspTbgWRdN/6w7LwKfGUzhZibEk3saNmJ/aW2E8F01TLOpSoClBhlA/zFa12jhiZiT18tq0qzCcnLsI1KVIVco8b0/y69nE5z2e7KkisMQZ4DlnAAeT9brMUqeXNW4363U3doeaMqNuIbM2j4hig9MVm8uztanhb41YuJhieX0gzMq71u7BFqvlCBRYMjLA/zy/K3eu3YEn69PbpLn34GjFYyunIzeXUSsMFuw7n6s/ezXH7du+pHuOHC933HHHHRPSff76179m6623Hnf7CRORarUaRsX89re/5YQTTsAwDF772tfy3HPPTbS7HMD+Bx2MUjCzZw2GKzFiEZIL8bjAWEuvKF8jkny6+y/uIuKRj+gCqnwXTqd1pxMZSU0i+oPvhonuzuzDi45wnJYPKcjf0Pm4saGJXzgd0Th4+UoUtja82i6uSVObOGLhYuH4gtN2Y9sZquEkGXDHjKH2jxN8C47CFssfM637EFG4YmL7NXkcFwztdnRzJOHQIjMag426i41uJayO7G1PZ49NzmiN7qZGkSFdYYndx7AuhQt9Xbx07FnXTilYaU9iaWNymt4pYdCJm4BrusjKZr9HoBJx66NuiWG3jCMG62zPqqgQhu0Sf1g5n/vXzqGpLUxDY+vskNuma/L0+qn84dkFrHx+EpYtzN5/Vz525vl0Vboyj8mR45WC173udVjW+B0oBx544ISiVydMRBYsWMAtt9zC0qVL+c1vfhPqQlavXv2SqEr4z4iTTj6KKVvNpa88ykClfUhgs9/ouBoLfoKzTOe9Z/VQOrooxxspyTguZZcfHxkR0ShDe0XvYuuyF34bPw3PTRInHSqWAj5LGJm0PGQj3k+wbdQuUHctGtrCFgtXTLR4qb2iNXea2vLTkmfDy/URHzuZPlzHFuO0q0nEi/JwIgu/0EplnowYCdw1rphUdZEh3cV63UujzSKbmrMQyVQauKgUdSmyxu2h7ufTcNqm5vVI3KBbYp3uDnOrBCGzAWq6CBg4CbW0AKNukWG3zIhbZo3dk7q/SSICsKzu5RepB/ULRLDWCd3/5/L8LdN55ifboH9RRn5WYs2PpvDInbNZ12hpOp4ansaoXYxn07UtHlu7FX94dj7PbJiCKyaVmhcNtcf83TpcxRw5Xpl44IEHePjhh8Pff/azn3H88cfziU98gmZz02owTZiIfPrTn+a8885j7ty57Lvvvuy3336AZx3Za6+9NmkSr3Qopbjw/A+CKjG9fy1FK30z3aJntSiM0pGMBM/oqHtGwMstIvFF1mvY9pf2CMhI22koDFMwzCyRp99FVEbSpoy8SPbbNP7Q7gSq9CbH9lwCafeGkWquqLnFjuPYkYJ23r9pF05AVpKWIY3CISBA8WOcMG9/HIHOpSFeBV5bTJSCId3lCTzb1VGh5ZJpzaSVYTW4AoO6i+XN/pDwpGegqInFcrf14uFixIiICNRCQpNwvwgsbw6Ec1jV7KeuC+E1bOoI2Yiet2vw7Mqp1J4r0/2gy+RbXCb/yqX7bwLLTfTqAsZGC3OkgDFosM3jo5iN1siumDyyfiaDjQrDoyUeWLwtdy+Zz/ODk2hVKiYkItvPeuklQsyRY0vjve99L0888QQAzzzzDKeccgpdXV38+Mc/5vzzz9+kPicsVn3zm9/MgQceyIoVK9hjjz3C7YcddhhvetObNmkSOaC3t4cz33Um/+/r/8HWAyt5du024UIpgN3t/Ww2PR7gZIj3dWR9MXQkg0Z7YwHhANGfgTYZy1tILXbxTkQrT+SqsslIsL9dttXOEw/2tIiKS3aUcjieL8YMarqAwtZC0YhHjETPIoBGUdcFym3yhjgR10vg2oguvF5eEgvBxRKNaQXF7eIXObDyRIWtHlmJn7cjCltZBO8RtjbDuTWxWO920200qCg7IZL1UtEn9SVZ+Ueetyfz98bW7FJaxg6l1RQNN7w2riied/pjriBXDGpSDMOjm2KF319HTJCWeFSjWGO3SIygWNKYwvaVVThVg5WrB7DWKNQwGCMKYxjUMKiqwpEeHKAbQXWbzFu4gCP3P4zdd9gtVpfj5t//hlv/66dMW1pjxYLuSL4ag7+vnw7LS4xu7Ia51dRXOSAic2aO38edI8crBU888QR77rkn4FUCPvjgg/nhD3/IPffcwymnnMK111474T43qejdjBkzGBkZ4bbbbuPggw+mUqmwzz775PH1LxD7vmYP7r1vPx75y5+Y1ruO1cPTEMCpqFhYr9UAlJ/gzEfglolt8HUlMddIwkUQSD8ChNV4s1bkJBL1S8J+FaCCcvPZoabe4qzDGjHtoP3usjOPttKxCyYNt3VO0LK0uHgEoWS6IQkJ5tAK782yhoQnSlNbFA0ncx4iiqY2qBguI24xw9ISTBbKRoNe3cxkZ2n3lEYrM5v0iQ5r52gMHG2EIaygGNUlGhToMesUVFCThdDtEj+7tLuo6rPaxxuzKBg2fWaDKUaVgnJZbE+JZ6hVgYVJUZMi3arpj+N9iRzfYhP0vcruiwlGxYHas2Uefnw7mku9ccvR+ZUUpf4yPdv0MnnKZKZNncJ2s+ey/6sObFvC/JgDXs8vf3oL05+rsWK7LjAV4iicmoVdLaDXlzwblKOgED//Ss1FF81xhyrmyPFKgoigtfesuf322znmmGMAmD17NmvXru10aFtMmIisW7eOk08+mTvuuAOlFE8++STbbbcdZ555JpMmTeILX/jCJk0kh4ez3nsa55z3OFNYz2izi5FmN04lsRCJYNV890TFc5PoErFFVPDcM9Ip70iWRYTxkpFg5xi+EfH6SleR9YnKGJEZjvZqjhStdmE9ylfxtibbckO0fhfXwDS0d3aRuThiYAZ5NDJJhnd8Qbk0XStd80U8a0AdC4dGexISnrmX5r1s2Km+km4bwcy4dp5FQmFiRvxvDbEwpRmzgDiYbHS7qCibotGkqaKZ7qIziv7m6T2aUsBQgovBoNuDoRQjuuTVXskgUa6fgbaqS3QbTWoR94+DERJWpWBlc8A7v7UG7mNF3CcL0FRgwKTti2w3fzfmzpzH3Olz2Wb61mHOoomgUCgwsOsCRh54gkl/0yztnYQ4fhSSC4Y/Nz1sYUy2Y8eWaw5mX54vJEeOLLzmNa/h0ksv5fDDD+fOO+/k+uuvB2Dx4sVMnz59k/qcsEbknHPOoVAosGTJktgD4i1veQu33nrrJk0iRwumaXLBxz6MFoNZ/auRksa0waoLVk0wa4JVB6sO5Q2COSpYo+kCeAo8U0DWHR7T2iEtjtGOJwQLd9L5Hx0/vTmyVzpqQKJts7J4JklXJz7kG2houqZXAVZHNSoK2wXbNUNRqwheJI1rYGsDR3vJv7SfBCz4ONrEES//qItF0x2L1wuG8urURDURrZlIOH5AD3TiRnnCVm/8KOzkFyDSa00KPGtPY62dzlMTXJvo70ubU2JWkrVOb7h3rdvLWrs3nQPGdzON6LKnD9HF2CiBO2jdcA9Df+ml+d/dNH/Sg/u3IpW+Egcdsx9XXnklnzvvOt7zpvdwxGv/hR3mbb9JJCTAZ08/C7dQYO7SjRjNFtOMpn6Rqolbb11LwxVKDZfugf5NHjdHjpczrr32Wh544AE+8IEP8MlPfpIFCxYAXjqP/ffff5P6nLBF5Le//S2/+c1v2GabbWLbt99++zx8dzNh61kzWLDrvjzz9z/TXawx5LaKCEaXJV0AQ/xqqFohgeQgMGZo2ognOq/+knxF9gvptSYhGa6ejAn6G7ycIPE3e6U0bazq6fngnSNCtp5ElJfJbax+xED7K6grQc+gMDAME9eBghEkvI8P5Ei8yF2ULASzdMTEEKet5iXM2InyytHrQsQyIqGGJdqv9lPNt+bh3WQNiXT1fk2ajARrIorVzT5EoN96PlH4L35Oo7rEerebsmpZCUZ0mYa2KBmOF6XjlpleHIodFxAmB5OaFFKRMuuW97J60RRWPTkJtEIVYcE+c/jXQ45jl/m7vihu3a5KhYP+5XD++D+/Zs669TwzbWrr++zDcAycmgkCZsWlXHdQwJQpUzb7fHLk+GfGM888w3bbbcfuu+8ei5oJ8PnPf36TK1JPmIiMjo5mvqWsX7/+n6Lq7T8Ljv/Xg7nm73+mYjViRCSANsAph2oIDA34QtZAO2lo0DYtghJFkmyoNvuC/WGbxIIfFoWLkBOSFVqTphPB8AUOWrzqrtFjEzYAvKgTw3NFZJCRsSlIq+e4q0P5M7cwEDQGWtpoWiIF7NIkxMNY+UKStXVcFDVdoGK0q74bt4hoUZHfBdfXxwTwdCxpItLQVhiKu87uZlpxJDKLuNZjSXNK+HMLwlqnl62LG6i6JWq6zKhbosts4KmAol8wYUi3tBVOzWDF3VNZ/3A/oBiY08PBBxzM4a896h/yvHj7Mcdyz51/YM669Tw/MOCLfCNwAVHougki9A96WWBnTU/XucqR45WM3Xffnblz53Lcccdx/PHHs++++8b2l8vlNkeOjQm7Zg466CC++93vhr8rpdBac9VVV3HooYdu8kRyxLFgu7kIii6zntonCuwulWkeCGrSBOuTApRDfLUeyy2T7BPlP7CDDhNHpKrqZpAFHaElRmJBDl004olXY8cGmg6PjGSTDjWma2Z8UClXSGxfGBGTcd3HPD47OkVjUNMFEElF27TGDVwy0QVfpd0zYma6uoad1gNiha/PSM5JBNY63YzqrIeJYp3TExGxwiq7N5xsnIApqrqEaGHtX/p57FvzWP/wAFvN7eKCT57LlZ/4Am943Rv/YS8tpmnyphNOwBJh/tq1XjbhyD2ytGb6uhF2fXY1B9//PDv/fQMAO8ye+w+ZX46XJ/7whz9w7LHHMmvWLJRS3HLLLbH9IsKnP/1pZs6cSaVS4fDDD+fJJ58cs99ly5bx9re/nSlTplCpVNhtt91YtGhRuH/VqlW84x3vYNasWXR1dXHUUUel+l25ciWnnXYaM2bMoLu7m1e/+tX85Cc/GXPstWvXcsUVV7B69WqOO+44Zs6cybvf/W5+8YtfUK+n16mJYMIWkauuuorDDjuMRYsW0Ww2Of/88/nb3/7G+vXrueeee17QZHK0YJompa5J6NGNRE0WAtiVBAlJ6iVUNHLEIyXapVX4ruPK3EZAGggx2h07pu4ksEZkhev69V9EUB3dSApXwDSi/Xl7om6K5PFjhedG4fpVfVNCUgWOQGEMVpN0pST7yDwGg1FdpIc6ZsZJVN0iIBQS7qegunBkhDC6JxjLFRWSB/DCatc73UwpjPpz8okI8HxgDcmgQy4m65weXN/lst7pYU5pHSaCE7miIjC4vJvld2xFbXWZUg+ceMYbeP1+x26xqLqjDjiAX/7612y9di3Pd0/yXC/VUabURuiv18K3MentZutXLeDQha9l/9332iJzzfHywOjoKHvssQf/9m//llm59qqrruLLX/4y3/nOd5g3bx4XXXQRRx55JH//+9/bWhY2bNjAAQccwKGHHsqvf/1rpk2bxpNPPsmkSZMAj9wcf/zxFAoFfvazn9HX18c111zD4Ycfzt///ne6uz0B9umnn87GjRv5+c9/ztSpU/nhD3/IySefzKJFizrmAiuXyxx77LEce+yxiAh/+tOf+PnPf84FF1zAW9/6Vg4//HCOO+44jj322AkXrlWyCZXkBgcHue666/jrX//KyMgIr371qzn77LOZOXPmRLt60TE0NER/fz+Dg4P/dJlfL/v8V1jy1MM8NjKXpr+YOCWFWyS2qomC5qQWCQmhwI18p7UJuiBIKX7LBVpF8JRk2slEiVeBN4somIIy4tsMK+NrpQTT9LKtZq9JQqkQdy0kiYtCU4hU6A1qxijll6tPEgg0hopuloyQYcGKbDODviL7FaBFKJtp1020bo1lOBRUupqtQigYbluSAmCg6TPriG9ZEf9TdS0ERb9RS+hqhKJyqelCGA5rKpc+q/V2MmiX2eh0Eb0CZWWza/cylIKyamKhWW5PYrk9KZxHyUjnTCmrBtWIxWROcS1bWUOsXDXAqiWTGF1WYfT5Ck7NAkN4zcE78Pbjz6LyEkiR/sDf/87XvvQlNCq0AmnDwJg6ld12fxUnvO4QZm01McX/P/Oz5eWG4F4cfeu7KXQXxz5ggrBHm/z6qG9s0r1WSnHzzTdz/PHHAx5hmDVrFh/96Ec577zzAG9NnT59Ot/+9rc55ZRTMvv5+Mc/zj333MNdd92Vuf+JJ55gxx135JFHHmHXXXcFQGvNjBkzuPzyy3nXu94FQE9PD9dffz2nnXZaeOyUKVO48sorwzYTxZNPPsnPf/5zfvazn3HvvfdyzTXXcPbZZ4/7+E3KI9Lf38+nPvWpTTk0xwSwy047suSph+ky6jR1EbcAbiljBW/3pp0gFIbr6UUkyywQxldm9K1894xWiKFTZEfp8YlFPX1H+wRmimSl2nQHafmKJ9pQZM8h65Q6j+FZNaJWkZYKxsQV8fOORE4qcRYa5eleYls7w7snQpN4vHXLFaRoiEU5JoZVsWRq4LlJtHgF60RgxC2nRq9LgY1OF5MKVZQvsl1pjx0lUtNFxBWcNQXs5QX+tmIuf1lRwGn4jxElDMwostXuRY4//J3M3/alkyL91bvsQmXB9gw+/zzztt+e1+27L/vvtSeFQjqvSo4c7TA0FC/BUSqVJuxmXLx4MStXruTwww8Pt/X397Nw4UL+9Kc/tSUiP//5zznyyCM56aSTuPPOO9l6660566yzePe73w1Ao+Hpm6IWFcMwKJVK3H333SHJ2H///fnRj37EG97wBgYGBrjpppuo1+sccsghEzqPKLbffns++tGP8tGPfpR169axfv36CR2/SUTkrrvu4mtf+xrPPPMMP/7xj9l666353ve+x7x58zjwwAM3pcscGdh371dx6y//my6zznrdFxGnppFJAzKam7ZCLEGK8WZe1naJb2zTR3wIP2qnc7MQWquURSJqrUj1n/BAtSvCpjOMNdkVYsaer/hHhvLLSFSIm4ie8TQsCZeJGJgZ0SudoEVRzMioFk0Z35ACBdGYkSRxOiUEVjS0RdmwqUdEqkksb/YzYFUxFCxuTo4lGFNIjKyJC81ni1QfrWA/X/SSgAEYQt9Mzdz5M3nNLoew2477vqBw2xcb137svC09hRwvMh5fNQ2za9NFk+3gVj0r4+zZ8bT/F198MZdccsmE+lq5ciVAKufG9OnTw31ZeOaZZ7j++us599xz+cQnPsH999/Phz70IYrFImeccQY77bQT2267LRdeeCFf+9rX6O7u5otf/CLPP/88K1asCPu56aabeMtb3sKUKVOwLIuuri5uvvnmMAx3PFi+fDl33303q1evDpObgWf9+eAHPzjhqLMJE5Gf/OQnnHbaaZx66qk88MADIQsbHBzk8ssv53/+538m2mWONpg1cwaiLCpWHbuQLU4NkRSHZmyL7YyttJGlrI0gNYQmXPFbFgPlizSy5xKdpHYNTMvNbKaUJKwVnSwPcUggWB0HIxpPIrUgIkVEQotEMIvoHLPIjiQsO1HXTUZjT+uBkSlmjVo8BE/v0W02McLonlT1IJraomLajDjtH8g1XWLILWMgrHfj+UX8bnE2mtT/Xqb+WAVdNUAJPdsY7Lzjq9hr533ZbYc9KRY3vxk8R46XKpYuXRpzzfwjI0W11rzmNa/h8ssvB2CvvfbikUce4YYbbuCMM86gUCjw05/+lDPPPJPJkydjmiaHH344Rx99NFEFxkUXXcTGjRu5/fbbmTp1Krfccgsnn3wyd911F7vtNrYV89vf/jbvfe97KRaLTJkyJab9CojIRDFhInLppZdyww03cPrpp3PjjTeG2w844AAuvfTSCU8gR3sopejqm4JsXJ16740io8aZp2nIsh4owXANdCNqFRE/Mam3mIsrXlbWQNwa7dfPW+IfFt8fTDBGJILV2Lcz+CfSthhe5EzHYw2JDhgVrHZc/EmG8WbNA7QOBo38oSWsJe3gihFaLlqWlTQB8qw5XnG75P2Nh+sG7RWjbpFus+GJezNcTxqD6voio0u6oK6goaCB/6/yt8GTei5GUSNFUEVBFQWKglnU6HUm9jLvC1LqN3n1UXtwzOtOYOrkiYnQcuR4OaGvr+8F64FmzPBCw1etWhXTVa5atSqs4ZKFmTNnsssuu8S27bzzzrGIl7333puHHnqIwcFBms0m06ZNY+HChbzmNa8B4Omnn+YrX/lKTEeyxx57cNddd/HVr36VG264Ycz5X3TRRXz605/mwgsvbFtiYaKYMBF5/PHHOfjgg1Pb+/v72bhx4+aYU44ICt2TMQZX0UuVIbpbAo/x+ELSgoqwZo3SIGGiMpU4zCMkygGUeHlIjFYfCAlxanIc37JhSOi2iNdQSZMKpfw3/CyNa0dS0UKg7dg88EhAkgZ4Kd09DYarAWWgpJU3ReGFG2v8vCeJ+cRcHhKEvnonnsxz4uh2lQcVVbdEt9mIWGQEXTewny5iP15kcGW29kEVFVbFxOgSlNlEOV04NQd3o8b1iz47AIYw+1XTOOqQN/DqXRdutgdOjhyvdMybN48ZM2bwu9/9LiQeQ0ND3Hvvvbz//e9ve9wBBxzA448/Htv2xBNPMGfOnFTb/n5P8/Xkk0+yaNEiPvvZzwJQrVYBUn/PpmnGXCydUK1WOeWUUzbrM2HCRGTGjBk89dRTzJ07N7b97rvvZrvttttc88rh4+x3nMTll/07M9z1VJsVUH4Ei1KI8rKrepExySUzsUnwfAARd4ppg1tsWSky1zzxanOgQYzA5JCOOom7egzEFQyzjaYjJW6Vtj8lSUhSjxGFFgPbEQylMQ2/vgk6U8Qq0tJ2tIsb0yiUSDzsVPluG3FxxAipijdP8Wrj+M0tpVGiMaNRPpGudFiR0DveFUXBaO0T0o6f4DJ7lpESXdRxlhSpP1bGfq7oZS21YParprP7Xrsxf/rOTOqdTF9PH5VKpWPmQ601jUaDWq1GuVx+Ses9cuR4KWNkZISnnnoq/H3x4sU89NBDTJ48mW233ZaPfOQjXHrppWy//fZh+O6sWbPCyBpoVbT/wAc+AHjlVfbff38uv/xyTj75ZO677z6+/vWv8/Wvfz085sc//jHTpk1j22235eGHH+bDH/4wxx9/PEcccQQAO+20EwsWLOC9730vV199NVOmTOGWW27htttu45e//OW4zu3MM8/kxz/+MR//+Mc3w5XyMGEi8u53v5sPf/jD/Od//idKKZYvX86f/vQnzjvvPC666KLNNrEcHubO2ZreKbNg3XIqbpNm6E/xQxDr0JiUZqYxj4LQWuET5ES5eIXxxoKAcpXntkmtZdGBgl8V2jEwLZ0iOOkMp9mEJVmNdmyoUMPRshMo3xsSjYLxLkirtozCcZUf0RMN442eW1QPorC1FbNGBHQiSqXqboGiYSPipixC3gxaZibPiqKoOl6W13aZVoNrJ6sM9JMFNjzd47lbEKbM62O/1+7P4QuP3qTKsYZhUKlU8qqzOXK8QCxatCiW4PPcc88F4IwzzuDb3/42559/PqOjo7znPe9h48aNHHjggdx6662xiJenn346Vs12n3324eabb+bCCy/k3//935k3bx7XXnstp556athmxYoVnHvuuaHb5/TTT4+ty4VCgf/5n//h4x//OMceeywjIyMsWLCA73znO/zrv/7ruM7tiiuu4JhjjuHWW29lt912S0WeXXPNNRO7WDDxPCIiwuWXX84VV1wRmnlKpRLnnXdeaP55KeHlEOt/95/u53vf/n+saUxiRSPuo3cL0BgwSNYzE8OvyAstawiQlTjTLQliyrjz7IqpoZzxrp7IHaKUYFja04ckFlbDaNWaCd0yeITE9NOUW2a82J1nDdAUOpYzEApm2sSoEAzfMtLJfWNE5mIZbowQKLRPTlpjeftbxMHw85YE+72cJJqS6aJwKZk64ZLxYCrPwtLUFiKKHquBqcQXsvpC2SbIXwvIkxYMeRevMrXErq/ZiWMOOI6ZW8XrP+V4cfFyeLa8XBDci+2+84kXLWrmmTMuz+81nk7005/+NDvuuCPTp09PiVX/93//d8J9Tsgi4rou99xzD2effTYf+9jHeOqppxgZGWGXXXahpydd2TPH5sEBr30NP/jhfzFZBlnVmBKGWgpevRnDBsOJV+CVqJFCJ9/V41COL0wdB0LeOg7FpojvAgktIBIe03LPRASX/j5XmxjKxdUqRWAEA2lTD8brIvssW9EuyTo4necfbZhlyYnaNbLmEuhFqo4CTApGHVcbiSyjgq1N7IhpatQp0mM1wvsoowr96xKsMzG7DXY8cHv+9aBjWDBn+y2WsTRHjhyvPHzhC1/gP//zP3nHO96x2fqcEBExTZMjjjiCRx99lIGBgZSCN8eLA6UUBxxwMHfd8WsmFYZYZw8Anj4E01MpWFVoRol6hIio5LZk/0EdmU4IvS5+XRgbz4JiSodVPUp//JBWf0IeMRBMM8v9IrSnTuOLWml3CvEYmIw2Ue1Iag7ZOUhi8pjIMUFSMUH5kTGK9Q2TbqtJIVKsTkTFSAj45MUtUjFtZINC/08ZRg32OWJv/u1N78rFozly5NgiKJVKHHDAAZu1zwk/zV71qlfxzDPPbNZJ5Bgbbz7+KASTKcWNBKGe0QRnhgtGg7gWJKoN6QCFQrUTTGf1EUTaaAVN5ZGSdqt77FjlhbD6lgmts7UhQeRJUPE20cPYnCmzgULEaCtMTU43sKwk22viczJodw6+hkdaJCRAXReoulaoUdFZ8dd4OUT0MgP9swrUDI592xt414nvyUlIjhw5thg+/OEPc911123WPjcpj0igB9l7773DQjoBXun+sxcL5XKZ+TvtzjOPPUiPWWXQ6CYhWKBQg4avZRVF6JIZF4JI004+nCQEj4xoBYb2rCPJJqJQSauCCAoD5ScMy7aIZKNT1Izn5pnICWSOEJlTlgWktS3MixKbQfB/SelAooJZV0xGXcOv6ZJ9l8ynDOTuMljwzrPewWt32+8FnFeOHDlyvHDcd999/O///i+//OUv2XXXXVNi1Z/+9KcT7nPCRCRQ1h533HEx37T4YY6uO7HU1jnGj7efcjyfufhBphY3st5Ka3KUBlzCuzoR94XRBL/Ia0dIoPOICGABzzJSIGVjS1sgJFa4TcQPRo1ONkI2oq6SSAO0pCNLVPSYDJIAnavjttoEFX4Jk5C5EeuNV9lGh21SZygkwpZaM1TEXTIjThFTQdmMFJkTsP6qKP6fgepWvPusM9l7wT4d55wjR44c/wgMDAxkVhR+IZgwEbnjjjs26wRyjB+zZkzHxcIyHc944cbiVFFAaZ1XidcwQCwVFzB0cJ8U1wMboDaLjAq7EXhJQOMkJOjc9RbtuBEgaqHwI1eSVoQEGYnvbqV9j+YUaYlPM0+no26lXVZV71L5sTk+0REER0M0tZmn+UgLZj1LD9mTas0shPbDgG1PgkvRdEFD4Y8G1uMGpWllPnXuJ9lq8lbtTiZHjhw5/qH41re+tdn7nDARed3rXrfZJ5FjfFi5cjWWctjgdmO1MTwZQGGtZ91o9oAu+pYOA9BgjUJQFV45niVENcH0S/VaoxpnLO9a2wR8yrPIJASsIWHooBRtkZG09SMgDvF9nV0wnSrsZoldg1xvYTE58S6YoYwYAQrgonwribRmI+BiokVTMNJEJdq/CDTdoNKuouZaKCV03eW5ZAbmTeLTH7qI7q646zNHjhw5Xm6YMBH561//mrldKUW5XGbbbbf9hxYCeiXhT39+EICqXRlTZuyF9AqGqzBrIJa3ABerMDrVk3QE63pUqGr5xc2cniz9w3i0F8ovjJdMACIYY7hEAjdOVvr4ZCht0C0h4YiSBRWGGQd0JZrMTKPAJzfRVsn5aTEQ8banRakmWlzSJee8iBdXBCt1vq3fG9qKqU8Eg4ZjUh7xLDGf/mBOQnLkyPHSwKtf/Wp+97vfMWnSpHG1P/DAA/nRj37E1ltvPa72EyYie+65Z8e8BYVCgbe85S187Wtfi2WJy/HC8fDfHgWgUS3BGGlblBAWrlP41g/b21ccDdLC+211/GezqVCDYPfJ2PGuWRDlCS38Y0UEM8Mlkzgo/Ddm+RAwTB2zJkSPaOd18vKNtPoLSQ5gKE+pERKxDvPy3DAqRVKC6J12hzpieOndI8cF83DE8C0u8ZEcTEZnCwMrXX7/4N284cAj208sR44cOf5BeOihh/jLX/7C5MmTx92+0WiMu/8JE5Gbb76ZCy64gI997GPsu+++gKei/cIXvsDFF1+M4zh8/OMf51Of+hRXX331RLvP0QGrVj6P7Vo4TQvlSljALgsC8URl0lp4DcfzoIQNI2tsQEoMV1HcCHavtFLAT4SQaBVaRToWyAvsApn6EG9DFgkJdgYRLJ1zi0SyfKgWdQmphHS21ggGrtaYhsSsL5pORfYUtjZQYXbWYJ4KW6dFOForHFE420D//S6/+/PvcyKSI0eOlwwOO+wwxpuIfaJJFidMRC677DK+9KUvceSRrYfkbrvtxjbbbMNFF13EfffdR3d3Nx/96EdzIrIZMTIygmuPUG32oFCYNcHp6bz0mk4rzbuKaEoMm3DlVokQ35h1RBSFIbC7BNPxzRsRrYeoVhSNmH6CtYh7BBcoaYzMpGfi13bJmnn8y95e79FybnS+EkGCsax9Xi9Zrp8oNArRntvEikX9dLKoeKSjaLRqzTRCXUgctjZwMaGkEQvqo+N/m8iRI0eOFxOLFy+e8DHbbDP+khMTJiIPP/xwZtnhOXPm8PDDDwOe+2bFihUT7TpHB9x3v6fNqTU9n4rZFBwtWatrCOUQiiuiBWgVtIrdJYWn0XwieNEiVhUoRsNVWv0EC70o0F0qscYqsAUKIEZ8wTZMHboqYgXhGGtxT5xHOK2W7iOJILV7J2gMlNa0zxWmEEwUTmJbOtNqbGzfMlI0HJpitW2tUWjx3GaGA1tvndeNyZEjx0sDWWv+5sSEUzTutNNOfO5zn6PZbIbbbNvmc5/7HDvttBMAy5YtY/r06Ztvljm4+55FAFSbXmVUBZgNSRoPYggIR9L9IniLHeKJVpPHJPs0UB3HAYXb1SZ7mmNijxRxaiauo3AdhXbxs6YaGdlT21WLSaIzAUjPMPLTJuY7U2iUUi0LyjhnoTFwxPATnGXDdk1cbVDY6N2Q7WbP3bRJ5siRI8c/GSZMRL761a/yy1/+km222YbDDz+cww8/nG222YZf/vKXXH/99QA888wznHXWWWP2dckll6CUin0CMgPw3ve+l/nz51OpVJg2bRpvfOMbeeyxxyY65X96LHluGUuf/TuOa9JwiuF2q55gGBkwHFLp2xVgBu6ZDMTaix/i2zZk12vTLqJG/AQn4hp+iK5vR4lYMnSMjIjfzp9LpvsmK4qFGEHohKwmZkaIbmpMlZy3p8ltN6b4+1wNTd3+T23ULnhnpKH/Mc/isscOu459Ijly5MjxMsCEXTP7778/ixcv5gc/+AFPPPEEACeddBJve9vb6O3tBeC0004bd3+77rort99+e2tCVmtKe++9N6eeeirbbrst69ev55JLLuGII45g8eLFmOY4y8X+k6PRaHDVFV9AGcKap3qhCz8kRqEEDFvQxfYrqOGQSTeVE9eNxPbp+GJdGAVdbsNb/ORmRlPhltNaEKWUpyMBxFUoS1CGpNSlrTwiKiaICtw/ySRkE8kRkhn4E3M/RVK6t3HvmKpNjhMUWtxYtv1g+p4VJAwdwlROqu+qXaCpTQxg8hNNutZo5rxmR3aZt2P2CebIkSPHywwTJiIAvb29vO9979s8E7AsZsyYkbnvPe95T/jz3LlzufTSS9ljjz149tlnmT9/fuYxjUYjFjY0NDS0Wea5JSAiXPGZL2LLKIPPlmkuU5R6RnEml3HLJiiFVReaxU6dZG9WtLdyhNvFq+qrNO2TmPncw2gq3AwdhoTbFOIYiOl6OoyEliTQhhh+AvXAn5RlDekQLOS1CBOgtY5JEZTIvvaRL8F4fn9tWthaeeHJ/jm4onASdWa0eEnLKqbD2tEuKgUbZUBTW4BQ2uAy9ZEmxkCZj55+ducTzJEjR46XETapjOf3vvc9DjzwQGbNmsVzzz0HwBe/+EV+9rOfTbivJ598klmzZrHddttx6qmnsmTJksx2o6OjfOtb32LevHnMnj27bX9XXHEF/f394adT25c6fvidn7Bi3WJq6woMPVf2F3wHa9SlsNFG2dqzbDiCiUtJNehSNfqMESYZQ0w1NjDdXMd0Wc9kGaJXqpSkiSEa6bCwBkREaTDr/sYsIhLxDBnaiMQERzuLtle+crR9yKtXrdYnJRkTHFuRoRIkpPNR7UODIUhyNpZwVvluGttV1F0TR8zMMbUYDDZKOFisHullxcZ+RupF3KZi5r11EDjrXe/LEwLmyJHjFYUJE5Hrr7+ec889l6OPPpoNGzaERe4mTZrEtddeO6G+Fi5cyLe//W1uvfVWrr/+ehYvXsxBBx3E8PBw2OY//uM/6Onpoaenh1//+tfcdtttFIvtTQAXXnghg4OD4Wfp0qUTPcWXBP7v/r9y5z234dQM1j7q52QHlK29vBeuUBi0KQzbTNJD7FZ8il2Kz7JjcQnzC8uYW1jB7MJqtjHWsq2sZr4sZydZwu7yDHvLE+ztPsGe1WfYob6MbZprmewMU9bNIE85AIWR1nKaaT2JemIUKCcuBA0L5EUP0WqMhd1zZegMF0lYy6UjxqNEbVlcOllDDJUurJcaTbyIl6pTpKELfjbW7LZaFDVf4yMoXDGoNopM+ptDaUiz52GvZbftdx7H/HPkyPFSxvDwMB/5yEeYM2cOlUqF/fffn/vvv7/jMY1Gg09+8pPMmTOHUqnE3Llz+c///M/MtjfeeCNKKY4//vi2/b3vfe9DKTXhdTkLkyZNYvLkyeP6bAom7Jq57rrr+MY3vsHxxx/P5z73uXD7a17zGs4777wJ9XX00UeHP+++++4sXLiQOXPmcNNNN3HmmWcCcOqpp/Iv//IvrFixgquvvpqTTz6Ze+65p23W1lKp9E//Rrl27Xq+/rWv/f/2zjs+qirvw8+5d2bSC4GEJEBCkw7SBEGUIoplEdQVBRREbK8FFWGFVQRUQHBFdFFUZLGLuth3QbHAAgoWBEWlCtICSEvPlHvP+8edmUxLICGYRM7z+Ywy55577jlJyHz5VSTw+8Y4pKdULwrL9g82y4Fh97hpGHMQE50WbTuRlJhAQkI8yYkJ1ElKJCkpnsLiYvbt/50Dhw5z5MhRcnNzOXT0KLaSAup58oFS4WcgKNKi2CXrUeCJD3guwQEYEeJkhUcgo0LyhEPQIsSHBK0hJHab6Y8ZCVzseCXiffefUC0dSRnWDm+sSoR8mPLqhQQsa1ViDXQrSSsGpsQdWGjFn3ZDyh4nnhidW68ceQIbVygUNZ0bb7yRjRs38sorr5CZmcmrr75K//79+fnnn8ssez5kyBAOHDjAggULaN68OTk5OZhm+L8Ad+7cybhx4zj33HPLfP67777LmjVryMzMrJLzBIqZw4cP88gjjzBgwAB69OgBwFdffcXHH3/MpEmTKrV+hYXIjh076NSpU9h4VFQUhYWFldqEj+TkZFq0aMG2bdv8Yz4XyxlnnMHZZ59NnTp1ePfddxk6dOhJPasm8+Rjz4LuoWBXFO7C4G+RVTHVQHqratVNz0W3mVw84AoGX1HxSpxOp5OtO3byw6at/PrbLnJy9sHRQ7Qt3MMREc9vtvq4hD3YNeMNUA3zfoTUIAkTLRJ0exkRst6JNs26HskSYUoRVDI90v0nmtJrRaJEWkv6S9IHxphEXENGeqJVWyRIjAhweWzewFXfvdaf4/PdxBQbZHZpi1Z2EROFQlEDCI05jPQP3+LiYhYvXsz777/PeeedB1gZoh9++CHz5s3jkUceCVt36dKlrFixgl9//dVvVWjcuHHYPMMwGD58OFOnTmXlypUcO3YsbM7evXu58847+fjjj7n00ksredJgRo4s/UfSlVdeyUMPPcQdd9zhHxszZgxz587l008/5Z577qnw+hUWIk2aNGH9+vVhBU6WLl1K69YnZ1YuKChg+/btZWbdSGl9SFSkhn1tpFOnM1n62S4Ss5zoDsnR7TGYbutDSgDC5UFG24mOKyEhpZCEuPoMuvzCSj0rKiqKdq1a0q5VaZbGgd8PMXvefNi7kyR3IXv1euw2U/DX4IhQKNXSHwI8WAXMkOARpZYUnzvEFGAL/mD36xYh/QXFZDkZLOVxovM1Ed4d19eDVxfBIqSCWyBUjHhMzdtpNxRJvQPFAFx8Tp8KP0WhUARj7I9FnoIeZ2aJ9YspNOZw8uTJTJkyJWjM4/FgGEaY1T4mJoZVq1ZFXP+DDz6ga9euzJo1i1deeYW4uDguu+wyHn74YWJiYvzzHnroIdLS0hg9ejQrV64M36dpct111zF+/Hjatj01JQA+/vhjZs6cGTZ+0UUXMWHChEqtWWEhMnbsWG6//XZKSkqQUvL111/zxhtvMGPGDF544YUKrTVu3DgGDhxIdnY2+/btY/Lkyei6ztChQ/n111958803ufDCC0lNTWXPnj08+uijxMTEcMkll1R027WKK665lG49O/Hk7HmQfpCYum6O/hpDYY4DEAiPgcAktcFRpBTcfc+tFa7tXx71U+sx88GJrPjqa15//Q2yXL9TLy+XzVH1ORYVF9FA4KsjIjxWxVbhFgSnt0jQwZ0fhSPJGaHsu8SmB5shQyuslh9Yygm5bnzPiryOhi2gHHsg4fEqx5Molr3FkFDosnufatUdcTodSK+JKfVACe4ona6t253g3hUKRXWxe/duEhMT/e8jhQEkJCTQo0cPHn74YVq3bk39+vV54403+Oqrr2jevHnEdX/99VdWrVpFdHQ07777LocOHeK2227j8OHDLFy4EIBVq1axYMEC1q9fX+b+Zs6cic1mY8yYMSd30HKoW7cu77//Pvfee2/Q+Pvvv0/dunUrtWaFhciNN95ITEwMDzzwAEVFRQwbNozMzEyefPJJrrnmmgqttWfPHoYOHcrhw4dJTU2lV69erFmzhtTUVNxuNytXrmTOnDkcPXqU+vXrc9555/Hll1+SlpZW0W3XOhpmZTLriYdY8uFnvP/hu9RtWUR8upMjW+JwF+okp+Zij/JwVpfeNGxUNX7AUHr36EaPLh157pVF/PD1l3Q6sJvNyfXJia8TPtlr8NAKBVqxhmmTGHGR4kUEriI7UfGukLogkYJTg+NEIpZp97mCpDxOc73SG8oWLOEl548nfsI3EorApktchnXN5bJR4nQghIfYYoPYQg/6GQ2VW0ahqAUkJiYGCZGyeOWVV7jhhhto0KABuq7TuXNnhg4dynfffRdxvmlaVZtfe+01kpKSAJg9ezZ//etfeeaZZ/B4PFx33XXMnz+fevXqRVzju+++48knn2TdunVV+g/TUKZOncqNN97I8uXL6d69OwBr165l6dKlzJ8/v1JrCnmi7fQiUFRUREFBQY0WBnl5eSQlJZGbm3tCP0A1kcLCQubOfoHtu38GoCDHQXy6C4cjgSefmvmHFHf7aduv/OOJOTg8TrZGp3PQkQxYH7+G3QpANWOwrCDCspC4k0xvXpZlDfF/TttMtCgDe5SvwJfEbotkiQiO1dAD3ClSloaeaAFzIbJw8BUtE0grYBZCMmZ8qbqBmTKyTCFiyuBBjxleNk1KcBsaLtOG6S3vXlAQjcfQQULmrgLa/HqYgdcO4y/n9g5/iKJG82f43fJnwfe9yJ75CNopcc2U8Nt9D1T4e11YWEheXh4ZGRlcffXVFBQU8J///Cds3siRI1m9enVQfOQvv/xCmzZt2LJlC4WFhXTq1Cnod70vkFXTNDZv3syHH37I2LFjg/5RYxgGmqbRqFEjdu7cWYmTR2bt2rU89dRT/PLLLwC0bt2aMWPG+IVJRalUQTMfsbGxxMbGnswSihMgLi6O+ybdxY8bfuH5ZxcgMvOREm6+efQfVmG2bfOmTHvwAaY8PI3mJfsRHsEhPckSAr5/7dsttwxY8SJ6ibcHTWhZDWGVcXe7dOwOIyQwNBifhcIX02Gly/oXAsD0lnz3mWV80jqyi6VUfJTaMAI7AftqnZT9tYgk3aUUpfuTVj0Ut6HjNnUCl/UYmn/r9Y4VIwX063JW2Q9TKBS1lri4OOLi4jh69Cgff/wxs2bNijjvnHPO4e2336agoID4eCtbccuWLWiaRsOGDRFC+JvK+njggQfIz8/nySefpFGjRlx33XX0798/aM6AAQO47rrrGDVqVJWeq3v37rz22mtVtt4JCZFOnTqdsKln3bp1J7UhRdm0P7M1Tz49i/lPv4Gm2ehw5h9bcyKzfhoP3v93pj4ynWaeHCRwyJbkvy68hch8H+KaU2CUlkAJQ5oaHrfEHuMpw7NhFQrTsGrOGwHjYTNF+DXL2Cf9P7siJEA1sD9OoGPE9L4/XrEzn/gwTXB7e8lYY9ZNhmkpMF+HXre7VJEJU5JyrBhbWl0l5hWKPxkff/wxUkpatmzJtm3bGD9+PK1atfILgokTJ7J3715efvllAIYNG8bDDz/MqFGjmDp1KocOHWL8+PHccMMN/mDVdu2C48iSk5ODxuvWrRsWo2G320lPT6dly6ptGbF9+3YWLlzIr7/+ypw5c0hLS2PJkiVkZWVVKkj2hBzTgwcPZtCgQQwaNIgBAwawfft2oqKi6NOnD3369CE6Oprt27czYEDF00cVFUPTNG65czg33X51tTy/YUZ9Jk28D2mLorknhxSjNJ0tpKo5puZ1yQQRbE4wDR2XS8eQAZaGEIuDlCB9ebCRREiEewKuWPEjwoxwp4a0egvja2In/fso32Pps3x4TKucuyl1zICKqsGWGwu3u1T3J+U6sZmSuo1rb+VfhUIRmdzcXG6//XZatWrFiBEj6NWrFx9//DF2uxW4npOTE1RFPD4+nmXLlnHs2DG6du3K8OHDGThwIE899VR1HaFMVqxYQfv27Vm7di2LFy+moKAAgA0bNjB58uRKrVnhGJEbb7yRjIwMHn744aDxyZMns3v37jIrwVUXyo97ati5aw/THp2JMFxsoQF5ZjwIKE4HI9r6MPbEGphxoT9eEuwmwl46rts82O2mv/S7L1bDV6pd1wzKi+XUvAXIQudYcSQmmrCyYXyxIeH3h6fyasLEFjbfeu8ygn1Nhilwm7aAeyVuQ8Njlioz05TkF5Sahxrszqf1tiP8dfT1XNCtR9mHU9RY1O+WmkNNjRH5M9KjRw+uuuoqxo4dS0JCAhs2bKBp06Z8/fXXXHHFFezZs6fCa1Y4VP/tt99mxIgRYePXXnstixcvrvAGFLWTxlkNue9v4zARNJBH0EzQDMAA4QbhAlEswBV6p4gYL+IPpEBgSg3D1L3xFGV12/OtFi4uLGuFz9tj/YiHBpeWjzf41vQGuXoDWTXhWyd4rdC1recHm4cMj2ad0wSjRCe2wHI0ndkkcjqfQqFQ1ER+/PFHLr/88rDxtLQ0Dh06VKk1KyxEYmJiWL16ddh4eWXXFX9OmjfOJqV+FvGiGDtuK4tFs/rSCAlCCPQSHVEkIns6vFYPKUUZ/VkEbo/thIuU+dwqMsCF47vVLPMZYatYWTPCkjihDfRkBEETOuYr8x64L7fH5i31bk2ILXFjCkG9SubdKxQKRXWQnJxMTk5O2Pj3339fZvn641HhrJm7776b//u//2PdunV069YNsFJ5/vWvf1W6zryi9nLuOWfz0Tu/kSwKOKDVQXMJzCiCwjk0j4YslJix3nReE78IsRCYpkDXI6oV3B4Nhz1y1z3h/xP4i5kE3176DF/ISkRh43MLBdwaod5IRCESlrbr24fwpvDqeDx66X6EJLbEjUiIVfVDFApFreKaa67hvvvu4+2330YIgWmarF69mnHjxkX0lpwIFRYiEyZMoGnTpjz55JO8+uqrgJVDvHDhQoYMGVKpTShqLxf07skH77xNHS2fHL0OthJw+YoNBn6omwKtQENGmcjo8JgM09TQyqhqapoCaRJWsMxbxwyPKdCEVTgs8KGhdgopNUxhooVl6JRRLySkoFqkANRIFVatdGOJ0617LTGW0PKhCYOYEjex9Wtu/R2FQqGIxPTp07n99ttp1KgRhmHQpk0bDMNg2LBhPPDAA5Vas1J1RIYMGaJEhwKA6OhotPh6xOcfxCY9uE07egkY0d7P8QAEAtwa0tCQerAYMQ2BsEd+hkTgMTVsgem33qwVwxuzIcMlQoiw8LlmREDdkYCrkYuiBlVaLbW6BO8tZABTQrFTx2PqaEJ6RUjpvBiXBw2om5oa+cAKhUJRQ3E4HMyfP59JkyaxceNGCgoK6NSpE2eccUal1zypgmYKBUDDhg3Ys+kgUaaJ6QTNbcWouhNK5/hjPIVXjAgQPneLN1A1ktUj8GbTFGje4mdWum9wLEhoefawZfxzBKYMdetEml1aAwQiu2XM4OKsmFJQ7LZ5pY7lDnK7gnOY44rdAGRlnJrS/AqFQnGqSU9Pp7i4mGbNmmGznZyUOCEHdUpKSoWiYbOysvjtt98qvSlF7aJOYhzgTYMFNBOij0HiDtCcXh0RUmMEKTBdWkDdDjDM8n8cfQGgMkSEgKVAAgNEI+uRwBwb4S0RLyKql0jiI1Lmjf96gAgJLqomkEbAuSTEFnsAaJmVHfmgCoVCUUMpKipi9OjRxMbG0rZtW389lDvvvJNHH320UmuekIw5duwYS5Ys8TfjOR6HDx/GMIzjT1T8KYhyWEEhgb1bfH1gTJ1wuRsQQCrdGugSdIk0yzFnYH2ouz3CWy8kKAcYvLYLwwRdKzvl10oVlgH3lU+pFSVyaXfD1ECEixDfyoZLR5ogDW/mkBQk5lo5zW0aNz3u8xUKhaImMXHiRDZs2MDy5cu56KKL/OP9+/dnypQpTJgwocJrnrA9ZeTIkRVeXHF6EBXlAEqFiAQMB5Skgiwj7qMUAYZAmhJDWi3shJDY7CGixmtxEEJDShkUKxIYy2GJEYFdD++kC6Xps0EVYKUIXpPgkvE+94wZZOmwGtpJr0sp1BICIE0w3F5TkM+A44I6eU5kYvxpXxhJoVDUPt577z3efPNNzj777KDWL23btmX79u2VWvOEhIivy59CEQlNs2IgNEykANMORgwIKbAVSjzxHN/4IAXSo2Pq1lSP20S3lQan+guJSRChFU9lcGyJRMPpMdEEVkM9fBYN6S1PFhxiKhEYUkPHtDSKlFY9lIBni4B0XI+pYXgLlgkBxa7Qrn4WHme4OIkx3ES7PaS2aHicL4hCoVDUPH7//XfS0sIz/goLC0+4J10oqoiB4qTZu/cAYLkpTDuWXvDGZ2qmwFYUekcEH4dXZ5R++Gu43ToejxbkEolcdD0c09TxGDqG1PBIgdsUGBIMqfsb1AWvY80p8ThwmToeU7PqfxgahqlR4tZweTSchg0joKeM6e2oG+q2kVJgGuHPqZNfAkCH1m3K2LlCoVDUXLp27cp//vMf/3uf+HjhhRfo0aNy7SpU1ozipNmzczuGFDiLoiHaazkICE7VPAKtRGJGezvehgauehGaDHGbCEwTTNOGpnki1PkI+l/psCwd82XJWPMsdWSYGlKLUMvEq54MqaNLw+/qkVLgMrWwpF8pwWPoCKHh9mjYbaV1UEwjPK0XAWmHi5BA3y5nhX8BFAqFooYzffp0Lr74Yn7++Wc8Hg9PPvkkP//8M19++SUrVqyo1JrKIqI4KfLy8nEWHqHAHYcwBJq3t4wMkbi6E3Cb6GG9Z0oRuhkuUPxuEf8sv/UhUl0PX/ZKUGnX8Cd5G9KVvofADBnhdb0ELhshs0bg7Ydj3eP26FYsS7GGdItSNWQCLg09V1DvaDGOtDRSU1LK2JtCoVDUXHr16sX69evxeDy0b9+eTz75hLS0NL766iu6dOlSqTWVRURxUny6/CuEgHyPlcJrL5KURHktIl4kgCawufWIa/hnlaEbTFPDMEx0PbRIWcRVAtwkElMK9Ail2j2mhi3AKhJaIdWQAr2cbBlrX8Lb3M6HwOW0we4YpEuDFAMMKxhXIKh/LA9dSrp2rdxfVoVCoagJNGvWjPnz51fZekqIKE6Kb75bB0C+2xIiuhuEUZquIvHGi3g/r2VZgkOLIDL8NUYEpml109M1aQkbCYYp0ELKgPjmW1gZLYQGt3rxmBp23fTe5yvnXnqv73qkjru+rJmgzSLwHItCc+loHoHtmMCILX12em4+Eri8b7+I+1EoFIragGEYvPvuu/zyyy8AtGnThkGDBlW6sFml7tq+fTsLFy5k+/btPPnkk6SlpbFkyRKysrJo27ZtpTaiqH2YpsnhA7sp8UThCfDF2IskJb6CYaE/YUJAiQYxwW6YSJkwpaKitO6HYVqxHAIwTB0hJLpmontjSEKLjgXWAQnZSJBVJJLYMNEwpYxYyMxyy+j+PwMYxTpmSWm3YN0lMO0SaQebYVC3oJCo9PokqbRdhUJRS/npp5+47LLL2L9/Py1btgRg5syZpKam8uGHH9KuXbsKr1nhGJEVK1bQvn171q5dyzvvvENBQQEAGzZsYPLkyRXegKL2suGHnxHS43fL+LAVA4a0REjo57smEYYGzmBrghbYeVdGsGz4o08FhqFhmJYIsIqcaThdOiUuPUIZ9vJ9OYapedN6I+MxNX/9EJ8VJK/YQX6xPUigSBM8x6LCHmkrEmBCWl4BmoRuZ6kgVYXiz4xhGEyaNIkmTZoQExNDs2bNePjhh5Fl+XiB5cuXI4QIe+3fv98/Jz8/n7vvvpvs7GxiYmLo2bMn33zzTdA6UkoefPBBMjIyiImJoX///mzdurVKz3fjjTfStm1b9uzZw7p161i3bh27d++mQ4cO3HzzzZVas8JCZMKECTzyyCMsW7YMh8PhH+/Xrx9r1qyp1CYUtZP/LPsfAHnuYCHiiQPTQUQN4AupEG4NXKLUGxJQoCxYhPgIcPVEEBseU8dj2PwBo4FEsmj47nObWki5+PB7PYZOkdNOfkkUxS4HSA2PYcPlLn2WWWQvbc0baOmRAluRoH5uPlLA4D59y9iLQqH4MzBz5kzmzZvH3Llz+eWXX5g5cyazZs3in//853Hv3bx5Mzk5Of5XYL2OG2+8kWXLlvHKK6/w448/cuGFF9K/f3/27t3rnzNr1iyeeuopnn32WdauXUtcXBwDBgygpKSkys63fv16ZsyYQZ06dfxjderUYdq0aXz//feVWrPCQuTHH3/k8ssvDxtPS0urUD8aRe1GSsn2bZtwS50iI9oaw+or44oDtMgf7EGawKmBWeqWkRJM09ufRRL2QnrTYgONJxKkKfzixJQaTrcNwwy2VpSNKLfHjccQFLkceMyAomUChLerrtNls/bsLM1JliHBsQ6nSd2CIqIzMkiIjy9vMwqFopbz5ZdfMmjQIC699FIaN27MX//6Vy688EK+/vrr496blpZGenq6/6VZ/SwoLi5m8eLFzJo1i/POO4/mzZszZcoUmjdvzrx58wDrd/KcOXN44IEHGDRoEB06dODll19m3759vPfee1V2vhYtWnDgwIGw8YMHD9K8efNKrVlhIZKcnExOTk7Y+Pfff0+DBg0qtQlF7eOHnzbhwMlRIwFnokZxiqCkrkZJioY7TkNzElbgw8qeCWw7J6BYt6wgBphOHem2Id16yEvDdGsYRgTLhQz5v3dlt8eG21cMrRzvjJR4rSKl733/N0yre64gXMlYcSCWAHI67aAHzAn5W1W/KBcN6FrJ1DaFQlEzyMvLC3o5nc6wOT179uSzzz5jy5YtgBW2sGrVKi6++OLjrt+xY0cyMjK44IILWL16tX/c4/FgGAbR0dFB82NiYli1ahUAO3bsYP/+/fTv399/PSkpie7du/PVV19V6ryRmDFjBmPGjOHf//43e/bsYc+ePfz73//m7rvvZubMmUFfnxOlwsGq11xzDffddx9vv/221U7dNFm9ejXjxo1jxIgRFV1OUUt576NPADhiJKGb4LGXuk7QQPNYb8woAg0JhHlJJMhCHSls4DDLkMbC26guPHPFZwmxSrMHTzFMDVMKoh1mwHxfbxmBlBqm1w2kixKEEF7xIYKCV4WQpVYV31mE/6lWcz+nbrlmtOAz2kyDxvmH8djsXHV+6S8IhUJR9cTs1dCjqr48luGNaWvUqFHQ+OTJk5kyZUrQ2IQJE8jLy6NVq1bouo5hGEybNo3hw4eXuX5GRgbPPvssXbt2xel08sILL9CnTx/Wrl1L586dSUhIoEePHjz88MO0bt2a+vXr88Ybb/DVV1/5rRC+eJL69esHrV2/fv2gWJOT5S9/+QsAQ4YMKS0W6f0X3MCBA/3vhRAn3Py2wkJk+vTp3H777TRq1AjDMGjTpg2GYTBs2DAeeOCBii6nqIW43W727txCiemgWEYhpJWyK3URFO+hGSCKwYguHYsoNHyNYtyaJUbCxIpEejSk7q1c6ivu4cvj9fagCV/bEhymKXB6G9SVZuEEd5vJLYpBE+CwR6jgitcyEqKHfNswS2zgsqG7BEZUcHpydt5h7KZB378MJCYmppyvqkKhqOns3r07qFllVFRU2Jy33nqL1157jddff522bduyfv167r77bjIzM8tsHtuyZUt/BgpYVpXt27fzxBNP8MorrwDwyiuvcMMNN9CgQQN0Xadz584MHTqU7777ropPWT5ffPFFla9ZYSHicDiYP38+kyZNYuPGjRQUFNCpUyfOOOOMKt+comay9PNVCOnhiJGK71NXc4OhQ0hBUoQErRikwxsHYgKBdc2C9ICVYRJ03Z+5KzBcOrrdQOhWjxciBrWGYwbWmw99qMQbI2JZR5xuOw6bGy2iqLFEhj+XR0gMl4bp0iwNJK3KstKbARTjdtGg4AhaYhJDLh5w3H0qFIqaTWJi4nG7Zo8fP54JEyZwzTXXANC+fXt+++03ZsyYUaEu9t26dfO7XcAqIrZixQoKCwvJy8sjIyODq6++mqZNmwKQnp4OwIEDB8jIyPDfd+DAATp27HjCzz0evXv3rrK1fFTahpWVlcUll1zCkCFDlAg5zfhi+f+QEo4YpX8hhQGY0t/szocE0K14EIFAc0NQPEeojvBoZXeyQ2C4dQxfzIgrPEMmEiUltjLnSbyixvdeCpxue3Cwa8DzTRMMT+kFo9jb5c833cSqpgo0zT2IBowcPtwfdKZQKP7cFBUVhf1913W9wl3s169fHyQofMTFxZGRkcHRo0f5+OOPGTRoEABNmjQhPT2dzz77zD83Ly+PtWvXVroZXSSWLl0aJJCefvppOnbsyLBhwzh69Gil1jwhi8jYsWNPeMHZs2dXaiOK2sHhw0coOLKPPDMOT8CPj8CyinhiCfKW+GuJ+OcJhEsiowJuDCTQKiIJ96IgMD3CF72BYWroDg9C8/pJQtYzPAJparjd4HCU5a8MTwd2uW3YdQNdN/0pxR6XjuFN2bVHeXAX28LuFQiEAcklhaSWFJDQKJuzO55ZxnMVCsWfjYEDBzJt2jR/gc/vv/+e2bNnc8MNN/jnTJw4kb179/Lyyy8DMGfOHJo0aULbtm0pKSnhhRde4PPPP+eTTz7x3/Pxxx8jpaRly5Zs27aN8ePH06pVK0aNGgWAEIK7776bRx55hDPOOIMmTZowadIkMjMzGTx4cJWdb/z48cycOROwsmjHjh3LvffeyxdffMHYsWNZuHBhhdc8ISESmhu8bt06PB6P36e1ZcsWdF2vdMMbRe3hpUUfIgi2hvjQPL4S7iJYhPj6uXjLu2umwPRITFsZ5d49GmhmgOslxJxhCKTwhpZIgeG0odkNNEfwPJ8IATAMgWEIdC34mabf8hGqeARuw4bbY2K6NTxum38cwFVoB48ecf8CSbPcg5gIxtxwfYQDKhSKPyv//Oc/mTRpErfddhsHDx4kMzOTW265hQcffNA/Jycnh127dvnfu1wu7r33Xvbu3UtsbCwdOnTg008/pW/f0rpDubm5TJw4kT179pCSksKVV17JtGnTsNvt/jl/+9vfKCws5Oabb+bYsWP06tWLpUuXhmXbnAw7duygTZs2ACxevJiBAwcyffp01q1bxyWXXFKpNYUsr9xbBGbPns3y5ct56aWX/AVNjh49yqhRozj33HO59957K7WRU0VeXh5JSUnk5uYe17enKJ9fd+xi+owZSAQbnU2RIZ49CbjiwZUskF6PRej1wN4rnigTwmO9LGyG13EoQJRm00iJFdRqC+9NY4txoXldQ4EixNdQTwhJVFRpMKovRdcX1BqK4dbwOPXIjfaKdaQB2AETtGIN3WW5lc44fICs3KMYWU144f4JZRxQUdtRv1tqDr7vRasx09Gjqu5D14fhLGHTU39X32sgJSWFVatW0aZNG3r16sWIESO4+eab2blzJ23atKGoqKjCa1Y4WPXxxx/nk08+Cauq9sgjj3DhhRfWOCGiqBry8vOZNWs2Gia7jmRgxnvdIwEyVgD2InDWi9jcJWxIM4lQoQNLIZRoEBPBYuKxhIOIEEhiOHUMKdCijAARYu1MeNNdPB4du92wUouDKqoGW0QMl4an2PrrIUO7/noEmMIaK5LYi20It0RqUD8/j6zco9iT6/LMhPGRTqdQKBS1ll69ejF27FjOOeccvv76a958803A8ow0bNiwUmtWOIIuLy+P33//PWz8999/Jz8/v1KbUNRsDMNgypRZSLOYg0frUVIYa32Q27wvb1KKxBIX+glWE5YARsAbA3B5y79Lbz+akCqqyAiuGu8106VjerQQERKMx2O5aKRvvUjndfpEiFeoeLxxK775bu/6Jti9YkXzCJJynbQ6moNps/Pg+HvRdT3C6gqFQlF7mTt3LjabjX//+9/MmzfPX8h0yZIlXHTRRZVas8JC5PLLL2fUqFG88847/qpqixcvZvTo0VxxxRUVWmvKlClhTX5atWoFwJEjR7jzzjtp2bIlMTExZGVlMWbMGHJzcyu6ZcVJIKXk0ZlzKSw4yLH8RI4VWGZJzeP9VBaAbgkSdyyUJHvrhhxvXUCYAs0pwAO4QRii1NYhANMrAnx4AiwYoZYS0zcoI+iU0hLy0tApKYiipNARwR0j8ZToeEpC/UoC6RZQIhBFpcLEllf618fmMmhdsBcNyagbRpNWr+7xvwgKhUJRy8jKyuKjjz5iw4YNjB492j/+xBNP8NRTT1VqzQq7Zp599lnGjRvHsGHDcLvd1iI2G6NHj+axxx6r8Abatm3Lp59+Wrohm7Wlffv2sW/fPv7xj3/Qpk0bfvvtN2699Vb27dvHv//97wo/R1E5Xnn1HXb99jOFJTEcOFoPXyVRzVNqzABvsksUIATCTWlWTMgcKK2+Kry1OYQhQj/3S3ELsPlql0XWzb4+NOE3lw5ZIkRDGtYapqHjKhJouoluNxGaxPV7tCVOomTY/UiBnmtDeC0yhl0ivIGuwiNpUbCfWMNFs7N60rNLp4j7VCgUCkU4FRYisbGxPPPMMzz22GNs374dsAqtxMXFHefOMjZgs/kLsQTSrl07Fi9e7H/frFkzpk2bxrXXXovH4/ELFsWpY+vWnaxa+Qluj519h+rjL3lOgEXEi+EVIQC6r7hZQJ84C+kXIaWGDeEbLQMBTh3sIcEoQQ8PHjQNgaYHxJdIkGapCAlcyDR0TEO3rDIeDSG9xddCniFMy1rje55uCjSnVT02s+Qoqa58bKkZjB+t2hwoFApFRaj0p3lcXBwdOnQ46Q1s3bqVzMxMoqOj6dGjBzNmzCArKyviXF/EcnkixOl0BjUiqkjjHUUwH360DCFg/5F6mGZwvIPuLhUGpgamPbhaqfDGVATWEZGCIBESOJ9Aa0loto1hBYeGWSrwddYNLR4iQJcBc0QEESKD53vjUYTPDWMvI7VYWHPsJQbCLUjwlNC04CAeRwyz/jbO33tBoVAoFCdGhYVI3759y/1l+/nnn5/wWt27d+fFF1+kZcuW5OTkMHXqVM4991w2btxIQkJC0NxDhw7x8MMPc/PNN5e75owZM5g6deoJ70ERGSklW7ZsxO3RKXLGhJf88mD1gBECIzpCni7eviyekHoiYVMjuEHC1rF8K9LbVK50j1gCJex+genR0OzeZndhIiQAEytDh1KLh+bWkBhIR+R9RrvdnL1/u38rBoLxd48hKTG+7OcoFAqFIiIVFiKhNevdbjfr169n48aNFaqjDwS1Re7QoQPdu3cnOzubt956KygIJi8vj0svvZQ2bdqEdToMZeLEiUGVYPPy8sI6JiqOz/oNPyPNEvILkwDLzSICgkKE1+phOggr6x4U7iFBGiBtMrzzLtbHP1KGh3hIr5Dxj1uxJ6ZDWmIkKC4kAoawmtzZpOV2ifSTLoESr/9IWmLDX7PVpVkl631N7ALOKDVrC6Y9CjM2gaGXD6ZVs6blbEahUCgUZVFhIfLEE09EHJ8yZQoFBQUntZnk5GRatGjBtm3b/GP5+flcdNFFJCQk8O677wZVkYtEVFRUxI6Iiorx3/9a/QryC+MRJrjiBPZiablcvAiPxEgMtzb4BITEEinSG2xqlvXT5psfkJkrIokMCcIjkLppTQu1hvgwAadmFVzTpFWp1SGDO/tKrNiTAHUkgiUUmkcgTQmmJPqgjiteIjSBU9hwoxEdE8czs6aVcSiFQqH4c1CRjNh33nmnwutXWSeua6+9ln/9618ntUZBQQHbt2/3N/rJy8vjwgsvxOFw8MEHH1RpmVpF2RiGwa7fNuNy23G6HP4QVXeMFmTVEGYZiSyytM+M9IsPYblzDMItGbL0/yLUEmLdWvoHA0tolNE/ShhYrhYprJfPLeMSUKSX1gNxaWWuEbgtYQochzU0j0D3phrrJYJionDlHatwIyuFQqGobSQlJflfiYmJfPbZZ3z77bf+69999x2fffYZSUlJlVq/ylJPvvrqqwoLhXHjxjFw4ECys7PZt28fkydPRtd1hg4d6hchRUVFvPrqq+Tl5fkDT1NTU1WxqFPI6i+/BTzkFdTBnyljSqRNwxOjYSsyrXAPM7JvxApeJchH4xMzwsSyMgS0bvEVFzuROE/hzaKRUUbkoFdXWT8X3qJkhd4MnFBrSiTx4xUtutOaay+y3mtuKBLRJMpiDh78nfT0+sffuEKhUNRSAhvZ3XfffQwZMoRnn33W/zlsGAa33XZbpcvfV1iIhJpopJTk5OTw7bffMmnSpAqttWfPHoYOHcrhw4dJTU2lV69erFmzhtTUVJYvX87atWsBaN68edB9O3bsoHHjxhXduuIEWbZsBWC5ZcBrGTC8Fg5dWGKk2ETzQEgyDRLwxBPZ1haYsuuxYj2kLiFy77gyEVIgPVppSq8ErcSycISnDEfAKKMoWuhBTLDnBqTtWsVMsBdInB4HOGDdD79wiRIiCoXiNOFf//oXq1atCjIG6LrO2LFj6dmzZ6XqiVVYiCQmJgZlzWiaRsuWLXnooYe48MILK7TWokWLyrzWp08fKtiPT1EFHDp0hP05O3C5o3B7rHicUOuHtAmK6mkYsSL8w7y8D/cIFgxpYP0UlnVfhHGJtEwongCBYIBmaJjCDAueDVvQkJYLx+Z1LQkQLrAVCTzxVhyM5hQIp0RzB/yse0xaFOaAR8dpWl+bTVt2cUnFfuwVCoWi1uLxeNi0aRMtW7YMGt+0aVOlXdUVFiIvvvhipR6kqPksXbqcd955G00zOZYfnD6tmTKoLYyQAld8uEoQEvRiMGIJEhGlYSAyKJVX4A0IjVRfpAyEFAiPbokICZoh/PEewi2QQgZbZAKCYAEwhCUw3N79aKA5wV7kjWMR3pRht8BWbAXZCmHSumQvdSgEb1qvKQXXXX3piW1aoVAo/gSMGjWK0aNHs337drp16wbA2rVrefTRRxk1alSl1qywEGnatCnffPMNdesG99I4duwYnTt35tdff63URhTVQ2FhISuXr+WDDz/FkIcxTY0Dh1LJL7aEiE8bBKbuSh08cdaHuekgsoDw1fwIFCNaQJGwwPgR6RUjJxL2U1agq1dsCASaL803MEPGVzVNWtVQfW+FN3bElw1kcwsMb3yLvchq4ochaeXcSx2jkIzmrTmr/Zm88977ZDdrTWqq6imjUChOH/7xj3+Qnp7O448/Tk5ODgAZGRmMHz+ee++9t1JrVliI7Ny5E8MwwsadTid79+6t1CYUfxyGYfD9dz+ycvkaduzcToknzx8kWlgQw4HD9TAMG9KbJe0TIMIbUSp1gStBWHU9wn8MLDSvLggUI6FWihD8YqS8PK4IIkSLlMLrS/P1xZAEWkQECDfhoscXbuKtjyIM0L3nyyo5RIqnkPgGjXnw3jFomsalF/UtZ6MKhULx50TTNP72t7/xt7/9zZ9AUtkgVR8nLEQ++OAD/58//vjjoDQdwzD47LPPVABpLeDOm/6GYbfqvZgeQckxO0WFMRSaCbjcNtBL4z6MKA3pkWguaQ2Z4EoWoJX2nAksz+5760vx9YkRqcvjxG148bkXI1VhLaOuiJ+gDB0rDkTqhAsb00rJlYHWGukNxvUNuSGq0CTBKCbBLCbTcxgZFc+s+/+GplVZxrtCoVDUak5WgPg44d+qgwcPZvDgwQghGDlypP/94MGDueaaa1i2bBmPP/54lWxKcWr4/tsfMOwFFB+ysf/7RPZ8mcyhnxPILUzC7bb73RTgtYAIgWnXcMVrmDq4o62smSBCKpwKCPqp8gkYv3ekHAKfH/aMCES0hgSspbtE8L3S6hosEGFragHZwGfk7ues4i20ce2mkecQEo0bRl2vUsYVCsVxady4MUKIsNftt98ecf6LL74YNje0FMaUKVNo1aoVcXFx1KlTh/79+/uzSgP5z3/+Q/fu3YmJiaFOnToMHjy4ys934MABrrvuOjIzM7HZbOi6HvSqDCdsEfFFwzZp0oRvvvmGevXqVeqBiurj7UXvAXBsWxRup24JDZuGtJX+8AhppWT7ioqZOrgSNFzJOmaIbJVe44kMtYoEZrsKQBdoJhie4PohkRAIpOGNFwkNMg18tve/gdVQpQgthhYQL+J9pvBm2ghZum+rb07pOnU8hZhodOjek55dOtGhTcvjVvRVKBQKgG+++SYofGHjxo1ccMEFXHXVVWXek5iYyObNm/3vQ/u5tWjRgrlz59K0aVOKi4t54oknuPDCC9m2bRupqakALF68mJtuuonp06fTr18/PB4PGzdurOLTwfXXX8+uXbuYNGkSGRkZVdLos8IxIjt27Djphyr+eHL2HeBQ3l5Kjui4C3X/d96MtodXE/N1zgXcccJ/TcgydIHPKiICxIf3/4F1PTQpMGWo+SQcqwOuV7RISuNMAgRPJLdQxLVMgTCkP8NGeAL27NuX2xsLowEGeNCx6YI7b7iu/MUVCsVpRWg390gtRXzCwMejjz5Ks2bN6N27d5nrCiFIT08v8/qwYcOC3s+ePZsFCxbwww8/cP755+PxeLjrrrt47LHHgvq0tWnT5rhnqiirVq1i5cqVYX3nToYTEiJPPfUUN998M9HR0Tz11FPlzh0zZkyVbExRtby68C2EgPxd3r803s65Zqw9vKSpBEyJO0YgteAAEFNItJDudUFWEe9LagTHYeANAvUO+GvEhAgM33OEryQ83kqsvnUD5ukl3qwdf0AsEZWScAukQ1r1Q3x+IxPsuVYFWOGx6ojoHm8gq/f8CoWidhGXY2KzV33bBY/bWjO0gerkyZPLbcTqcrl49dVXGTt2bLmWg4KCArKzszFNk86dOzN9+nTatm1b5prPP/88SUlJnHnmmQCsW7eOvXv3omkanTp1Yv/+/XTs2JHHHnuMdu3aVfC05dOoUaMqr/F1QkLkiSeeYPjw4URHR5fZ9A4sVaeESM2juLiYrTt+xlOiUXLE64aRIB06RAi+FIChCczocMEBwS6QwF4zhsOquRExSNRHqOgow70TNCXA1yOwRINvTHdZ+zHtZWfc+FJ6ZUDars0pcBSCYZMIj+U6EgboHkkMLmISVVquQqEIZvfu3UEBmsdrsPree+9x7Ngxrr/++jLntGzZkn/961906NCB3Nxc/vGPf9CzZ09++uknGjZs6J/30Ucfcc0111BUVERGRgbLli3zh0j4ymZMmTKF2bNn07hxYx5//HH69OnDli1bSElJOYlTBzNnzhwmTJjAc889V2UJKickRALdMco1U/t469X3EbpJ/u5o/B/xUmLElB33oHkkmGaYUBHSKzR8De98BgYdjChKf6JkyP99fzRLy7r7REgY5YltGSBMAvakuwRSkwFN9gLvsRbUnN5gVQ/YC61LuksExYvEU4SOSUxCcjmbUCgUpyOJiYkVyhRZsGABF198MZmZmWXO6dGjBz169PC/79mzJ61bt+a5557j4Ycf9o/37duX9evXc+jQIebPn8+QIUNYu3YtaWlp/hjO+++/nyuvvBKw+sM0bNiQt99+m1tuuaWiRy2Tq6++mqKiIpo1a0ZsbGxY/NyRI0cqvGaFY0Qeeughxo0bR2xsbNB4cXExjz32GA8++GCFN6E4dZimyZdrVmFIQWFO6Q+MFBIZfZwATF8dkACEtCwQUvc2tqO0yV1QDEmI1cN6qDdWI0I/GL9XRZYhTvwbj/DeFzdiegNdA2JKhAc0jwAX2F1W2q7mDHiGBCFNMjxHSTILSZTFGGhc8ZeLy9uFQqFQlMtvv/3Gp59+yjvvvFOh++x2O506dWLbtm1B43FxcTRv3pzmzZtz9tlnc8YZZ7BgwQImTpzo71gfGBMSFRVF06ZN2bVr18kfJoA5c+ZU6XpQCSEydepUbr311jAhUlRUxNSpU5UQqWF89t8vwOamYIcDGZDuakY7wmND8FZNjdKQNiuNLPRz31eBVILVtE5QWl3VG+QaUYRgue40U0bM0C1dtJxr5WTRBO5PSstlIwJSdX3xJsLjrZbqJcZ0coZ7H3HSiYmA6ESGD72asztF9s8qFArFibBw4ULS0tK49NKKtYEwDIMff/yRSy65pNx5pmnidDoB6NKlC1FRUWzevJlevXoB4Ha72blzJ9nZ2ZU7QBmMHDmySteDSggRKWXEoJsNGzZUqR9KUTX89z9LkaYk/9fg8qJmXKlvUwLSLjAclgAJuhCKifVJr4lgEULkWI9AJNLq1nsC2S5lIY4ThyZMge4MX16YWI3x3KUbTfccJdv1OwJJ6y49uOOGa7HZKvxXQqFQKIIwTZOFCxcycuTIsN8pI0aMoEGDBsyYMQOwvAxnn302zZs359ixYzz22GP89ttv3HjjjYDVhmPatGlcdtllZGRkcOjQIZ5++mn27t3rTwlOTEzk1ltvZfLkyTRq1Ijs7Gx/F9zy0oZPlLy8PL9LKjRzKJTKFDk74d+6derU8RdbadGiRZAYMQyDgoICbr311gpvQHHq+OXHzRSZeRTuMjFynYh4m1U7RBPIKJslDBwaRpTmr5YaRFlD0ooJMR3BWSzgtUiEVkX1YgYGsUYQI/6U3LKQltAo97onoJaIb6pp/dFWbAWkmtGQ7jlKE9dBpD2aW265mS7tlQVEoVBUDZ9++im7du3ihhtuCLu2a9euoArNR48e5aabbmL//v3UqVOHLl268OWXX/rdLLqus2nTJl566SUOHTpE3bp1Oeuss1i5cmVQZs1jjz2GzWbjuuuuo7i4mO7du/P5559Tp06dkz5PnTp1yMnJIS0tjeTk5IjGCJ+RIlILmOMh5Anm4bz00ktIKbnhhhuYM2dOUIl3h8NB48aNgwJuagp5eXkkJSWRm5tbZeVoawuTxz/M/rw95Cx14yrQEYkJlkiIjcYd78CVGg3lVMJzJmmY9vAqpIYDnMnCihEJzfwVIB0hY8jSdN7AmJPQLJfj/SRKrK65kfAafLQS/NVchQm2IjCiLQFiLyj17sRTRBtjFy7TzlOzZ5CYkBB5XYWiDE7n3y01Dd/3ostVj2CzRx//hgricZfw3dsPnLbf6xUrVnDOOedgs9lYvnx5uanI5dVLKYsTtoj4/EJNmjShZ8+eqtJkDae4uJicY3vw5EtcRyXopvU5HxMFmoa9yIPHDaYWUswskEDLQmD2i88lU9Y9oZVWA4NTZRl/Dry3ku1chMvKitE9gAGYoLsBj0QLqMEqDEmRiGWPJ41GtoNMfeQf/OPRKVVSIVChUCj+bPTu3ZsdO3bQpEkT+vTpU+XrV9ghHqh2SkpKcLlcQddPR7VYE3E4HMTq8RQnFZDYRiPvZxOi7P50XIkVtCmcEsOBJUbK+xwWVvyH1L3FykIFRCAm/u62pubNYvEWKZOBNdhDhIjmsqwYRhRlihEZUtbdj8eyeGgRrCqaDI4t0VwSmwcOkUQShXBsP0ePHlUxTgqFQlEGzZo1Izs7m759+9KvXz/69OkTVOfkZKjwvz2Lioq44447SEtL8zfgCXwpaga6rvPwrCnYjWjqnGkjoV0UOAKqqNo0EFYNDd0pEZEqicrS//lSdoOsJwLrJ0gPefluFzLofeCaEd8Lq0aIr+S6kBGum8H32golsfskjnwZLkK8Qak2J9gLTX89Ed20zq0ZAl0YGKZGcnJy+PkVCoVCAcDnn3/OyJEj+fXXX7npppvIzs7mjDPO4JZbbmHRokUcOHCg0mtXWIiMHz+ezz//nHnz5hEVFcULL7zA1KlTyczM5OWXX670RhRVT0JiAg/NnAJOGyntJXGpJf5rgV10BZaVQHOXflh7Z1kCRPOWeveKECNaWO6PspCAr79L0FPKmOv7o3e+Zgg0t0BzCTSXZSnRnKAXglZsxYFoTu9YEdiLwZ5f1toCzSXRPaAXS/QSiTDAJjyk2Q4Ro5XgJDYoeEyhUCgUwfTp04cpU6awfPlyjh49yrJlyxg6dCi//PIL119/PZmZmWWWpT8eFXbNfPjhh7z88sv06dOHUaNGce6559K8eXOys7N57bXXGD58eKU2ojg1pNStw98fnMAjD8+gbrN8pCEoOhKF1EMqpgKaYRJjdxId48Ie4yLHlkIRMcHxHoDpzfzVPFZJ90jZL4ad4AyZ8nJ2I1Rh9U331QHRSrDSct3Wur4Kqj5BZHNZZd6lvfQZmlsiTGufALpLkkghKTG5JOkFCAEeaeP88/uX8xVUKBQKRSDR0dH069ePXr160bdvX5YsWcJzzz3Hpk2bKrVehYXIkSNHaNq0KWDFg/jKufbq1Yv/+7//q9QmFKeW7KaNuH/yRGZMe5R6Z+RxcHMiBR47UQ4n0Q4n0VHWy2F3BXlebHjYTHAxHBkgPIQ39kKGuF9Mb0xIMKX9YgLf+sdCAmIDy66Dt/6H9z7NDUaE6qz2InAlSH8qsu4CpMQuPKTY8kix5xKluZESDFsy1/z1L/Tr3RO9nMwhhUKhUFi4XC7WrFnDF198wfLly1m7di2NGjXivPPOY+7cuZXKmIFKCJGmTZuyY8cOsrKyaNWqFW+99RbdunXjww8/VH72Gkzjxo24996x/OPxf5DWMo9U8tECgiqEiCK5TiMaN25C585teXfJUti7gxhKKKY0HS60IJkwvIGrIWm5AmF1aIxkBIk05quCGogv6NUIjg0RWGLEdBAmXmzF4Im1NhjtcZGhHSYlNg/Na/1Iy2rL/40aQmZm2S23FQqFQhFMv379WLt2LU2aNKF3797ccsstvP766/7y8idDhYXIqFGj2LBhA71792bChAkMHDiQuXPn4na7mT179klvSHHqaNGqGbfdejtz584nLjGe7MaNad++FV27tKdOneSgudGx0Tz91BPU5yg7sX7QpAZCK+03A5Yo0D1eV4zwCpWAQmJhwaqRkFYciL8HTYiLRvOEO3Y0w+rJF4ruBrvLRQPPEVL1XIQAl57AlYMu5aLzz1WVUxUKhaISrFy5koyMDH/GTO/evalbt2q6lFf4t/I999zj/3P//v3ZtGkT3333Hc2bN6dDhw5VsinFqaNjl3a8sPDJ487r0LYV9oQU6uYfZQ+peLCVumBMQgJRrYqmgeXegdK6IOXhFSGaFEiPJXD89T68zxLugGUD1rdcL8GLpYsjZDgPoQmISqrHsCF/5eyzOh33vAqFQqEom2PHjrFy5UqWL1/OzJkzGTp0KC1atKB3795+YZKamlqptU86VSA7O5srrriClJQUbr755pNdTlGDuOD889GQpJJraQoNMLyxG1rwS/i63ZqWu8ZypwjLKuJrhke4+0VzW1ky1rUIKbvlIGRpOXcbHppre2igHcKjRzN05GienPWIEiEKhUJRBcTFxXHRRRfx6KOPsnbtWg4dOsSsWbOIjY1l1qxZNGzYkHbt2lVq7SrLWTx8+DALFiyoquUUNYC/XNgXqTtI4ygCaXWuNaxXRLEQUmJESBBFYM8T2HPBnge2AtALvD1hPF7REkikdXVf5GrwM8C6P54i2mi/kSSKyGjSimfnzKJPz26qUqpCoVCcIuLi4khJSSElJYU6depgs9n45ZdfKrWWcpgrykTXdaLrNUIc2E6sKKGQGOtCWRaLkEAO4QGb22cqEf7+L5oE6a0P4okXYXJYc+K3ruhOiWZYa7kTQtaXJg3kYTK1w0ihMWDglVzxlwur7PwKhUKhsDBNk2+//Zbly5fzxRdfsHr1agoLC2nQoAF9+/bl6aefpm/fvpVaWwkRRbkk2AVOwCkCegtFMDRIvEGsvikebwCqAAyv2PAHfwjsRRLNbQW6uhK9acHeeBHh7aLryLXiQHy3GVHSasIHxJvFNHXnEIsLpx7DpPvuoVnj4FRjhUKhUFQNycnJFBYWkp6eTt++fXniiSfo06cPzZo1O+m1lRBRlMvRQwdxSRseb/qLP1YkQhEzv/fE7W02571Bd1kdcIOQYHOC225VdTVsAnueZf0wdbAXegVJwC1RueCsY5DtPkS6cRSJoG23c7ht5DCVDaNQKBSnkMcee4y+ffvSokWLKl/7hH97X3HFFeVeP3bs2MnuRVHD8Hg8eEryKZKxVqaMDUxRToVUQkQI1lTdbRU983fs9QaZCjdgSvQSsBdI7MWlQkdGKIqWQCFtCw4QjRt7Ygp3/d8tnNG0cVUeWaFQKBQRuOWWW07Z2icsRJKSko57fcSIESe9IUXNwOVy8fzLbyKQFGrRYPPmu+gicoyIt96HCKnt4Wte5y9A5sM7T3cCMQKbU/pDSXxpv761dAwa6L9TT8/FkIKz+vTnhmuuVP1hFAqF4k/ACQuRhQsXnsp9KGoIbrebN975iFUrPkcYLtxS54hILJ1QljFEhosQf6EyvP1eTPyBqT63i+4EI1qW1hwRpfOEhCSRT5btIHbhITa5PmPH3EqjBplVdl6FQqGoavbu3ct9993HkiVLKCoqonnz5ixcuJCuXbse997Vq1fTu3dv2rVrx/r16/3jhmEwZcoUXn31Vfbv309mZibXX389DzzwgD9DcMqUKSxatIjdu3fjcDjo0qUL06ZNo3v37qfqqFWCcqwrAOuH/M33/svyzz9FeErwSJ39MpXfZTJSWOrBbwjxAO4ImiR0QJQ2nBN4Y0W8DfN89T98ganCBFuJxBMrEAJi3CVk6EdI0fMxhU6/AYMYMvhilZKrUChqNEePHuWcc87xN4NLTU1l69at1KlT57j3Hjt2jBEjRnD++edz4MCBoGszZ85k3rx5vPTSS7Rt25Zvv/2WUaNGkZSUxJgxYwBo0aIFc+fOpWnTphQXF/PEE09w4YUXsm3btkoXG/sjqFYhMmXKFKZOnRo01rJlS38Hv+eff57XX3+ddevWkZ+fz9GjR1U/m1NAzv6DTJj8MNG4MKTOflmP32UdTJ/5QgKmRGhW2VMdK8XWDEikkRCcGYPlltECrCSaOyBoNWBcd4JuGiRSQpRZTIK9iARbMQBRyRk8+Lc7qVdFpYQVCoXiVDJz5kwaNWoU5EVo0qTJCd176623MmzYMHRd57333gu69uWXXzJo0CAuvfRSABo3bswbb7zB119/7Z8zbNiwoHtmz57NggUL+OGHHzj//PMreaJTT7U72du2bUtOTo7/tWrVKv+1oqIiLrroIv7+979X4w7/3Bw6fIQpj8wgGhf7zLr8aDblgKxbKkK8CDOsmnpYI7rACqrIgI65XjTTsn4A2KWHJC2fTNtBWovfONOxjSaJe8h0HCaOEhJSsrj51jt4auYUJUIUCkWNIC8vL+jldDrD5nzwwQd07dqVq666irS0NDp16sT8+fOPu/bChQv59ddfmTx5csTrPXv25LPPPmPLli0AbNiwgVWrVnHxxRdHnO9yuXj++edJSkrizDPPrMAp/3iq3TVjs9lIT4/cCfXuu+8GYPny5X/chk4j8vLzmfTQdHAXsVPU56CebAWWhlY7xVsB1VaaMSMkiNC03MAmdyJYiAhpkiiLSCwoIIkiYnGBN3jVIzViktJp1eIMup/VkfZtW2K3B5hbFAqF4gRI2FmMzXaCPSIqgMdTAkCjRo2CxidPnsyUKVOCxn799VfmzZvH2LFj+fvf/84333zDmDFjcDgcjBw5MuL6W7duZcKECaxcubLMUgQTJkwgLy+PVq1aoes6hmEwbdo0hg8fHjTvo48+4pprrqGoqIiMjAyWLVtGvXr1KnnyP4ZqFyJbt24lMzOT6OhoevTowYwZM8jKyqr0ek6nM0il5uXlVcU2/3QUFRXz98nTMEvy2WPW46Dd8l9KwNClVQE1oJS7CKiy7kMzwTTwN8PzWUWsnjSSGMNNsiwgSRaSKIvQvIu5hJ2E+o05s11renXrTHZWQ5UBo1Aoajy7d+8mMbE0eD8qKipsjmmadO3alenTpwPQqVMnNm7cyLPPPhtRiBiGwbBhw5g6dWq5NTreeustXnvtNV5//XXatm3L+vXrufvuu8nMzAxat2/fvqxfv55Dhw4xf/58hgwZwtq1a0lLSzuZo59SqlWIdO/enRdffJGWLVuSk5PD1KlTOffcc9m4cSMJCQmVWnPGjBlhcSeKYJxOJxOnTMddeJT9Rgo5WojrQwikDoYm/S4W3Q12Ca6Qv3e6CzzRgABhShI8xdQz80hxFhItLZOIRBCdnErr1m24sHdPmjbOUkGnCoWi1pGYmBgkRCKRkZFBmzZtgsZat27N4sWLI87Pz8/n22+/5fvvv+eOO+4ALDEjpcRms/HJJ5/Qr18/xo8fz4QJE7jmmmsAaN++Pb/99hszZswIEiJxcXE0b96c5s2bc/bZZ3PGGWewYMECJk6ceDJHP6VUqxAJ9G116NCB7t27k52dzVtvvcXo0aMrtebEiRMZO3as/31eXl6YOe10xuPxcP9DMynJPchBI5l9Zr2yfwqEFZwqkAi8dUI8EmkrFRFCgs1j0NBzmFRXHlHSCgIx9CjSG7emZ7eunNe9CzExMaf+cAqFQlHNnHPOOWzevDlobMuWLWRnR25BkZiYyI8//hg09swzz/D555/z73//2x/oWlRUFGY51nUd0wytmxCMaZoRY1lqEtXumgkkOTmZFi1asG3btkqvERUVFdFcprB+ICdN+wf5h/Zy2Exkj5kGCGSodSKweKop/bVABOAoBGdi6YQY08kZRXuJkS48tiiatu3IoAv70bJZU2X1UCgUpx333HMPPXv2ZPr06QwZMoSvv/6a559/nueff94/Z+LEiezdu5eXX34ZTdNo165d0BppaWlER0cHjQ8cOJBp06aRlZVF27Zt+f7775k9ezY33HADAIWFhUybNo3LLruMjIwMDh06xNNPP83evXu56qqr/pjDV5IaJUQKCgrYvn071113XXVv5U+DlJJNW7bz0dLP2bxpI7rp5KgZz29GOpatg/Iqtgc1nQPLKqK7wHBAXSOPJkYOGpJuvftzw9C/KvGhUChOa8466yzeffddJk6cyEMPPUSTJk2YM2dOUFBpTk4Ou3btqtC6//znP5k0aRK33XYbBw8eJDMzk1tuuYUHH3wQsKwjmzZt4qWXXuLQoUPUrVuXs846i5UrV9K2bdsqPWNVI6QMDUH84xg3bhwDBw4kOzubffv2MXnyZNavX8/PP/9Mamoq+/fvZ//+/Xz77bfcdNNN/O9//yMhIYGsrCxSUlJO6Bl5eXkkJSWRm5t7XN/enwUpJT/9soX/frKc7Vt/AY9Vk8MldY6YieSY9ZDe9FwpwIiKIB4EYEqr+VzYNUn9uINkGMcwdQejb7yRszvX7PQwhaKqOR1/t9RUfN+LPmfdj80W2mHz5PF4Slj+zTT1vT5FVKtFZM+ePQwdOpTDhw+TmppKr169WLNmjb8C3LPPPhsUeHreeecBVr719ddfXx1brrGYpslX32zgf6vXsHPbJjCsdDO3aeOYJ5lcTwJ5egymXQRVCLGay/mavARjc0Y2lmTbcqhj5ONIqsfk+8aqOh8KhUKhqDTVKkQWLVpU7vUpU6aE5WgrSsnLy+fzz7/k22+/5+DvuxFYgaIu08YxTx1yPfEUmd6UFgCbACEwdemvFSIFEUWI8Eh/8bFQ4rRi3CKaZ2Y8hK7rkScpFAqFQnEC1KgYEcXxkVLy0NR5HDy0E7cr168hnG4Hx0QKx2Q8RUa0leESeq9vIECMyDJCOoJiQ7yxqZrLxFZkoqeZxMfEKxGiUCgUipNGCZFaxisv/Zt9ORswDEFRcSyF3ldxQhTuKGGl1gowbaB5ZNlxqEJg2GTkIv+mDCjVLtFLJLYiw29FEUKil1H9T6FQKBSKiqA+TWoRP23cwqovP8PttrF7dwMMqYEGpk1g2kRpPxhhvUy7QHMHiJEA0SEBqYvIQSCawHBIHIcMbMWmv1Ou705Nk9hsqgS7QqFQKE4eVVe7FjH/+YWAZP+BNEypWRpCghGtWU3oJMF12L1iRAYO4BMhlJu2iwl6mAixrCEAdiVEFAqFQlEFKCFSi9A0HSnB5QoWAUa0rwIqiNCgjwAxIsWJihCJkGBGSOv1CxHVlE6hUCgUVYASIrWItm3bomkQG2PVBTF1cCXpwVkvZoSyMD4xIuXxRQjeUu4SDLtAasEBrZpXiDgcjpM8jUKhUCgUSojUKs6/oBcAMfEluON0PPE2r5kjwPlSRnk6KUCeSESQ1xriEy8IAZolSACEsPoaqGBVhUKhUFQFSojUIrKzG+Ix7MTFFflFhQAI6Hkkyut/ZFK+NURKtIDaIaYjYLKwxIjw/sQcO1xUob0rFAqFQhEJJURqEUII0jKycegeHLob8GXJCL/7xIoTiXy/ZgLlFfQ3g++VugiuMyIEptcNVDdZddNVKBQKxcmjhEgtY8AF5wKQ4CgELOEhdfyuEyBinIgQ1jy/UAmdEmIN8S9lF0FzbEesVrxul7vSZ1AoFAqFwocSIrWMnmd3RiKI9woRiVeECIGpe4ugRnDPGDaCglqF/2bvy4xgSZFg+JJjTInjUDFasbW4y+2qukMpFAqF4rRFCZFahsPhID4pnTh7MZowrO+gT2D44jhChIipAVqEVNyAl71EWpVYffdKq6S7cBkIwyTq9yJ0l4k0rXXcbmURUSgUCsXJo4RILaR9+3YIAXH24mCXDJQqi4BMmuNlywhDIqRAM8DmlNi8oiR2XwlxOS6ifstFc1kKRXqFisdTRkc8hUKhUCgqgBIitZALz7fSeOMdhchQS4fXKuKLAQl1yYTHhuDvsusPeJWge8Be5K2sqmmIEhe4DZCWGPF4lEVEoVAoFCePEiK1kAaZ6aDHkOAoQurhgammEFYAq/C6acrLlIHSIFVRWg5euE2/m0ZG2bxjHoTTjTQFHkNZRBQKhUJx8ighUkuJik/Drntw6FbQqL+Eu/CVcPdaRnxiJFCQBLajMQKa4gUYTvSSgEATm47UvT8qpsQ0wFlccgpOpVAoFLWbKVOmIIQIerVq1eqE7l20aBFCCAYPHuwfc7vd3HfffbRv3564uDgyMzMZMWIE+/btC7p33bp1XHDBBSQnJ1O3bl1uvvlmCgoKqvJopwwlRGopl192PgBJtkJLcHhFhxSlfw4qASIDrCOBJdtDDBu+mBO9xAiaJx2lgSaaBppQPzoKhUIRibZt25KTk+N/rVq16rj37Ny5k3HjxnHuuecGjRcVFbFu3TomTZrEunXreOedd9i8eTOXXXaZf86+ffvo378/zZs3Z+3atSxdupSffvqJ66+/vqqPdkpQdbprKb26d+aNV14kUS/kgFkXiNzQzopblf44EWF6DSLeOBJhhCwsBCCxFRtBlhPpsCGLXQhNotklSYnJp+hkCoVCUTPJy8sLeh8VFUVUVFTYPJvNRnp6+gmvaxgGw4cPZ+rUqaxcuZJjx475ryUlJbFs2bKg+XPnzqVbt27s2rWLrKwsPvroI+x2O08//TSaZv0j8dlnn6VDhw5s27aN5s2bV+CUfzxKiNRS7HY7CfUy4Pe96BgY6KUZM8fDY4JDQzhNROgNAqQpEa7wwBLpsGETlisoMSHx5A+hUCgUVYhtRw427RQ05DSt33uNGjUKGp48eTJTpkwJm75161YyMzOJjo6mR48ezJgxg6ysrDKXf+ihh0hLS2P06NGsXLnyuNvJzc1FCEFycjIATqcTh8PhFyEAMTFW9etVq1bVeCGi7Ou1mDM7dEAISNCKSgublYeU6C4Te6FEc0kchRLNKcMqsepOM1zPCIGMsqNFWbEjyXWSq+gUCoVCUTvYvXs3ubm5/tfEiRPD5nTv3p0XX3yRpUuXMm/ePHbs2MG5555Lfn5+xDVXrVrFggULmD9//gntoaSkhPvuu4+hQ4eSmGj9g7Bfv37s37+fxx57DJfLxdGjR5kwYQIAOTk5lTztH4cSIrWYAX3PASBRKyhXhAgp0dwSe5FEd4NuQFSuiWaCZkh0J0GCRHcakRNtdA092vpjvdQ6VXsYhUKhqOEkJiYGvSK5ZS6++GKuuuoqOnTowIABA/jvf//LsWPHeOutt8Lm5ufnc9111zF//nzq1at33Oe73W6GDBmClJJ58+b5x9u2bctLL73E448/TmxsLOnp6TRp0oT69esHWUlqKso1U4upn5aKFhVPYkmhVcBMRPDLSImtxPq/8L4PLeUusGJHhBOkJtGLjDI9PFqsdSU9I63qDqJQKBR/UpKTk2nRogXbtm0Lu7Z9+3Z27tzJwIED/WOmaVmdbTYbmzdvplmzZkCpCPntt9/4/PPP/dYQH8OGDWPYsGEcOHCAuLg4hBDMnj2bpk2bnsLTVQ01XyopyiWrSXPswiBWOMMvSsDwZsz4xiL0ofEhsDr02gvcVoBrBPRoazyjUf2T2bZCoVCcFhQUFLB9+3YyMjLCrrVq1Yoff/yR9evX+1+XXXYZffv2Zf369f6YFJ8I2bp1K59++il169Yt83n169cnPj6eN998k+joaC644IJTdraqQllEajndu3Zm56b1JFJAsScKYVjWDc1XOwRAC86aCSPUmiJBSNMa07Sg67rdWiAt7fhmRIVCoTjdGDduHAMHDiQ7O5t9+/YxefJkdF1n6NChAIwYMYIGDRowY8YMoqOjadeuXdD9vgBU37jb7eavf/0r69at46OPPsIwDPbv3w9ASkoKDocVnDt37lx69uxJfHw8y5YtY/z48Tz66KP+9WoySojUcnr37MqiV18kWRZypDiCSpbgL1kmZYQS794xTYCU2Aut+iFCetN+DcMSIVKCJrA7DKQJsbGxp/poCoVCUevYs2cPQ4cO5fDhw6SmptKrVy/WrFlDamoqALt27apQ3MbevXv54IMPAOjYsWPQtS+++II+ffoA8PXXXzN58mQKCgpo1aoVzz33HNddd12VnOlUo4RILUfXdZLrNYDfd6MLA0PqpRelRHi8Rc50AWZ4dq9ebOLIc+NMsmMvsXrLCK/wEAHrAAhpEB3vJMaRUisCoBQKheKPZtGiReVeX758ebnXX3zxxaD3jRs3LtNVHsjLL7983Dk1FfVp8ieg45lnIgTE2wpLByV+4eELTo3kltFLPAgJUfke/zxZRqRqTJILoUGXrp2rcvsKhUKhOI1RQuRPwIX9egKQYC8MGvcLD4mVmhsoqqVEeEx0p9faYUi/5aMsYpOt/jIXXdqnCnatUCgUCoUSIn8K6tWrix6VSIKtCJ/aCGxmJ6Q1HOhq0dwSW2FI4TLDuldGSgMWktgkJ3Y9gbS01FNyDoVCoVCcfigh8iehabMW2DSDGL0EzIBaIVKCIUt7ykiJXmKiuySaJ9hXo/msIhF0SHS8C80madmi9Sk9h0KhUChOL5QQ+ZPQr3cPABJshThy3TiOuYg64iTqsIuoQ06EKcFjYisy0byiRBjBrhgJ3uqq4UokNtmqUzLg4t6n8BQKhUKhON1QQuRPwpntW2OYGolaAbYSE80jrWqpgO4xid1fjP2IK8hSEipEBF6rSIgO0e0GsXVKkIadM1o0+yOOo1AoFIrTBCVE/iTouk5avYbERLnQdaP0gtdNI0ywGRK92GMFqR4riVzGXRIUtKrbDdJbHsHmMOnY/mwrtVehUCgUiipCCZE/EV06dwQgNr7YPyaM0jgQS5BI7MdKsBV5ECXuMoqcWf/THQbprY5gjzY4t8cF3HbXtaf4BAqFQqE43ahWITJlyhSEEEGvVq1a+a+XlJRw++23U7duXeLj47nyyis5cOBANe64ZnP+BVY33th4K802KGgVrEDVQhfCbSI0gVbgRD9ahHD6BEmAJcThIaPlEWwOg77nXcy11//1DzyJQqFQKE4Xqt0i0rZtW3JycvyvVatW+a/dc889fPjhh7z99tusWLGCffv2ccUVV1Tjbms2ycnJOOxJxMSVgDTRjOCsGOE2g9J6wbKQ6H5B4rHKvEd7yGh1BN1hcGH/gVwzfPAfeg6FQqFQnD5Ue4l3m81Genp62Hhubi4LFizg9ddfp1+/fgAsXLiQ1q1bs2bNGs4+++yI6zmdTpzO0k60eXl5p2bjNZTmzVvw8y/fEB3jwpVnD7omzABhElK8zCdIoqMLqNe8EKFLLjj/Mv465NI/YtsKhUKhOE2pdovI1q1byczMpGnTpgwfPpxdu3YB8N133+F2u+nfv79/bqtWrcjKyuKrr74qc70ZM2aQlJTkf/naKJ8u9PVWWY2NLbIGpATDROQXQ2DdECOkmBmSxKxiUtsWIIRgxPAbuOpqJUIUCoVCcWqpViHSvXt3XnzxRZYuXcq8efPYsWMH5557Lvn5+ezfvx+HwxHWwrh+/fr+FsiRmDhxIrm5uf7X7t27T/Epahbt2rXC9GjEJjqtFN38YrSDxxCFJcEZLwGpu0KX1GtTSHLTEuwilkmTJtGrd/dq2L1CoVAoTjeq1TVz8cUX+//coUMHunfvTnZ2Nm+99RYxMTGVWjMqKoqoqKiq2mKtQ9M0hCuOqNh87MeOYBZb4kOaJhgG6Lo3MNXEFmMSU9dNfIYTe5xJalIj7p96b6W/9gqFQqFQVJRqd80EkpycTIsWLdi2bRvp6em4XC6OHTsWNOfAgQMRY0oUpVx6qVX9NCbZ5Y8FEUIgPE5ikp2kZOfRoHsemd3zqNO8GFuMSZtmXXh45v1KhCgUCoXiD6Xag1UDKSgoYPv27Vx33XV06dIFu93OZ599xpVXXgnA5s2b2bVrFz169KjmndZs+l7Sm/98/iHRdUpwHXUT2yiKqBSTqGQT4ZWemnTQsH5jevTqRs/zziI6Orp6N61QKBSK05JqtYiMGzeOFStWsHPnTr788ksuv/xydF1n6NChJCUlMXr0aMaOHcsXX3zBd999x6hRo+jRo0eZGTMKi8TERKJlLHHZDjL66iQ19+BIMIm31aFJWnv+PuF+nnnuKe5/6F76XXiuEiEKhUJxinj00UcRQnD33XeXOeedd96ha9euJCcnExcXR8eOHXnllVeC5oTW3PK9HnvsMf+cxo0bh11/9NFHT9XRqoxqtYjs2bOHoUOHcvjwYVJTU+nVqxdr1qwhNdVqM//EE0+gaRpXXnklTqeTAQMG8Mwzz1TnlmsNl/7lUpYt/ZQmjZvQ64KetDuzLbquV/e2FAqF4rThm2++4bnnnqNDhw7lzktJSeH++++nVatWOBwOPvroI0aNGkVaWhoDBgwAICcnJ+ieJUuWMHr0aL/HwMdDDz3ETTfd5H+fkJBQRac5dVSrEFm0aFG516Ojo3n66ad5+umn/6Ad/XkYMOgCBgy6oLq3oVAoFKclBQUFDB8+nPnz5/PII4+UO7dPnz5B7++66y5eeuklVq1a5RciobGR77//Pn379qVp06ZB4wkJCbUujrJGxYicCqQ3WPN0K2ymUChOLb7fKVKGNmtSVBce6QLz+PMqtS7hnyPlZWnefvvtXHrppfTv3/+4QiQQKSWff/45mzdvZubMmRHnHDhwgP/85z+89NJLYdceffRRHn74YbKyshg2bBj33HMPNlvN/qiv2burAvLz8wFOu8JmCoXijyE/P5+kpKTq3sZpjcPhID09neX7Xz1lz4iPjw/7HJk8eTJTpkwJm7to0SLWrVvHN998c8Lr5+bm0qBBA5xOJ7qu88wzz3DBBZGt2i+99BIJCQlhLU/GjBlD586dSUlJ4csvv2TixInk5OQwe/bsE95HdfCnFyKZmZns3r2bhISEGt/CPi8vj0aNGrF7924SExOreztVgjpT7UCdqeJIKcnPzyczM7PK11ZUjOjoaHbs2IHL5Tplz5BShn2GRLKG7N69m7vuuotly5ZVKBEgISGB9evXU1BQwGeffcbYsWNp2rRpmNsG4F//+hfDhw8PW3/s2LH+P3fo0AGHw8Ett9zCjBkzanR9LSGVXbHGkJeXR1JSErm5uX+qDwN1ppqPOpNCUTW89957/gxQH4ZhIIRA0zS/xeN43HjjjezevZuPP/44aHzlypWcd955rF+/njPPPLPcNX766SfatWvHpk2baNmyZeUO9Afwp7eIKBQKhULxR3H++efz448/Bo2NGjWKVq1acd99951w9qJpmkENXH0sWLCALl26HFeEAKxfvx5N00hLSzuxzVcTSogoFAqFQlFFJCQk0K5du6CxuLg46tat6x8fMWIEDRo0YMaMGYDVrLVr1640a9YMp9PJf//7X1555RXmzZsXtE5eXh5vv/02jz/+eNhzv/rqK9auXUvfvn1JSEjgq6++4p577uHaa6+lTp06p+i0VYMSIjWIqKgoJk+eXKN9eRVFnal2oM6kUPxx7Nq1C00rrSdaWFjIbbfdxp49e4iJiaFVq1a8+uqrXH311UH3LVq0CCklQ4cODVszKiqKRYsWMWXKFJxOJ02aNOGee+4JihupqagYEYVCoVAoFNVGjWp6p1AoFAqF4vRCCRGFQqFQKBTVhhIiCoVCoVAoqg0lRBQKhUKhUFQbSojUIJ5++mkaN25MdHQ03bt35+uvv67uLVWaGTNmcNZZZ5GQkEBaWhqDBw9m8+bN1b2tKuNEWnvXBvbu3cu1115L3bp1iYmJoX379nz77bfVva1KYxgGkyZNokmTJsTExNCsWTMefvhh1Q9GoajBKCFSQ3jzzTcZO3YskydPZt26dZx55pkMGDCAgwcPVvfWKsWKFSu4/fbbWbNmDcuWLcPtdnPhhRdSWFhY3Vs7aU60tXdN5+jRo5xzzjnY7XaWLFnCzz//zOOPP17jaw6Ux8yZM5k3bx5z587ll19+YebMmcyaNYt//vOf1b01hUJRBip9t4bQvXt3zjrrLObOnQtYVfUaNWrEnXfeyYQJE6p5dyfP77//TlpaGitWrOC8886r7u1UmoKCAjp37swzzzzDI488QseOHZkzZ051b6tSTJgwgdWrV7Ny5crq3kqV8Ze//IX69euzYMEC/9iVV15JTEwMr7566hqiKRSKyqMsIjUAl8vFd999R//+/f1jmqbRv39/vvrqq2rcWdWRm5sLQEpKSjXv5OQIbO1d2/nggw/o2rUrV111FWlpaXTq1In58+dX97ZOip49e/LZZ5+xZcsWADZs2MCqVau4+OKLq3lnCoWiLFRl1RrAoUOHMAyD+vXrB43Xr1+fTZs2VdOuqg7TNLn77rs555xzwkof1yYq09q7JvPrr78yb948xo4dy9///ne++eYbxowZg8PhYOTIkdW9vUoxYcIE8vLyaNWqFbquYxgG06ZNY/jw4dW9NYVCUQZKiChOObfffjsbN25k1apV1b2VSlPZ1t41GdM06dq1K9OnTwegU6dObNy4kWeffbbWCpG33nqL1157jddff522bduyfv167r77bjIzM2vtmRSKPztKiNQA6tWrh67rHDhwIGj8wIEDpKenV9OuqoY77riDjz76iP/97380bNiwurdTab777jsOHjxI586d/WOGYfC///2PuXPnnnBr75pERkYGbdq0CRpr3bo1ixcvrqYdnTzjx49nwoQJXHPNNQC0b9+e3377jRkzZigholDUUFSMSA3A4XDQpUsXPvvsM/+YaZp89tln9OjRoxp3VnmklNxxxx28++67fP755zRp0qS6t3RS+Fp7r1+/3v/q2rUrw4cPZ/369bVOhACcc845YSnVW7ZsITs7u5p2dPIUFRUFNRMD0HUd0zSraUcKheJ4KItIDWHs2LGMHDmSrl270q1bN+bMmUNhYSGjRo2q7q1Vittvv53XX3+d999/n4SEBPbv3w9AUlISMTEx1by7inMirb1rG/fccw89e/Zk+vTpDBkyhK+//prnn3+e559/vrq3VmkGDhzItGnTyMrKom3btnz//ffMnj2bG264obq3plAoykIqagz//Oc/ZVZWlnQ4HLJbt25yzZo11b2lSgNEfC1cuLC6t1Zl9O7dW951113VvY2T4sMPP5Tt2rWTUVFRslWrVvL555+v7i2dFHl5efKuu+6SWVlZMjo6WjZt2lTef//90ul0VvfWFApFGag6IgqFQqFQKKoNFSOiUCgUCoWi2lBCRKFQKBQKRbWhhIhCoVAoFIpqQwkRhUKhUCgU1YYSIgqFQqFQKKoNJUQUCoVCoVBUG0qIKBQKhUKhqDaUEFEoFAqFQlFtKCGiqFZ27tyJEIL169efkvWFELz33nuVvn/58uUIIRBCMHjw4HLn9unTh7vvvrvSz6pJ9OnTx3/uU/W9USgUClBC5LTm+uuvP+6H66mmUaNG5OTk+Pu1+D74jx07Vq37CmXz5s28+OKL1b2NP4x33nmHr7/+urq3oVAoTgNU0ztFtaLrOunp6dW9jeOSlpZGcnJydW8Dl8uFw+E45c9JSUkhLy/vlD9HoVAolEVEUSYrVqygW7duREVFkZGRwYQJE/B4PP7rffr0YcyYMfztb38jJSWF9PR0pkyZErTGpk2b6NWrF9HR0bRp04ZPP/00yF0S6JrZuXMnffv2BaBOnToIIbj++usBaNy4MXPmzAlau2PHjkHP27p1K+edd57/WcuWLQs70+7duxkyZAjJycmkpKQwaNAgdu7cWeGvTWFhISNGjCA+Pp6MjAwef/zxsDlOp5Nx48bRoEED4uLi6N69O8uXLw+aM3/+fBo1akRsbCyXX345s2fPDhI8U6ZMoWPHjrzwwgs0adKE6OhoAI4dO8aNN95IamoqiYmJ9OvXjw0bNgSt/f7779O5c2eio6Np2rQpU6dO9X//pJRMmTKFrKwsoqKiyMzMZMyYMRX+OigUCsXJooSIIiJ79+7lkksu4ayzzmLDhg3MmzePBQsW8MgjjwTNe+mll4iLi2Pt2rXMmjWLhx56yC8ADMNg8ODBxMbGsnbtWp5//nnuv//+Mp/ZqFEjFi9eDFiukJycHJ588skT2q9pmlxxxRU4HA7Wrl3Ls88+y3333Rc0x+12M2DAABISEli5ciWrV68mPj6eiy66CJfLVZEvD+PHj2fFihW8//77fPLJJyxfvpx169YFzbnjjjv46quvWLRoET/88ANXXXUVF110EVu3bgVg9erV3Hrrrdx1112sX7+eCy64gGnTpoU9a9u2bSxevJh33nnHH69x1VVXcfDgQZYsWcJ3331H586dOf/88zly5AgAK1euZMSIEdx11138/PPPPPfcc7z44ov+9RcvXswTTzzBc889x9atW3nvvfdo3759hb4GCoVCUSVUc/dfRTUycuRIOWjQoIjX/v73v8uWLVtK0zT9Y08//bSMj4+XhmFIKaXs3bu37NWrV9B9Z511lrzvvvuklFIuWbJE2mw2mZOT47++bNkyCch3331XSinljh07JCC///57KaWUX3zxhQTk0aNHg9bNzs6WTzzxRNDYmWeeKSdPniyllPLjjz+WNptN7t271399yZIlQc965ZVXws7kdDplTEyM/PjjjyN+HSLtJz8/XzocDvnWW2/5xw4fPixjYmLkXXfdJaWU8rfffpO6rgftR0opzz//fDlx4kQppZRXX321vPTSS4OuDx8+XCYlJfnfT548Wdrtdnnw4EH/2MqVK2ViYqIsKSkJurdZs2byueee8z9n+vTpQddfeeUVmZGRIaWU8vHHH5ctWrSQLpcr4rmlDP/eKBQKxalAxYgoIvLLL7/Qo0cPhBD+sXPOOYeCggL27NlDVlYWAB06dAi6LyMjg4MHDwKWVaNRo0ZBMSDdunU7Zftt1KgRmZmZ/rEePXoEzdmwYQPbtm0jISEhaLykpITt27ef8LO2b9+Oy+Wie/fu/rGUlBRatmzpf//jjz9iGAYtWrQIutfpdFK3bl3A+vpcfvnlQde7devGRx99FDSWnZ1Nampq0DkKCgr86/goLi72n2PDhg2sXr06yMJiGAYlJSUUFRVx1VVXMWfOHJo2bcpFF13EJZdcwsCBA7HZ1K8EhULxx6J+6yhOCrvdHvReCIFpmlX+HE3TkFIGjbnd7gqtUVBQQJcuXXjttdfCrgV+0FcFBQUF6LrOd999h67rQdfi4+MrtFZcXFzY2hkZGWHxJoA/vqSgoICpU6dyxRVXhM2Jjo6mUaNGbN68mU8//ZRly5Zx22238dhjj7FixYqw76lCoVCcSpQQUUSkdevWLF68GCml3yqyevVqEhISaNiw4Qmt0bJlS3bv3s2BAweoX78+AN9880259/gyQgzDCBpPTU0lJyfH/z4vL48dO3YE7Xf37t3k5OSQkZEBwJo1a4LW6Ny5M2+++SZpaWkkJiae0Bki0axZM+x2O2vXrvVbho4ePcqWLVvo3bs3AJ06dcIwDA4ePMi5554bcZ2WLVuGfT2O9/XxnWP//v3YbDYaN25c5pzNmzfTvHnzMteJiYlh4MCBDBw4kNtvv51WrVrx448/0rlz5+PuQaFQKKoKFax6mpObm8v69euDXrt37+a2225j9+7d3HnnnWzatIn333+fyZMnM3bsWDTtxH5sLrjgApo1a8bIkSP54YcfWL16NQ888ABAkMsnkOzsbIQQfPTRR/z+++8UFBQA0K9fP1555RVWrlzJjz/+yMiRI4MsDf3796dFixaMHDmSDRs2sHLlyrDA2OHDh1OvXj0GDRrEypUr2bFjB8uXL2fMmDHs2bPnhL9m8fHxjB49mvHjx/P555+zceNGrr/++qCvS4sWLRg+fDgjRozgnXfeYceOHXz99dfMmDGD//znPwDceeed/Pe//2X27Nls3bqV5557jiVLlpT5tQk8a48ePRg8eDCffPIJO3fu5Msvv+T+++/n22+/BeDBBx/k5ZdfZurUqfz000/88ssvLFq0yP/1f/HFF1mwYAEbN27k119/5dVXXyUmJobs7OwT/jooFApFlVDdQSqK6mPkyJESCHuNHj1aSinl8uXL5VlnnSUdDodMT0+X9913n3S73f77e/fu7Q/O9DFo0CA5cuRI//tffvlFnnPOOdLhcMhWrVrJDz/8UAJy6dKlUsrIAZEPPfSQTE9Pl0II/1q5ubny6quvlomJibJRo0byxRdfDApWlVLKzZs3y169ekmHwyFbtGghly5dGhSsKqWUOTk5csSIEbJevXoyKipKNm3aVN50000yNzc34teorODZ/Px8ee2118rY2FhZv359OWvWrLCvh8vlkg8++KBs3LixtNvtMiMjQ15++eXyhx9+8M95/vnnZYMGDWRMTIwcPHiwfOSRR2R6err/+uTJk+WZZ54Ztq+8vDx55513yszMTGm322WjRo3k8OHD5a5du/xzli5dKnv27CljYmJkYmKi7Natm3z++eellFK+++67snv37jIxMVHGxcXJs88+W3766adBz1DBqgqF4o9ASBnieFcoTiGrV6+mV69ebNu2jWbNmlX3do7L8uXL6du3L0ePHv1DCprddNNNbNq0iZUrV57yZx2PnTt30qRJE77//ns6duxY3dtRKBR/UlSMiOKU8u677xIfH88ZZ5zBtm3buOuuuzjnnHNqhQgJpGHDhgwcOJA33nijStf9xz/+wQUXXEBcXBxLlizhpZde4plnnqnSZ1SGiy++mP/973/VvQ2FQnEaoCwiilPKyy+/zCOPPMKuXbuoV68e/fv35/HHHw9LPa2pFBcXs3fvXsCKDanqcvRDhgxh+fLl5Ofn07RpU+68805uvfXWKn1GZdi7dy/FxcUAZGVl/SFl5RUKxemJEiIKhUKhUCiqDZU1o1AoFAqFotpQQkShUCgUCkW1oYSIQqFQKBSKakMJEYVCoVAoFNWGEiIKhUKhUCiqDSVEFAqFQqFQVBtKiCgUCoVCoag2lBBRKBQKhUJRbfw/h7Npg6e6WpsAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ds = mikeio.read(\"../tests/testdata/wind_north_sea.dfsu\")\n", "ws = ds[\"Wind speed\"][0]\n", "ws.plot();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Step 2. extract geometry" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "shapely.geometry.multipolygon.MultiPolygon" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "shp = ds.geometry.to_shapely()\n", "type(shp)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Geopandas does not like multipolygon, it should be a list of polygons" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "poly_list = [p for p in shp.geoms]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Step 3. Create a geopandas dataframe" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
wind_speed
09.530760
19.652719
29.806072
38.775489
411.013206
\n", "
" ], "text/plain": [ " wind_speed\n", "0 9.530760\n", "1 9.652719\n", "2 9.806072\n", "3 8.775489\n", "4 11.013206" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "df = pd.DataFrame({'wind_speed':ws.to_numpy()})\n", "df.head()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "gdf = gpd.GeoDataFrame(df,geometry=poly_list, crs=4326)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Step 4. Save to shapefile" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "gdf.to_file(\"wind_speed.shp\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Step 5...\n", "Do further work in QGIS\n", "\n", "![QGIS](../images/dfsu_qgis.png)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Would you prefer to have this workflow to be a method on the `mikeio.Dfsu` class?\n", "\n", "Post an issue on [GitHub](https://github.com/DHI/mikeio/issues) !" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Clean up" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "import os\n", "\n", "files = [\"wind_speed\"]\n", "\n", "exts = [\"cpg\",\"dbf\",\"shp\",\"shx\"]\n", "\n", "for file in files:\n", " for ext in exts:\n", " filename = f\"{file}.{ext}\"\n", " if os.path.exists(filename):\n", " os.remove(filename)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": ".venv", "language": "python", "name": "python3" }, "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.13.0" } }, "nbformat": 4, "nbformat_minor": 4 }