diff --git a/notebooks/Sentinel-1 wind.ipynb b/notebooks/Sentinel-1 wind.ipynb index 369f445..51c07a5 100644 --- a/notebooks/Sentinel-1 wind.ipynb +++ b/notebooks/Sentinel-1 wind.ipynb @@ -32,7 +32,7 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 25, "metadata": {}, "outputs": [], "source": [ @@ -41,12 +41,50 @@ "#fn = '/vagrant/lib/project/products/WIND_WIND_UNPROJECTED_S1A_IW_GRDM_1SDV_20190704T052723_20190704T052755_027963_03284A_AF7D.nc'\n", "fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc'\n", "fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165553_20190709T165625_028043_032AC4_4501.nc'\n", - "#fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190704T052723_20190704T052755_027963_03284A_AF7D.nc'" + "#fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190704T052723_20190704T052755_027963_03284A_AF7D.nc'\n", + "fn = '/home/vagrant/project/products/TMP_WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc'" + ] + }, + { + "cell_type": "code", + "execution_count": 26, + "metadata": {}, + "outputs": [], + "source": [ + "n = Nansat(fn)" + ] + }, + { + "cell_type": "code", + "execution_count": 29, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAADxCAYAAAAKooq5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOy9f4xk2XXf97lT1VXdNd2zvTu9O01ur9mMZqMRubFW9ICiISFgJDGiCCdUACGggiiMYICWQwGW4yCh/IcRyxKovwwksGSCgA3RgCNCgS2JcGgpNGFGIuKVtGSW9lAcYkdiU2zu9ix7hj3Tvd1d3a/m5o97v++c97q6u2a2Zram5h2gUFXvx333vbp1vvec7znnhhgjjTTSSCONNHIvcubN7kAjjTTSSCMPrzQg0kgjjTTSyD1LAyKNNNJII43cszQg0kgjjTTSyD1LAyKNNNJII43cszQg0kgjjTTSyD1LAyKNTKWEEP5pCOG1EMIVt+2JEMLnQggv5/fHjzn3/SGEr4cQroUQPvbget1IIyfLJI7rBkQamVb5DeD9tW0fAz4fY3wW+Hz+XpEQQgv4NeAngHcAPx1CeMf97WojjYwsv8GEjesGRBqZSokx/gFws7b5g8Cn8udPAT855NR3A9dijH8eYzwAPp3Pa6SRN10mcVy3x9FII42MQy6GEHdHPPZV+Cqw7zZ9Msb4yVNOuxBjfBUgxvhqCOGpIcc8DXzLfV8HfnDEbjXSyBGZ9nHdgEgjEyO7wN8Y8dj/FfZjjJfvQzfCkG1NbaBG7lmmfVw3INLIxEjgvg/I6yGEt+TZ2luA14Ycsw48476vAK/c3241Ms0y7eO64UQamRg5A8yN+LpH+Qzw4fz5w8DvDjnmT4BnQwhvDyF0gA/l8xpp5J5k2sd1AyKNTIwEYGbE16lthfCbwL8DvjeEsB5C+OvArwLvCyG8DLwvfyeE8NYQwmcBYowF8PPA7wNfA34rxvjVsd1kI4+cTPu4Dk0p+EYmRd4WQhw1eP1/gC/dJ99xI42MVaZ9XDecSCMTI5qxNdLINMm0j+sGRBqZGHkABGQjjTxwmfZxPc331shDJtM+Y2vk0ZRpH9cNiDQyMaIolkYamSaZ9nHdgEgjEyPTPmNr5NGUaR/XDYg0MlHSDMhGplGmeVxP87018pDJtM/YGnk0ZdrHdQMijUyMTHsUSyOPpkz7uJ7me2vkIZNpJyAbeTRl2sd1AyKNTIxMu9nfyKMp0z6uGxBpZGJk2s3+Rh5NmfZxPc331shDJtM+Y2vk0ZRpH9cNiDQyMTLtM7ZGHk2Z9nE9zffWyEMm0z5ja+TRlGkf1w2INDIxEpjuKJZGHk2Z9nHdgEgjEyMBmBl1RBb3syeNNDI+mfZx3YBIIxMjIUB7iv9sjTyaMu3jugGRRiZGQoCZ1pvdi0YaGa9M+7hu1lhvZGJEM7ZRXie3E743hPCSe90OIfxC7Zj3hhBuuWP+3v28t0YeXRnXuE5tTd7YbiyRRiZGQoCZ7htvJ8b4deD51GZoAd8GfnvIoX8YY/xrb/yKjTRyvIxrXMNkju0GRBqZHLk/AfU/CvxZjPGbY2+5kUZGkfuXKDIRY7txZzUyOaI/2yiv0eVDwG8es++vhhC+EkL41yGEd95jrxtp5GS5P+MaJmRsN5ZII5Mlo4/IpRDCi+77J2OMn/QHhBA6wH8J/OKQ878MvC3GuBNC+ADwO8Czd9/hRhoZQcY4rmGyxnYDIo1MjgRg9CiWzRjj5VOO+QngyzHG6/UdMcbb7vNnQwi/HkJYijFujtyDRhoZRcY/rmGCxnYDIo1Mjozfd/zTHGPuhxCWgesxxhhCeDfJtXtjrFdvpBG4X5zIxIztBkQamRwJwJiiWEIIPeB9wN9w234OIMb4CeCngL8ZQiiAPeBDMcY4nqs30oiTMY5rmLyxHZr/TSOTIpfPhvjiO0Y7NrzIl0Y0+xtp5E2VaR/XjSXSyOTItNfMbuTRlCkf11N8a408lDLF5SEaeYRlisd1AyKNTI5M+YytkUdUpnxcT/GtNfLQyZT/2Rp5RGXKx/UU31ojD52MOYqlkUYmQqZ8XDcg0sjkyJTP2Bp5RGXKx/UU31ojD51M+Z+tkUdUpnxcT/GtNfLQyd2Vh2ikkYdDpnxcNyDSyOTIlM/YGnlEZcrH9RTfWiMPnUw5AdnIIypTPq4bEGlkcmTKZ2yNPKIy5eN6im+tkYdOpvzP1sgjKlM+rqf41hp5KKUZkY1Mo0zxuJ7iW2vkoZMpj2Jp5BGVKR/XDYg0Mjky5WZ/I4+oTPm4nuJba+ShkymPYmnkEZUpH9cNiDQyOTLlM7ZGHlGZ8nE9xbfWyEMnU/5na+QRlSkf11N8a408dDLlf7ZGHlGZ8nE9xbfWyEMpY4piCSGsAdvAACjq61aHEALwvwEfAHaB/z7G+OXxXL2RRmoyxuisSRvbDYg0Mjky/hnbfxZj3Dxm308Az+bXDwL/OL830sh45f5YIhMzts/cr4ZDCO8PIXw9hHAthPCx+3WdRqZIFMUyyuuNyweBfxaTvAAshhDecmoXm3HdyN3Kgx3XcI9j+17lvoBICKEF/BoJEd8B/HQI4R3341qNTJFoxjbKC5ZCCC+610dqrUXg/w4hfGnIPoCngW+57+t52/Hda8Z1I/ci4x3XcB/G9huR++XOejdwLcb45wAhhE+T0PFP79P1GpkGuTuzf7PuC67JD8UYXwkhPAV8LoRwNcb4B7Wr1SWecs1mXDdy9zLecQ33Z2zfs9wvEBmGhBWfXEbQjwCcPXv2r1y6dOnEBr/NK3fdiTDG5/YYtxjQ4gx3iAQiodJ++n6HyJly+xnucIcA+dgI+Vssz7nDGXfenbxXxwdaDLhD4AwRyhaoXVf70l3f4rGx3fc45dUvbWzGGJ889oAxloeIMb6S318LIfw2CQD8H20deMZ9X4FTB9mp4xoezNgep2x86UmYnYEDkkulA9wh+Ski6Tc5AOZIVG4E+iTtMZP3SQpsmPbycQHYz+2E3HaHRPnOAof5nHY+X+203fVI5y6/6819VsfJxpdePX5sj7nsyX0a2/cs9wtETkXCGOMngU8CXL58Ob744ovHNvZ3+Xt3dfFWORLHJ/8pfwhAJ/9jBrQY5JHRYgBAj1126ZX759hlQJsBLTr0GeTH3WOXbRZoMWCPOebZZoeFsu0DOmWbi2wxoEWLQXm9ghbtvF/X98f8X3xg7Pc/Dvnl8PFvnnjAmAjIEMJZ4EyMcTt//s+BX6od9hng57M18YPArRjjq6c1PWTbkZnK/Rzb90M+/n2/lJ77FjBPUjlLJEW+kw+6ClwCFkmK/2r+PA+s5fcCeAn4sbzteeD3gPfk7fvAan6fBV4AnsjfF/N7AWzka+pvPJv3nYWfffHNf17D5OPhHxw/tsdIrN/HsX3Pcr9AZGxIeNyfbFCDdindFkVFwWpfelXBpc2gopBb7jyv9FsUDGgxz3a57YAOi2yxxSIdDmhRlMCg14A228zTY48BbVoUHNCl70CioMUePQpaFMzR5aAEE39MiwF9unTpc0CXNMWr3nufDt3KtPAhk0BSGG9cLgC/nSIdaQP/R4zx90IIPwcQY/wE8FlSCOQ10pz4Z0dod6wzvP+C/5MD/kr5u/qJwhy7AAxoc0CHXXp06LPFInv0WGSLFgO2WKRPpxyXGn3b/QUWu1vs0svjt81g0KLT6vPa9QssLm0xKFpsby0kBf4iScHPkgDkiyQwuJY7e5kEKOv5mDYJVPR9B1gGLuZt/gnt5+0b7vitvF+f1YfngU0SKAFsZbNmGQOXh03GN67h/o3te5b7BSJ/AjwbQng78G3gQ8B/c7eNjDJL0wwcKP+IPXYp8p+xnf9WC2xni8AUtGb8Or7LQflHXmC7/HP22GWOXQ7o0slKvMNB/oOm4xM0HNBjN3/q5j4dOCuh7c7vl8e0GNDN7fVdnzsZKNoM2GYBoLwvgZKeQQKyFr2sfB5KGZPZnzmL7x+y/RPucwQ+epdNj2VcAzzFX/CdP/sp5pc32dlc5LHlGzzVvV7u32aBGzeWONzvcKY94M5Oj5nFbQ5/5xw8R1KomyTl+wIJAK6QlO814BLc2lmGL2AAsE9S6l+Em6tn079/jWQlXCQp89Xc7qV8ziJmdVzJ58sC2aDqglrL+wVIO66f78Esih137hbpfnZym/skENvMx88G2I+wGYAb9/Ko33wZr5v2fo3te5b7Ep0VYyyAnwd+H/ga8Fsxxq/eTRv/C3+//NzK83RvSQg4OhzQoV9R4hJvYdQtFyn2Dv0SQHStTmWWn6wQKfw9eiywXVobun4ng8U2C5XrdjIw7TEHwAHdsi+7zNHOYCY3mMS7wur7vMXUztdIz6JfgsxDKXcXxfLAZRzjGmCxv8FBvwNFYGdzEYo257s3eK1/gQHtcvy89fwrzC9uc6dowU7gcONcUtKQFOxq/nyRNFNfJini1bx/naSgZRnIkpDiL0gKeyVvuwz8m3ycjl0nKfz5vF9treWXQKCg6pKSRTJPAhDvmrqK/Y4rGNDomnJ3FcD+LhDy+efv6jlPjEz4uH6jct+6HWP8LMmsumv5n/jlUhXWXVA2+26XLqp2qazTvqRUO+4z2R3UoeWUbLJGulkpp+O6HLDNQgkqdfdQK1sFyeLoZjAqOCiv12c78xsJfDrOvSUrIrmwkpWzl11i/fJ6HrQKd44HDzCLRdJjj20W+K/4Hbr06WfnmPoii8XzObonteOvc+Ceoed/BNg3OE+HgwqH43+TOXZLK7A9GPDLQ3/tmkz4H+mNjGuAuVs3WXxslx573Jpf5LGlLVrtAa+9/hTnz95gmwUW2C6t4t2dHuzM8tilDW6tX4CNkBT/LMkyuExS5HIbFfm7lLvS0RZJSv89mBWw4fbJrfRekiVzkQQAOyRQKvL2H87vi1StDVkRcmkVJLD4ybxf/dJ+WSHiXPbz+ywGgADlBOplaD/Lx8MvATFZKCsYH7Ocz1kkPR+5w9ruGgK5NmZNCejUP7ncfqx2X9rfrn0u0/3+ASfKhI/rNyL3LdnwXiUByOCIwpTYLN8sEx3ruZBOhVtI1kTiJlrlfq/cvYLvsXtEyVq7Sdnv0svKNHEnfWdhqB1PdvvPRfm9KN1U7Wx39GtK/6AEqkHp9gJz3VkfC3aZo8cuLYqyPwtsl/0ZBhQC13rfJeqrjpFltMtceT9z7NKp9UfPvp3deCOJzP5RXg+hPMVf8Mxj3+JxtthljseWEjFwvrXJ4tktvvntVW5cP89BftY9drlTtHhsdYNbV5ZhJySFdo2kxJZIik58wjKm2HYw99QKSRl6lxQkZbiJKTgpxEskALiW9wm0IF1/Pn/XuXKLeREgtUlKecu1q+tJqW/k/VxPfVqHarzCLnARinVSKFeOb7i2a1FcG1RdYQoKEJjIjSbw2sfAYT5vFxituGOvuf6LC9rBLKVRhvaUj+uJApG/zcdP3C+QkMLVZ82OwZQXmNUixauZdR04Os4VZiR1t1TqHljEZ+iYQdmDpJh3Mrl+QCe7pKrK2/fZiHtT6HJv+eO9ZeBdWwYuRQWEBJZyrx3QKS0aWSYi9we0sqI3ENU9eysIEqBon+dz9uhli66gx27FevO/x6kyxWb/f8xXWGCbPXpcH1xILtH2gLlucr/OsQtFi6cuXC8DNDY5z8zsAbdeWIbFmLiLVZKiFIiskxSfoog3MEWomfYGCTwgKT8pyYu5jZdICl2uJagCz1I+/ocxa+Na3i9rQqCiFxhIiQvZyNeU1bGJKeeLwOIFU+iLCoS7TioT9W1SPPEf2X3M98yiUJ/1jNQHXV+goPbVP9w+WUF6vhtYpJqere5V9z6KTPG4hgkCkb/Nx0+0PjzP4C0PAUirdJ8IOPqV97R/74hrRiJmRcpT9sICO/SzEpb14mfWFouV2vQWiQBF4snxXXql1dDNXIa3APq5Px5QfN8VsWX8Sq8EGCkhRfh4QFQbAg4fAebdUn06FctuWEhxuu5ceU7d9dfNLsaRQ64ffHmIByI/yP9TPj/9Jnv0aLWyq4/kVn32bX9aumZ36SViff1cUrovZKUq4nwZI743STP8DZICbbt9SySX13P58xI2C3+BrLxzG7IypFAFIpp1L+b3WSyKaidvex6b2V8iKfMtzPUlRax+LeZ70Gz+GhZWPE+OyoJkeTxF0uaHwLuA27B1I523jAHqWm5Xocq6JzCrTfcoQMBt36cKgHLLCTBmMdDaccedJlM6riUTASLfYmMoIVwn0uszem2vnysAaDvgqLtafEiwV34i2Xvs0afDXg67Paj9wokbmc/uLCPp5TrTdVLEVac8T5+T+ynxIuqrwgcSeKSZvVfuyYpIoOgJ+nouyjBXH1SjvHSvskjqgCGg1PZ+vn9tP/rMB+7eLJy6tFiK4ROEikzhjO1PeZ1v8PaSs7sxWGKhtc15brCUNVafLotscUCX89woJw6Ha+d47NJGUlxSWuuYcnyJqjKUW8eT3H6WLUWqnA+R7t46kDtMEWDiStawUNyt3MYqBkCybARgl10fxB/INSaXk0BgHwsXVs4Iisaay5+/AZwjWSY3gL3UbuGucRlzV4FZXLJ4ZA0JDNVvD47eBeZdeKvuOPL7cv58WtTYFI5rLxPT7WG5GsOkvl8gIAshuVT2SsUu5Sp3lMGLhcNK+cqdpe8CF0VKKd5e121nV46us8ViSYwr70PieZDUdpVI1zOQcvb3mawOCxPWfgGO3XPa5/uqEGZvWehzN7u1dA1ZJ7LEtN1LN4OmrMC6O7GFPWUP/LvdHqbtjpEpXHfhDmfKsbfANlutx9M7i/TYZY9eDgjvl5bk1uuL7GwswWLk1rVlm90XpOQ9WSEXSYrwixixrhm/XFRrGFEti0XPWBaArA9ZBeIULmFE9HtzW+JHRDgrp0Ngs5TbLLCQ37Y7TtaK2l2ufQdzqc2fz8evk1xZh6Q0iZukLMXbcC1HrAk4ZM0oOED5K94Npfb9WPMAo+g0AausEPXTE+3AqVFjUziuvUyEJQJVS8OLj0A6br+UsxS7J6L9MT6U14cCd/Of2Cwdy8GQVFVmizqfIeLe99EnGOpcH3WlvggMZIFUSe50HW/NeNK9Dkr1/vnnlK7ZLz97HmSYyPpoZ5AY0C7daAbQB0O4pnbZ7kFWkcddoyJTPGMTd5SCHjplbtIiW/TY5bsZVCRPrL4Ca4EzS6+btbFOshBE8EopK3dDbplrmMIUwU7eJvePV+jy+wtMpNj9TP0Kxn/IslgjAY3cPFLg4mpmXXvzGNjIWtH3i5hVUXIimBXFNRKIzOXXTD44T9J2MCtL7cpq0LvuVQCx6L4LPGSJ7Lv2VjGLZNYdJ8tuFJnicQ0TAiL1GlfDwKIaiTRw7qNBaSn4nAnt7+RQ1zl2K4ADSXnP5Vm62rdopnapJJUk6An0Yf1KkUtzpbUDzp2Tt0qkWE+O6DIwazNgN/9pPKge5Gxn76bS+YryqXMWAge1K5E7zdoxy8sDSvV3KDJR3y1Jd7n0TrMqh8qURbFEQjkBmCOF9kIaZ4+z5aLe2tx4/TzXbz3FztYCN194GvbhzhfOmltG7h6fxPc8ScfK16+5zzLJ6ngeszBewJSf3FJSlgIdEeMCDJ2rkF4p0U0symvLHadrCyCkqGXBaFbvXVr7JFfUPrAerc0lYP82yQyCBB6vkayRc/n7rkWAeatDQCXj14f9ClRE+AuwZGXoGPKz0fMSnySXWRnee8ipMmXj2stEgAjEI9aGL19SD2X171LCsiQUFSU3VzdHXWkW3eEA5TCo/WHRTLruQRksnJRoPc+iSqp3XJ+LClFdJ+PljvKAo+2eTFdQgZT7HnMlQKT+9B1f0aqAne+rOAtvjUjMldUuwVifJd3M9/jfo26d6dNedssYoFbB8liZwhlbO7sat1ngBkuobM02C2XAg37zna0F9tefgGuzSWFdxNxDOyRAWMNcKTsYNyHXjfIz0sXT8cskJfgekuLXjHsln7eChdbO52vqOQsoFKGlqCVdT0pY1okUs8DBu33Eqch6kgtJQLMKZYkTgM0I7XP5yxzwMsmNBbCbQaFnocwKOhBweAtF9yNLxeeRyIIDszL0HBQEIN5FfVM4NQC3OVGmcFx7mQgQSWHU3m1VJdQtmspmxsPCY/V+4JSqPx8szFeKTjN8i7w66uLxlo+5aDqVa6hPiQhP8JVKpRgp7iVZCe0jFkn1ORjQSLELTAQQug8L4TUrox7RpVkwULrPdF7hnrHIfFkVAim1vZuBrP6s+nTZzXXADuiyR4+9HDU2Uib9FEax3OEMu8xxgeu5/tVBaUHLBdmnwyvX38qTT7/GmaXXmXnuNrPLN6szar0/h7mvBARSPiK234MpyecxN4+O9USylKMS6y5hZLc4Fbm/2liI8YsYsKy6PpTEOJZRr+1r7rrenSXwmgeIFmW1EswiYA9DNIAe7Ny2e7nm2lMfpfy9BaIQZ1lJS+4c9Xl+yLkKm9bzU/ttwHGfQ2UKx7WXCQGRiC9t4kN6gTLqqY1F/YDxJf5YzcbTedVsas9ozJUk+kF2g/WpA5Lak9umzjOIH5jL2ev1Gb6I7bprx5PXcnX0qeaG+MirLgfl8X7m7/sqkrxuGXnwOsiWktqq99dXA1ACoVxquldxOR0OynDVei5I3ery7yfKFM7Y7nCGBbb5Fs9wgevc4Dxyj/bp8o1bq3zr1jPcWTvLd77yl1hY3GZQtNjfWjDrQDPjZYz03seiqsR9yDLw/AUYiPjkQg8mnjtQ+3LdXCUpzyvYjF3KWNaHQE4AIotgNR+7iuWXeDeSwGgZCwBoh6p7CfJDuEDJh7TzubPn7B59lvsmlsCo56P7Ut9knXhQW8WiudYwEt6T6b79ctspcb5TOK69TASI+OzUYXWnJEruS/uLMmxSihwoQQE8AFhmtUTRWDrHtlv+Q30fWLKeZJc5dlz+xzYLRxSmQK2+/cC5M+riw2SrRH5xxFWVuJhe2d86aNavrefji1TKfVUt7FjlZHrOBegBrdrvdLxcZ3clU/hnG8T0DJ7Klsgcu7xy6y28dP0H+LPXv4d2e0Bvfo/Hnt9gZiW5RXrzu3BtxniGK1RLbcjVom0ivsGUtchin5GuWfVFTJGDgYdyLfSML2LFEWexGli+vStUFbVeCtmVYtf1PJkta8SH6S4Bm7erOSOlBfIycNv6q+t6S8NHerXdZ3E5AgS54HT8VfddxLvOV9s+nLoEkF1wFv5QmcJx7WVCQIQj1ocXr0Dlm/ehsQCqgeULKiq81bfjXQlyM1XzK44qYEUm6Rxf5wrI+R7ViC3jJY6WbPffvSI/7r7rRHgqHz/nIGRQWmrpPHsuApe6lVW3qlKmuX+27fJ6VpzSnonPWPduQ3MlHrVGTpUp/LOdCXcqeUZbPE5vfo87+x2KokVn9oDdnTn6+x0Or57j1toyOy8+aT73Nkkhi9+Q8gMjxMGyrEW8ywIQOHjiWOcrRFeRXTtY6XeF6f4wZl1sUeVNXszH7uS+lZYDBh5Syh54xJ1AAikBSUm6n7PIsB1IivqQRKRnjkQuKE+eS7HrGrruldwf8UUr7jxltOt83H49X53rM9/1nE9zZcFUjmsvEwEiKnBQjw6SeFeSX/+j7tLxAKDz5LoRnwGWxe6jr+ogUM8OB3JORlXZm2vCXEd+Fj6gnYstenWf7mE3JzIOhtz3MH7GcyIeqKS8leHeJ5WqL0i5MCpzor55K0f3VA8h1vPTvjq35K8t7sRzWnVOadT6WbE12uthkj4drnOBXebK6L/ZxW0WH9tid2eOpcducLDfNYBYhCcufTspFSk5WRyyCjRbFjktRaz9CrPVDFyum8K16Wfn85jrzM/k5W5S+6uY+2yeBAI6VqT0Ggl81vK5S+4a3uoBu0cBkywaXJuLPRKAXMh8y3VbDKuIZqisYmCpvm5hICwXmqwPRXN5INE9+Oe1jrnDlGNSWmSvgYuMPE6mcVxLJgJEQMl4psSl6LrOSpD4cNxqORQL3ZV1IVdNNSFRvMhePreouJykAOvuNBHqShBLi0DZuiAeBNUvcTJe+fvFsEaNXPJAkgCoV+FAWrXnUOdFvKXjAQPMWum7+/B8Sd1KqgYTGKFfD1qQjFr6JJ6Bg9nRXg+LROBbt54pa4qtXV9lULQYFG02vvIfceGx19i8dZ7ObJ/ZlZslCNz8V08bUMjCaGOuKilbWShyBWlWLotBocGz7qU25RLaIini57As9nnSrF0WzjKmaMUH+DBdWRJqSwCmvJTCvTTj17uiwXRP7yG50laplhzhdr6/p9J5y8B8qFbX9ZFlsuIUDUbuk7csxNF4S2ANAzI9V/VNUWX7JADjBimL/niZxnHtZSIMqEio5FNIfGirgEQKTW6pOmmtCCS5x5TDIJeTLBLxKYPs+rL9B5XSIlL6dcUoxaxaWZ5n0MqD4k7qZU88We3vdRg3UlfYw1xiXhRmXO+TtVf9yVUyvv7svQUyPCigGiKszHl9F2CrLpR36x0nMUDRGnVec2fE495sCbz1sVfL0PPzF26wyHd5ee37oQ03Xj/P0mM3WP/mKuzMWHFDzzMo6ugKBhSyTpZJEViQlKP2F26bFOgatp7IKubSkfuGvG8lX3fV9UHk+2puZ81tk1Vw0R0vS0NWk9xM3uXVxlxGdY5CFo3uY5/UsSK7tQReK7X2Zt1n9UOurVm3XfesaxXufB9Crags7yYrAagAniUtLXO8TOe4Njn1zkII/zSE8FoI4Yrb9kQI4XMhhJfz++Nu3y+GEK6FEL4eQvjxUToRCeWM1ruiJJ7clVLzfnq5U3y1XqtBZaGqmg16XkDKzdYi6Zdrfeg4hWZ6i0JFFIdZEVYufg4tWStXlLgVr8w9WKZ+WYjtsOgmKWxvhXlLRa6sYe5B76rybR0XOl3/PGyNee8O9Jn41euePl+JITBot0d6jUMexNj2v80BXRb5Ln92/WLa2YadtSfZ7i8kAJGLRTNrP4PWzP5i3raCKepNrAAjVMuQeCWu5EQ/m34Oi2TyM/jn3bniThSmu08Ck0skd9ds/qzIJ690N7CorE1s8SxPdrfd8bJcdB0POnoQItMFmHpmCi32ZLo4DFap6QoAACAASURBVLUvl5/a9e6vNVwCIVVQkxWlZ1z26zbwdL2TFXnQ4/pByyjw+BvA+2vbPgZ8Psb4LPD5/J0QwjtIS4a+M5/z6yGEU301kaP5EV68O0tkubdChkUg+QKG6Zy9cnYNVfeXX85W1kc9qTHtq7qFdM1hCtt4m2oEkyr2ehCoHm8z/fozabmrHSeq02TlXSxyRADqz/fuQq/oD5xF47kn430sx6Se+Fl/vvXneJIMWq2RXidJCOGZEMK/DSF8LYTw1RDC3xpyzHuB/xr4LvA9IQStxTzmsR3Kz9uDhRTOu98pFdOT7/wLbm2ct1mwCO1L2CzcR2PJZeO5kheoFkD0EVRLGBAtY0q3cO+ehNeM3SffKSfCcyrkzwKMDVK9LnExAiUpcl0DquVTfH8U2SXXlq6hz1wHLsDOYQao22YdbLmXlP4G5o5T/8FqbS27tuu8yD7AoYUD+8TJEmhE9p/szoLxjGsYfWyHEG6FEF7Kr9PXGX8DciqIxBj/gFTxzMsHgU/lz58irV+m7Z+OMfZjjN8g/QTvPu0almxoYbtQdcMoy9yX1KgnCPp6UIrKkmJMGcKdcs7sZ9Fyt/goI+WktNFa7JZ7bRxCUqiqjSXRIlVe/D3V3Uy6F3+vdRdS3eKpWyE6PhV66ZRuM4Uqp8+9Sn9Fwqu+V5dqJJnPhNfL1/c6yG7CvZpFNirPUxdZpKO8TpEC+Dsxxu8jedg/mkGgLl8Afhz4sxjjL+VtYx3bdwh885vfw/rXn+XmtacZFG3OzB7AfOTM0uts35qHrZlkSexg64KsY7yD3D2rJIV5jaoj+r1UI7A0uxf4XMOUuUhuKUTvrhIJrYgtr/jFp4CBkZ/Vy0UmvsZbUApRVhCAlDIcBSgBo9pdAnZivvYFoAezOVdk8ZxZGiLPV/M5y9gSvRexvBC51Faxc/W8fKDBIsBMlYhXf0pQ69V+iOEyxnENo4/tP4wxPp9fvzRk/9jkXon1CzHGVwHye2a6eBr4ljtunWNsvRDCR0IIL4YQXtz9zuuVmbkP9/Wht3v0KgrWK9Buzg+RQpay1szdV5nt5UQ6SZph21K3YDkpil5K0U++NpQNHiPl267NAs99DAvhHRbWXOcgjrPO6jyJ71e9XV9WxT83yzVpM4yz8M9MOTBqx6LCvGU3PIx4GNczTCKhYtWd9DqxnRhfjTF+OX/eJjmtT/Y5mIx1bBcbt3li5TVmlm7DPhxeO8edzbMw26fVHrC/8QSswex78zxNXAiYe2YDm6lLSWvWvUlS0vLl15Pj5BIT/9DGeALN0pcw5SsORkr8JaoAomusYAUKr1BVxnKDqS1ZAN7SgKqrTNaTrAq1tQPMBtiPdp4ASIpf4KCwXLnRdI+brm2du4aBoKLMfB923Lk+bFrurDYkS6Q+vz4q4xrX8IbH9n2RcUdnhSHb4pBtxBg/GWO8HGO8fPbJOcd3VGfs3g3is8z12H1UlVemEiXQGUGvMF+b+acy3HNluG/VUjHFK3eQqv52a2CliKZh2e2qkyQZZm1U5yTFEStD91Pvl1fa3oUnC6MOHooY0/ZUj6tqqeg3sMKK9lz0ShFi3UpbXnwo9SiiAItRXsCSFHV+fWRYmyGEVeAHKJfEq8hfJa2X/rYQwjtP6d49je3w5OPcXH+KVrvISi1ZIKzPcrhhGdf7V54wYABTworIklKTopObR+4mHQeWMLiPEeviI/TuieI2tgLiNRJwCJgEVOtUixz6aCbxKrKYFBGmPuh6G26fwMvfjw9d9oBR5MevY7huz2PetS1OZdPtey6fI5AQoMpS0bPbcteUZbMIcFg1NgTWRcwHPkGZu3KM3I9xDaeP7RDCV0II/3qEsf2G5F6ZnOshhLfEGF8NIbyFFCwNaXg8445bAV45rTG5s6odOz0Rz8q3Wx7EAtvleXXgSPt3StTXWiC72cLRsR0XtgtHXUZ+FUTfz2EzCZ8oOUyG8Tm6Rj3E2D8bf94Ay8b39+GTBtVe3e0ElKHO6mc3570MI8O961BgKVDy3Ijccyp6OYo14gMsRpDNGOPlkw4IIcwD/wL4hRhjvUrel4G3kdTIF4DfIYXajHVsn2lFBkWb/Z0ezB9Ce8CdolVNYpPFoQzzVWwd9U2M63gBU7CaLYtknsXWEJE7aR1TqrJwPIEu4hiqfMgS1fXbvYvMH+8J/EXXtgcED1ZghSUFWvWw1n33kgW0GUla4jbsnwMupJUPZ4O5sJRE6KOsZrH1U2SVrWFL/0pkEen5yOqbJbnOdL8Cpg1Sf/zzOEHGPa5htLEdY9wJIXwAG9v3Re7VEvkM8OH8+cPA77rtHwohdEMIbyd1/I9HbbRTq+VU39crVyksSgCRa2qJG2UNK+/O8UlwysqGaoY7VJe1rc/29dnPwOsWAhwFO2+o+mivFoNKCZZhLiwpfJ/zUSfVqxZG6rM3jS1b/Wjyn7+uMtWNV6mG/NbrdVVzehTVdlBxKcoiUvrjqDIu33EIYYb0J/vnMcZ/Wd8fY7wdY1Rg6w4wE0JYYsxjO94JPLZyHdZmE/exMwv7XXPFFKRFpa6QFBy1902SCriKzeg9KS0+QrN8gYB3HSnhz0dfSckWtX3voerGAcslkWUh60D9UO5KO5+niC0Bh+c7NrC8FV0bdy9gFo2AjkAy+tyMfzFUQ3fVTz0P9c2792SpaJ+sPu86EyjLhSVgLw7tmev5FZDWf5/hNBkjJ3JXYzvG+FlsbN8XGSXE9zeBfwd8bwhhPYTw14FfBd4XQngZeF/+Tozxq8BvAX9KWoPtozHGU6eggTi0PEjbqS+welY+TFcKzy87KzfMUY7C8hcEIPqsTGJfsPAk5Vf/0WWR1En2+jmWhFgtG1I/zvapVEs1zFn7q3W/qtcVOKTAAbNC6q6wugUlV129Pd3zML5FojZTra29I/d9kozLdxxCCMA/Ab4WY/yHxxyz7Mb2JdL89IOMeWzf6be59eJydmXt5+ihYG4d8Qo+10PzBJHqcsF4kJAbx7u2PDEtpb9OUuptkgUAVt5EEOqtAZ0369pUuytUlfCs238RA4IV1496pJesAIGoAgrUL937vjsfoB2sDYGLXFhX3Xa5yNQXnyTpiX1vWflotjWOlk8B4ND6VgD7WkPkCX/QUBknJ3IXYzvkz+8m6flT1vC9dznVnRVj/Oljdv3oMcf/CvArb6RTBhq9coasgolV7qPKBdSjnoz/6Jfb+jnSSomC/rg6ObxHr7J9QKsSIjyo/fApJLZaoNCH9Na5lmHfj3P7+HszHmZ4mXxLgrSghDpPY9errozo3XPaV38+ilaruhQp3YDeraV1W+qFLIdJMvvv1cNakR8Cfgb4DyGEl/K2vwv8JYAY4yeAnwL+Msld9U3gf4wx/r/52PGN7fYdWDmEop2tkLR5/rnvsNN+0ma/mtGvUyXWNSOWe0UzYbmSBEJS3gqZ9W3ILaaZ9xZWYVf7xBeAgdWs+w5WD2sHq+x7Eat6K4BQpJh/r5PqAibdO1hYc90q2XHv+zFlqQtUBTw+ikouPLmgtlx7vlSMj8baoMrTyAKRZVNch/m3O+5kJic+3gTnVRgmYxzXMPrY/pshhIJUHfJDMcah/N04ZGx3Ng4R2e0Vo74rQ92SA81FJf+9lOfRxEE7v674rP3qTN1X661bHYlHmSv7W+dE6pFZXknX3Vb+Oj5TvpddcxKfuV/vv1lXRQkCHjjtuKLSv/p2SbUsy1FrxN+n+qF8E++S1GfPl5wkiYDsnHrcqe3E+EWGE+H+mH8E/KM3fLFTJJy5QyzasB5M8e/AzvqTtm6Fn9WDzY7FaUih61jN0i9SzX+Yx8hvKUe5ZmbdMStYPorcN1LA4mSuYMCka8hdpX4oz0Ig4PuxQTUUWGHJqyRuR4R3gQUEeJ4IzCIorZnMjQjU3PMsAUqg5TPg2+44AYUHa11Px+1TfX6bAG83ICrv+ZBUgPHkRanGNa5hssa2ZEJAJFbcHXUXj1/y1ifG1fkJr2i94u9V+IeinM23KEprw4hsc/9ou/qwy1wZ7TQsUVDuN7Xt96n/at9bPN7NlpTxUYvEan8NL4Piw3t95JkHVS+e7/A5MBJFePn1WPwfwVcP0H2on75vPu/nNIkMD054mCUWbWgXsDxjM/964p6U53OYO0UuK/EQ+qdK8Qlc1jjq05cCXcnnXsKKCGqGLXcSJHfXRcxtJqASuQ7mflvFFLN3P7Wx5Xj9OVLeF3Guq0yKg5Vz8VaR7sGD4SKwFUhembw0rlxmcmP5iC4BgPqtexUAKSdGkWI+Uk0ArmdAtn68FbgFiQ855DQ1Oo3j2suEgMhRkV+9rpy84q8SxNUFmzo1q6Yu/tgWAxbZKvMgvPKsA1o9Ksp/Py4fxK5ZDbX119Jn72Kzle+6DDjKnXgAqYvniOoWWrUCgJXVr7rJkh0i8QCSXIsqXulB0lx29eff659s8icZq9k/EdLqHDK3tMXOxhLlBPIaVdeVFOc+aa3xwr2k0AUUsjZEmIMpWilBMD5lnxR79p68vUIWY8S29qmNZcwSkJJWP7y7zHMoq5gLyVsZAku5q5aCKennsCKJs65d7wYjb9+KpEWpZuzeCqrAVrckZBF5V6Ai2Ap3HVlTi1QDA9okwPNWloCIp0jeojVOlukb114mpIqvWWf1Wba3KETZVmfA6df1RLXO69Jnjt3KLMASFo1obzFgi0UGtEqrRUvbqk1fBqTeT4msEFkEQEUx1wFI/fd99nkyEvWpHssBVsnYEist497X7FK7fq109VnSJa1WqH7IGqmLdzP678ru9/eh/f3u6eb8mDN7J0JCuENRtGC2b8AgS0ShtIsYn6BQWyleqEZlSUkrp8PnVFykCjpye/0YRkArgfBSPkdKViS1XFo+ekvKtB655COrxGPofjxfIYtCnIoPCrhSa0MupzawqRPyNgrgnLW57NrybcBRC24Ls9pkiW0x3MrD3QeYK1C/l9yApd46z0kyjePay4SASOJ86kRzv3QlVfkEr/xF5Jo1YkpeCrwaTXQ0T0I/4By7pXtM5d51jHIdekNCc71FkrYLSnzhwn4FuDRTT6SzKgsfjQarcx/e+lpgu3K/stYEnHVlXne7eUvP3FfVRaUECtpfB0VFcmndEr18uO/dyLT92SJnaLcHsDkLi/vGSVw8tKxvqNaaEqGrGbLCd6XkBCQCJYl3ySghUIpvDXNZPUc1ZFiuLj8D33HteGtEYa477jjv+vEuOhHubQysdI+ydDxYSTnPAuu+sGFM9bKYAb5hx+iaAgFtE3h67kSuPH2XlaRt4l685ePbEh8isCyf+4zrwPEybePay0SAiPDcK3UpRf8upVjPWPeZ7b5EfH2RKX2Wy0jrhJvrpe3aOSj7oT55ZdriKMkv4KjzJcqE1zHelSTFrXBjRaH56Ks6eHjArAOY7kPZ9FDNzi9q9wC4SsPVgey5G1lh9pv0h/wF7L4EuoKUkar4TuGMbfD6THJlLR5Ckfu9Drw0Y+GmmtXvY8S4FLBm5kqUEz8gV5bP4fDtrFENn/UKUoS9n7UrUkl5FTICCrddEUye6xC4CbwEgLIWFL21mF/XsMgpcTYegEqS/BymGbQ0IMCqWSFgIEXuxyoWJNDGgFqBCerzItUFpvTdR4+13asOdqU14lZbPEamcVx7mQhHneKo/cweqpnPYEo67UtqVpaBXxMEkguo7o7x65L7yKK01sMeWlujxYC9XAbFlwjxlo2tad4v+zLc9SOOoFpyXn2wVQOrC2/pXN93nV+NrDoKNJI6p6Q7l3LX/fXYqwxi70rTMbqWAU/XAYlxLd6aqoPwaRIJRwpXToXM9mFr1tb70Kz+KuYm0mxdCk2EM1iUFPm4NUwB10l1fRexvIpVohXh/lzeJstGfVrPr/fnfkj5+hk+rn8ipL3lIetlHXOviYgW+Czne5d7TNfxlogHuPJvkaOzfESYb7vIz0buuZV8HQGYdyXKAhN4+WTKi+756ffR/cm1pWcJjBKdNZXjOstEgMidTDz5ykwCEB8a23VWhpSVFHt9tT4pwuTjlxVQlOeqXEcHuaq8ck+rHgpctJ65F7Xh21U/65+9pSMRIA4vlVLd5l1ont/wbfnz9H342h/e0uiUgHKS26neroBj7pj4+AW2nVvs5PwXL3dZHuLhkDOk/JDFfbg2m5QaGOkrMFijOkOW8hI4yCIAU/re9aLjBCxqR64mnbeFAcgVrDyKB5Q1jDMQ4BWubfVf15BlpP7LotAMXm3VXUfz2JrusoY8ue+l7cht9YV8jriZJbd/BYtQE0jpHF3P8yHeGgIrIz9ba3fzMJVCKf+CpxdhnMpx7WQi3Fl3OFPO9KsuLblH6iVB+igc1ifWSbzSVZ6IFUdMSYTd2ox5WHa7v2Z9ll63HOpK0lsIAyxqy7t8AFTMsZ4s6RX2NgsoC79+DQMMKwLv++OV+DBF7p+3z5qt97P6G9gCYF0OWGC7BJZErnfL696NTKXZf4akoK5l7SvguMJRpVmQlKpIaOU0SLnNYtZLG6u8q3wIH9EkZe+TBjX79rklUrSKrhI4CDygun6IQEGAoMS+ekSZLKtr7vs8VYW/RYoaW8Gy6QU8QCpyCLDrPru+qy0R/XJB+b7qGfjILR+8oPuQtSKLSG5BAXqZ/FiLDOMwv46XqRzXTibCEomcydnpqZiiT8CTIrL1Qaz8uCnIojLzVTiruXCS+6mX1yNp5WPmXJtV68GytOs/rICo6yyheh2poxZBn3q+huWTpPYEEp4s17n19n2RxQGtsmTLsGRHiQc1ZZbXgVLney5ERSq9VeUXugIy/2JE/fB1W0b7g0xbPH1oDYgCiJdIeRRXMQtkBSu3voRVm5WbaB8jt8nnS1FLaa9ifIkUvL4LgKQIFQXmrQpdTy4kH/4rEJL1IgUq8FvPfVJ/wNxE3g12kVQjTG3oOlLu3u1F/rwYYEPVOtzy1HJ36X7Fkcii07V1zzu148T3XMQ4p1UsKEBRXwoaWHdtlQmfArWbjFI7a9rGtZcJAZHkWtllrpzNmiLeO6KYjuMe7LO14d1fqm1ls/SkaEXOV11FVaDylojnA3SOD+/t1cKKJT46SsUJ6/ukqI2LaZd91b35cxVIAImnENj23f35gpJ1nkjn+Kq9Vr6kmkjYYuA4qH5Zmj99r5Y1Sde5O2tkzOUhJkJa7QFFG1iNOVmO9K+7RIqG2sH4Ds2SNQR9rsOOtVnO6uW3l6uonnMiXmODKnegfA1tE9AsYhV2vfUhd9ciZgmtYhyH3FVyMy1RLbKo7VLgz2EWmSfFda7Ac0cd0EG30+ednlkK2iXwXKdqqYh7kZUhAFf/BF569rIO/e+g6+h7G1J+SMHolsh0jWsvE+HOkiVyQJc9epUZdcsp3rSt/m5rZvhV/NKMu7qe+nFrf/v1RsAvkWssjb7L/SQZVrlXq//puy+BUgcXta+w2LpZW11t0LLSvbtIjqxuDhUeoCiq6rrzqb3qNXSOzzT3a4fI0qsHI+jZ6bfwFodefpXFUWQazf6i34H5/VR0EZLylfLSbLpuAUhpKYpqnqR4BR4Ch8Jtk9Jvu+0+QdFHKynaawULKZbi9JaArBW5hhRt9RxVfkYuMYGC+uDDi3UvS6598TjzbptX6kBazVB5GAukMiOqqButqKOeyZI7X+AnYBHpL85Ez9tberK8fOiwf34lX9PL/TnnTh4u0ziuvUwEiABlhoFXwFBdv0LRVVphMJ3XrShOME6j6qYy3/6AFNUlxXqQZ+H1yrXeCpGiHzY7r+dfdOmXrjN/7XrorweCROLv1oZUUfIgAhNZQLJWfH/VvzoP03Pt1ger2vfRZuqP55b0y/j78c9I93GSBXaapCiWzkivh0VanUNmZg944tK384w4z1o16/Yz/CUMZCBZDOIRfDiqzpGyln9fyl7uJDC+QxZKHRTmMQvBWx9SlrJ2ZNkoaW/TteEjldQnEe2zWDTYvmtzA1vud91tI79fpOYn+QapzEhBxX20Hg0EJeq7rJWLmEWlvq1hpLzOrycX+kgutekTGtnLr9MtkWkb114mBkT6pPW6t5kvH2iam6cHq2xoMNeJXDlSqlJo3Zqyl7tKVYGhChA6dyuPolZNwfv26/uN96iGIqvP6u8wF5aU/lxeJ8UsoKLMRE/JfP2yHd2PoqNU+LCXo8k8P6Jn6kOfPfcxPDIs/XPnnNtK666pbz5vB7RIV6dyDYGrto2ypkgso/ROfz0sMjjMrsSiBQXML20lJXYJI7JVxwmqyXlSbmsYKMgCERjImrlIlWPwvMoSaR12+fxlBQh8lAioZEDNzAUMcpf5yDEp1jUs+mkNU7qyaq7k/XKh7WPL6Sp6TPfl3UhXSWT6PCQAWQU+R7nS4CKwElIJFe82W6RqQflnJZeWXFnqoyfhvcUmsJal5LmoMqz3NtQ4wrpM47j2MhEgEnNSkYDD51ODrSQIR4szdmvKyc+SfeiuuAvlMogbkEWiduXeMTfY0fwVKf50jX5lnwczr0gFbsPbs9InSjSsg1iBZfD7ZyTCXve56whIbfchvFXXVbs8p9qX5CIsqIKCd4/Zb9Eufzfdk2qb+b6M+geZOrP/MHC436HVHsAiOfEQWI7GjShxTgparzWqVoZm/22qri+5x3zI7Y77vIUlDErJysLx4b8CHSnZHXcsWH7KFlYsUdFZm65tn7MisBEYys1V51D0DGRhzZJqVu2DKem/nE5Su+Jz9t05G67P4m4kepZyvcny0NDccp/1rJeoVjFehmR5yI11DrIb/SSZunHtZJRFqZ4JIfzbEMLXQghfDSH8rbz9iRDC50IIL+f3x905vxhCuBZC+HoI4cdPvUYue1L310t5trBoKClc72rxcrROldWxart2wHiCFpZ858NoZYH469bXdLcVEr2VUHV5STww+P1SsN3auQINBRXU14X3gLlXRrfZcxDQ+Eg3KXofROCthGHck8C1zinVCzrq/nz9rbK/g1EtkQf3Z3sQYzsp0FkWWtuwBPPLmzn0NVhSnGbhchm9RHV2XGkLcxMtUp1py0oBU64i7tcx7kCuI08YFyQlr9n3DqaA1Y4Uvw851vHeAoCa0nXXw7WttjZc25AskAqZ/Wz+kBeA8oS+J+F9MqN3v6k/esb7bptAUcDnOR3vylO/N3IRyAqhrhWUh0vDiaTH+XdijN9Hiur+aAjhHcDHgM/HGJ8FPp+/k/d9CHgnKff110MIIz0d71oS9yHp55mx+BK/LrrEu2oEPN08T043YsR2N7Mw+i5GJvWhCliacYsX8PkYupb6Xo3+Mh5l4K7tz9Fx3Ro4CTjqz0fvyo05mu/SKp+XrArdk+eXrB/GFek869NB5Xn06dBxUVk+eMDfk2yw1E4i6XdbLkTzGHkT/mz3f2wPgB341vVnmF2+yc6/ejIprkuHpuBWMBJ6n1TJ1xO5UoB1t4rnMaCqWGV1CBwUobXm2pRil7UioBC5DRYWK/Lcu43kchOYeOWrMFsBpcCmwEJrdQ2fx7FEKrsuS6iNexBzqQFdW0peBHrh3n2Awo67rq6xiQHIIgnU9Ex1/4o62wc4zACpUiy3SaB2G/hPOEkeeRCJMb4aY/xy/rwNfI1UGe2DwKfyYZ8CfjJ//iDw6RhjP8b4DRKl9u7TrlNfATBtO6gdk777qrZzZUlyK3ooziAp46Prfwxoly4cHS+SWoAhsRS8ouRVqvtsxu7zV+xaBoY+J9+70EDuNou2Us0t3ZNvz79rTXkvVVA6GijgLTO91yPOBBxamTBZf7bcre5Z12hhC3XJMukOAeaTROUhRnmNQx7I2A7A0iF3Ns6mQoyr+awvZnJYs1zxClex7S9hRLRemk1LqfvcClkH2i4wketGyYhQTZgTeCjqyoezCmAuUg2VFXCsurZEZq9SBYlV1z9ZHuuYBdJ2ber66mOhh6gGzqWS8Dpn3x0rjkNAIhAQD6Qm9Bz8c5MLzwcUeFfc/Azsy/txCHybhMjnGYUTeZDj+kHLXXEiIYRV4AeAPwIuxBhfhfRnJBXXh/Qn/JY7bR0rx+nb+kgI4cUQwot733kd8GuGmNvF6lIVzr3TLal3lSRRlJZAxbKmzaXio6d0HW+5qGKtF7mEqpFfcj8dVIhvibdE7A7MxeaPkQvsINtMksIBjT+n+mwMSAR13vqpJiW2K/kn/pkMCz3u5fwc9SW1bAtqHX1OyXLpl0CawqYFhu3B6GVP3owZ2/0a29z8Tto4CzvXnjQldhFb50JKepbEkUj5X8ZWGlSJEM3uPZ+hNkVMw9Gqv4rSUsivFL5CfAUu4iTEc8iK8MSzZvPeZSSroZ7rorZlmQiIVjELAWzGv465x3DXbQeSJbIO7FVBThFeatODo/rugeGaa1uhvJ6Ub1N1fZVWkSaQbVJ47x5wnVHLnjyylogkhDAP/AvgF2KMJ1UcG7Z045H1fWOMn4wxXo4xXu4+uVA+xjpxqxLjYDNkRSAd0HEOp/QD7OXVB+uzZLAwYbN4LO+iRTWvROfaZ8sM9/vq67T7vns3l5+9e5DxbdXfAXaZK5W8V972vPqV+1C7/fxs6ln0+tyu9KON5zVkvQkcVaDxgG5eZ8XXxTKQtn4V1PmQzv7JYZC+b+P4s4UQ3p95i2shhI8N2R9CCP973n8F+Cz3YWzzxJM8+bZXoR2ZWbmdCHVZHF/AeIyt/PkKSfEJZER4K1dkCVPYPgRVloe2eX/+JlbRV4pU7UqRKmIKqlFaug5Uw4qVPLiGWQyyBARQakPg4UNli3w/Kxho+fYFKh4oAZgzZa778VaV54J0b3JfLWNLCq+76yrI4BoWmeZzbySVoT6TXzc5DURgvMT6XY7tfx9CeNdIDd+jjAQiIYQZEoD88xjjv8ybr4cQ3pL3vwVjl9aBZ9zpK8Aro1zHFKT9Wgod1Yx2m4XSbVQnjD0QePfYsPWND7Id4XMx0nazRnx/Wgz4Y70p9QAAIABJREFULot5hl44l5ot3uRf/vw6kNTdQf6ztwj02ZeXr7fjEzPr5rBfqAsSIMlqGW5NaMGvomy7Wj2gun6LztG7t6T8vRWtFoP26UNN1ZxHeZ0kmaf4NeAngHcAP535DC8/QWJsv4+kcs7er7G9dWMRrqUore/7nv/PlOsqpug3SItHySLxvISUtvIcfGa1FC/uWB95JWUu19NKPq4+O38vlq9xFbMAFA5cd6u13fY1qotreaUO5jYSya4+6j4FFLJCfM5GRdr2Vg9DFjgIGLyFoSAEue3kclM/wEBRVpUsEz2rSpDDNskqeppkpB4xRisyrnENdz22nwU+AvzjUxt+AzJKdFYA/gnwtRjjP3S7PgN8OH/+MPC7bvuHQgjdEMLbSTfyx6ddR35zc/8oYbBd4Rp67JUuK1kGXqkld01S7GkG3XauHgsb1rna1s28SD1EV/ka+gyWryFRf8XF2P1Y8t8wwLDvqYfK+VCxSF1r2ODybQgsPDkP1fIlHhCl8Ie3W5gLCnP3idvwFQL8tXyf/PUE9L3X7xy5Vl3GGE//buBajPHPY4wHwKdJfIaXDwL/jDS2XwD2BRyMeWwfbpyDixGKFl/76rtMoa+S6k7JvbOO+eY9N6LktksYMKxhvABUQ4OvuH3iCsRFyCIAs05mMRePFKhvD4zEl0L2kWCyhjw3IuCSm03kt3JflFfi+6O2fELfLCRj7zZJec8Y6ErpC9CUJyPrRiChtgR6spg86KgNWU716LiKFXKOxIVcJ2XUn14Kfox5IiOP7ZjkBWDRje2xyyiWyA8BPwP8SAjhpfz6APCrwPtCCC8D78vfiTF+Ffgt4E+B3wM+GmM81SGuWe9RDqBw5Hl1NT/54aXEkyK3AoECFO+uEXEugnxYXoksGV1D3z0IgRHgdUDROT7c1+7neHcWVLPL64UOfUiyRNZDvaCjWS5HC0X6a/rKyeq/v29viaQ+dI/ck/q154IOeuwZEPYPCFWv17FyF2b/kniH/PqIa2YU7uJpkhr5GeBHgLcCXxj72L5DIsgBdmY4s5T4vzI7+yqmwK5iWdzLmDIVkKxjyt6H98qXLxfWJUwZ1yOxpGxlrfj8EVkpYImQnrcAs5xEWAsQdtx3HVcnv69glsdLVMN4V6gCpobqPiQSew64QSW0diO/CyzrEW1qR2AgwNGz3aodL2tFoLiFWU/PuT5xHdglgcnJpLpkTOMaRh/bp3J345JToS/G+EWG+4IBfvSYc34F+JVROxEJpcKBaikRuY8UJbTLXKmgPO/QrvwUvuhgVXEKWNJ3K8ioHJF2hhgVJvRuJrnIDugemfW33XFmFXWP5JWk7UZuq7/1+099s6iuYZaMv1eBMFBxK6W7mSvvvYNZffX26t8FzQm420f6JX5E26vgliv+9ndHskIgjYNhrsdjZDPGePmYfaNwFwH4DzHGABBC+DzwP8cYv5T3j2VsW/Z0gGtw5+rZtO2HD+H3ZozYhkSkr2HuK+9GkkjZ+cgiSEp7BVOCnvxWJJd4l5V8nf1ae96lcxWLxFI/fJLfMlbc0Lu21AcfTiyg0fMAA6wd1ybuWWxBAosZEoDczAeqjlbO11A7HgB0Pd07GFj4XJLF3EefdOjDpAU2/jcCkvUhILnJKCsbjmlcw+hj+7RjxiYTkbEukVK0zyKPU/7FvHMt1RWazhnQyol3SbmKePY8hlxUoJpQ/bLkiVxQC2yXfIzm9r4KsNw8dX7DLBGzeupiiYnVzPEWlmB4nNWiz1b3q1sDTbNutM0staqy9xaapMgAo/L1Slr04Ov7Xs/Yb5MsqPLe2q20KuwIlsgYfcejcBf3zN3dlWg9kTWSIhI5vj5jyvMixkf4PItlqiv9iXCXW0khsluYwn4BU5w7+Zj3YoqxIFkE4jCkTD0IeULdWwSexPa8hcDGk9Fyy+l8hdXq/rzbylsL4isAZmdIrqI59y4JwK7ldej5eHJfVpMH1Xp2vSwyWYZtqoCmPJrSurmdrlsWhrzAKLWzxsWJMEljO8tEgEggVpSwFJMAQN9v5DhEFSv0vINWIgRz3Uh1eoLZu3U69EuF18t1qGQlJNDZKwEl7TdS2/MC3jJQu/XILv9eF3E4JyUYetASEHgg8cf4Nnxl4+P64ZMGfaa7LDE9Z1/pt95HPwHw/Etn/5CZVJf+VBmj7/hPgGdDCG8PIXRICYKfqR3zGeC/y5Es7wFuKax3rBJJiugLJEX2ApYjofwEgL+GzZBFMEuJXcWIYilHfZZbawtTlFCtW7WG1efyeSkCFCligcEqBiSeIJ91ffMhtlB1r+H6p30iwr1yF8jUw3nX3DGcI1kjT+Sd1+3Bzveq19IzEFCVyYsYCAo8r2FWicDYt6W+eittFtLa7z5x9vSowzFzIpMztrOM1Ov7LZFQzm5leagESI+9srigXwd9QDtbC91sfcwhHgJUIr5TupSkAKVgda5fV10gILeVXwjLAMIiwrT9JHAAc1/Vl6sV8Q2Ua3oIHHz7Ahd9T4DTzcvQ+rpYZnl0nSXha2MZABmwDnNt+TIm2r/NQvkMvGUjq80XiGwPBnT2D2kVENsQTs81LPv3RiXGWIQQfh74faAF/NMY41dDCD+X93+CFNL7AZI62QV+9g1feJgckJTQe/P3FzClKutDFoPyOZZJCzgJOJQvohm2ckd81rh3CSm6Sm6nNjbTlqJdx5TkvtsH1TVJpJSlVAUU4mpWa32Darjsfv78HNU1zb3Sl0Uga2UVeDGm3JDikMQ7KKT2BvAy8LStK6LgAd+GrA4PCjvuGPVN9/y8++6t5iuYW6vw+0Smj8aLjGNcw4SN7SwTASLJEhExbqXOB1jJdisD3y1nzAd0WeS7vMYF5vO63pbv0EZrtCegSTkn3now90xSmF367Ob1TBKX0C7zVFoUFTeZsrOlSDVbFz8DSZkLoCRWhsTcQZ689u0IQPwMRdyHrCJZH306ZTSaz6z3HI0HPSl733Y7PyMf5FAl2Kv8kp61riGOpNvPQFvAzG2gS9Wvf4yMcy3qGONnSX8mv+0T7nMEPjqWi50kM1gYrJTlFZLSkqtJyk3Wh4jmK6RiLC9SdS3JbSO/vQ//lfUgwl2Kb4dkibTdPoHYJslNJtfNSt7nXVZrGLBIoYpn8Upa9yCQEJDp3kXi72B8hqLGtH8HQKVPZqA9A8U6prDfAaqdJ/eZriuX3Jr7DBbR5SO3dI6eP5ilNO+OkdttFih2SZbIHFbgwPLHhsm411ifmLGdZWLcWd6og6Rc59hlm4Vydu4iGFhkC3EVVSVZUFe6IsGV6Z5WUUxKdpuF0oVWLb5YZACxmlNSpipGaMS2uXPMmqkmTdYJ/2HiF5Ayhd8ur9djr7yur+gLZFC1jP56m5JqHoyVJxHns0u1xpXqgfUdJ6Ln49t/fJCmdAv9bVrFHTr7d5h5He7mvzNm3/FkiHffQFLWWyQl/WMkRbVKNanwJfddnIFmw941IwtDyk6uJinKS/maUtg+N6Qe5ioQkKtNYCGrQ4DQxlxwElkly64tf/+6HlSXnlWY7xrmVpPLid0MYLezNfI0yXV0Pr9vU3Ft+ZwZWU8bGECt5Ofho7hkkejZ4vqG6/NWXqJ3fzc3fDt3eiEfcHrZk6kb104mAkSgmskNNtNVbkaHA+by6hgCDoX79vIMWFFEHnBkDajkh0qc1xW6OAO50NI2m7kra1t98S4kH42l/fWMcUk3k9b++n69dz2LDgfZ2umjMFzdl6wTFUT0ZfLrpVuqwQHV3A31X6siigNSuLCR/BYkIDec+q5CmLutHt3+AYN2i92zs8aDDEj/u7Onj4EUxdId6fXQyICksOSmmcVyQ+R28u4szYjXsOKCqyQlu4i5uKQo5baRG0whqYuuLQGLlLksEbDaW4vus0TKdRZLgtRsfomq4lVfNbOXQlcUFxjAKMBgDSP25ZZayW3N9jJYKfJpm8SLLLjGcjWaxWCWgywyEeYr+Xpr+TS5u/ScfISa3IPqu57Z7HmS5dODeVlDet22fhwjUzmunUwEiIhYVz0sm+WmqCm5d3xy4W7mQCSt7MIxYnxQbktuqkQwG0FuoAFKZNwtj5UiPcgABKZ866v+pX398h40u/dEtFxwUvii/XWfHowEXJ7H8RWHdW7bKXa7ztEcFE/O111rvsS7+llQXSXSt+fLnnj3o2Ru55De6/vsnyW5seRief2UQcCU1hg6gyUHrmGWySapDrAUuPbNYgS4J6cvYUrNk8QCGUge8DVMKftw1jWsZpcIZ4Xzglk8qxiAgXEwslLkKtIMX2CkKDNPbEPV5abjBEbibQRscsttYcUO22Az/ddI1sdM3hbMchAA6j7U7oZ75gK3LarWoSLg1qnmysjCKpMOr8POIQk0niDlrwB8k5NkKse1k4kAEbAFjiReaS2wXbpVRLxvs1BxM2lJWn8u+BDU1LYUsQ8NlnL1me3DFKmilQRqw4ov+rpS6of65QGy4+Yf9eM9GGh4HWTiWqS6Sj+qvXrJE1kSw7gatWt9MhfiAQrrrSZC+oq8OhYsUKHFgN1uj+3HZpl9HWZfT4R67szIMnV/tjskJfQeTOlexABASnMFI8efo7qok9wyPtxWM3ofdisrxUdK6bNAQzN+z2XIipClIpCT20ptyU2mCCewmb93sQlwBFCycmSZ6FglKCrEdhkDInoYeS3rw2eGr+RrHVZdVAIT3RNu26q7F/E/AhdvBcqNpWe7Q7KMuAD8exKKaonem5yWsQ5TOK6dTASIBCwZrZ54B2QFajN9zdKl/OQ68hnb3o0kd009a1tEu3ebKZxXrh3BisCjVeub8irEL+izjzQDKu8+VFhKvh4h5fmVOsnu781bBVo9pcp5HE0E9EmMvg86Xs/SV1XWU6gvHKaIsvK3GwySFZK7ePgEyaUzojtr6nzHM1h2tngIKd16dJQUuo8QkuvIk9o6VyCx5Y4Bsx4EHGtYKRW5nLyob35mLrCQpbTjPut8T6DPu2PkShJY+cgtkdQ+Wkp8itrchxQbfS4/h20SkFyksr56sZu++8ipZdemgFj9FijIomu747xF48OW9UzKwJBnScCxC7yLZBGtcpJM5bh2MhHRWWdcMqWiizxHIrdUXYFLkclNlKyCfhmlJWtBM3e/uJJAaIHt8sdbzFMn76LRDN6KPqb269naddEd+Hbqyjy1d1A5X1aVQpe1TWDnj/MEvC9H4rPX9Rx8DkqVE/FrzVejtST+2Vm+SgKzHntsM0/JJbWgk2PnixaJG2lxV3kiUyVSbp7UlgJVTod4DTAXli/VsYQpelkWItVXsdmyIqakiK9hgHAlt+fdYQqNlRUka+Yq1cir0p2DgYsn/XUf3uL00Vo77tw1zCLQ/clC8dbVUjAw5BwUkaSwnyC5tBZIFX2phg3LzSaAVv89T1O47R4EBeL6jdYx16IsNXzI8bexSr7Hy1SOaycTYonE0i2k8u4SEb+amcsi8QDT4YBFvlueZ5FQFhLslXlyj5lrRpFe2q+yKgIab+GIN6hzDYos85aC7+dcDp31oOM/H2SSHCgBxPcxHV+Uz2mYywso+2DPb3jxRt+OAUnqu+7F3GG2/oq+e66ktLLyErit/CdtD+DwLrhClYcY5fXQyIAUxrtPItTBlJgnbz3AaHYsZS9SHaorD4r30Dki631o7hXMKpAV4S2aRarXlGtNCtZzDQILH80lLkF9kQLfd59XsCgyT9AruVHfr1INq9WzKhQVtYZlrWdOBCwqzGfT49r3Fg8YWOp4WSty4wngdO+6nzZY5vwetqbI6RnrUzeunUwEiJzhTmV2DJZPAeaykcvpPJtOUafj5MsX6byD1iix44x/sLU35HJS9roUpZaA7XLA44Otcrt3o0l599ijvjLZUbDwC1sNKvsVUGD7ilOtC2/ReKAQaCohM7V/ULlfD6p+DXkvcmd1SstuuM9WUXMARauVI7Rsf1vRSSOAyVSa/R2MsJXyXac6W1aoK5jFoGfoI6skP4y5vDR7F5hImXv3kXJKfG6Ezr2KgYkHJL3WsbL1Eh2/WGtTwKOQYrnavNtL7XpLR5aWhqBCjEtRJNQqiVy/QJnBXuwa9+FBWFaWnoMPPVa/BIxyYan/AqQtzLJB55/LH9SnNtVyLEdlKse1k4kAEU+MS+GqIqy2K8R2ka2KYk35E0mJzbFbWc5VbfmZuWbVu/RKV5KKOkICrPOv3yzhaECLfisp5rkaca+cE9WVqi9J26op7vS5qMzsdbyW5/VRX55yU2+GBQboKF+1WNfynIi5wKqRaR6Q6hwOcKRP4nV8lWH/effsrJHqAH1GJtfHWB5iMmQb4yUUqSQXiaKHNCsXSQ3VdUakEBUdJRDwICTLQLN9EcVS1ipnIstEM3HN4kWCS6kqSVCurk13vHep+X6A5arsu/ZWsUACnaP+6v69FbOazy1dZedJyvvbpMio6B5Uz8DjIsYh6fx991mgJUsFqjWzFNiwhYG6AL4NlqfSJgHIdu7TKMT6lI1rJxMBIme4wyJbrkZVv5zhKrQ1KfsEJzsslNnZm5yvKeNuqWh9hrpEyXudHD4sYFAIcIc+22fnOaDDQn+7jKhqDxIgeXLdh/96XkaJit6FBmaRWJKfkeneElGV4HqlYJ+34d1QVdJc7ruivI7vx9FVFqvHptwT64uuu8BOBczUnrVV0O0f0O926L2+T9GCfrY+4llGYt+mMhTyLFQq866TuAol1ylCSDkez2GulfdiClrHSPmL75CSVNiuZtiXMeWnsFWF4UqhF+7lLRHytTSTl+WgiCvdj5S8FLS3mrxyFqDJgtnAAEX3I2vEu+pkxZTuotV88byuCK4UvCwwgaJ3a+naq+6etE8ckAcQPXO1hZ6LCkI+QQKzQ6ym1/EylePayUSASCCWXsH5HM7rOZIDuplYV2kUS8Trsce3eKac3VddO50j0Vx+v1w1e/TSzLqfrBj59ne7yRpqDwYs3DostytCqZ7Nrhm79useVNrEVxQGyuTJLrbgU7KKEn/il9717qzUZlXR+6FYn9H4sOJOxeJLsCAXYD3cuJ6Q6Z+rorJ8JNh2d6E8ZtCGQfsMRSsR7PsjRmdN3Z/tDMnykL9d1gdUq+d6F4qS/NaB/5ZqToZmzKtYBJYsE++OUUjrKtWaUOIqlt3xihpTHwUoni+RC06Ao+PlovKuMNEX6669TRKArWFRYx5E1C8wgCsthpsk62Od5MrKhLqAZB5L1hQ4yd1Wz0+R5Ya7HwGnIrdW3Hc9I20rLe5zqR+co5pcclSmclw7GWVlw9kQwh+HEL4SQvhqCOHv5+1PhBA+F0J4Ob8/7s75xby+79dDCD9+2jXaFFzgOm/lFS7wGhe4ToeDXKtqgNxdShL0BHSLgrfySlkcUJaFzpeLSO6eTg71BVuhUIr4eveppOJb3tVT0CoGhAJahRH2uo6IfR+iqxm/TywUSPrSLnLZGeGeFLr1z9xe9RUOvfXhP+u8rstJ8RZUvSaWL2jpgwvqVpRqhQnM5F70x+keds+mf2mruMOgDQezZxLZPoI8SN/xgxjb7JNqX61hLqw2VvDPv7x7SBFH/4aUYyLy2YOIkhillDXTB5tVq933UlWscoFJfMitSHDPt0ihymJYcW14ol5tX8zt6vhlqi4r3ct+7bP6Xs5dImnGX5AU9nVKEmQemJ8xl5/uQ89uEytU6fkbT7zDUQtljeqKiyLn25BAjNyPgkSui5Q5Xh51TqQP/EiM8ftJQ//9ubzwx4DPxxifBT6fv5PX+/0Q8E5STu6v53WBT5T6+hSKwoKkyOdyHSsBghSxZs9L3GCXuVIxK+zUXEGm0JPbzHz4c7mCrqKy2oNByXdAWhMDoLN/J/fHIrd0nMhugYrcVb3M00haDCrfBSy9bJVIrHijRUF5Yl39rpP0Xvw9ygpSW8YHJQtP1/a1y/w9paKOcxXg8NYMVOtyCUi0tvrB7OlD7Q5neMDlIe7/2D6DzWRFUK9SLT0iC0W8hJScZucv5HMU2XQRU4iyKsBIeylEWRmeTJfCl/Wj2f9W7bP6JYJc5yq6a2vItfYxF5qsFB13DVtzRK4lz8N4AALYVF2sgqS4ZzBSewaKb1h+h84V1yPX2vOYReeJdW9F6b7kKhR4ynKBqmVDzA9F66zfxBbKGi5vwrh+oHLqPzuv06thrXrMkbSO76fy9k8BP5k/fxD4dIyxH2P8Bmn4vPuka4ScJ6LZ/CJbdDhgge3yJYXmF22SW0ZJdirvnjgFK5VS0Cqr8WrW7TmQA7r02CtLqRStpBgXBtt0+we0igEyNAa0U5FBp7gX2SqLNdYBQlJfDdAXcBQYSYzSt1DheqJji6LCf9S5EX2Wq8oXVlTbfk17AYv2iyMZuF568PRuvPrCYJLts/P0u51skYw2y3qQZv+DGNuc4ahLREp5HVNwPsxXs2Epw2WstLp3JUkxK5ppiaovX9sEBj6qSiR9AaVLSG3K8mhTTVAUuS63F1SJcOWt+MgxqEZJ7bgXGMfilTbA/AXMCoHEReyRgOR2epdb8DIGgFANMvChy7hnoOsJIOSy8/kk9TyXTbBFA99F4kNOzhGRPNLuLIAQQiuE8BIpvu5zMcY/Ai5ooZP8ripkd72+bySU/IZyOM6zWVH6A1rssFAi9pFMaaf4VCZFBRpVpdeT2QsYEPiIqQ59Fm7ts3hrJ7m1ipxN303hqi0KtrqL9Pq75fxBGfMDbFVFv16H5xX26FW4lPrsQyBULwF/kHkIb40YD+RzU8zaUAl9WTo+PLme6S4w8GCkTH4PSIUb7IPa54IU6XZAl/ZgUD6H7bP+H37yOHjQf7b7PbY5g0UhSZFJgUn5yU0Exm14a2TLbVvjaCKhlLbO0TF1glkK+yrmihK3IMV+2fVVICaLZ5VquRCoJumtYuVLoMr1eBeWIqhkwSgKTEq7DFMOudyIop+Um3EB5s/b81F0F1gEmZ6H3Ghq07uyrlENLhBA6rno5d1eQPrZtSzuKqdFZz3ynAhAjHEQY3ye9FO/O4Tw3AmHj7S+bwjhI1qM/tZ3zN+vfI8uKYGwXptKjztlSi9UFkqC6ixes+dFtso2kwtoWORSyuO4wVLpx3/85j6D9hnO3rjDYRd2z56h9/p+eY4I/nR+taaUAMoy3Q9KsPJuIiX0+ZyQYSKLIYU1W+SWCHswLqRuAdXbldVUT47UdexzUV5DItCUKFRaVoqk3zLwb1HQb51uqqfAzQfrO77fY5tb37FZuABDBDqYVSKA2MQWQlrFXCsSuW6kKH0mus8F0U/qCfv3YKAgELpKUrAXqa7gt5O36Xxq/RDwzbrzpLRFrstykiL34kOEFYXlgw+8C43zJABZy6+5o/cpS+pq7dnWEyB1TVlSO26bwEa/lzgnkfZl7orK0bcZbY31hhMpJca4BXyB5A++HkJ4C0B+fy0fts4I6/vGGD8ZY7wcY7z82JOdcgavma+U3Ft4pVR6kJT1Dc7ndUbEO6SVD8WjyPrwM2+VNPHKOpHuaeXEXt/WX28V0Hv9DkUr8SFxNpfvIPn4WwwqkUiee/Di+Q9xB37Gr36adZCewxyycvoVi8a3K4LbnFnt0qIQb+GtCs+fyLqTS9AnIqqvkj629kqd2G9h5fgl/pnrOL+U8Mky1mVE70ru19hm4clqtBWYolaYr+dDLmF+/Q3MjSUlJr+9yOd2PgbMmtGMWspvLR+35tqVYvWJfrIuLuX3NuZu88S3Hr8si0WqhLvPENc5K7ndNdcvuffkcpoFisOq9aC2Z0Von0+d9daDDzi4iJVkWaNqaYBZGuKkZDXNuuMUTHDN3atAml0sS/0myTI6bYncN29cPwgZJTrryRDCYv48R1pK5yppHd8P58M+DPxu/vwZ4EMhhG4I4e2kimV/PEpn9jK9DEn5zLNdRmhJOe3SK5eFhUSY7zJX5pPIVaTZuHJC1K53W6XtcxU+orzvPOjmdtIAUfJc6d7KIu7FK06zcixEVtfzloK4Dp2vGb2UrgcFHSOFr7XhLfFvt2K1iSzv0C8DCXyIscKKU98OKucJxM3iMqvHl0KxflspGs+X1N1jp8mDNvsfyNi+g0UMLbntCuWVm0cuqxcwf7zCf6X0xUn4Wfs+tthS4c4Vee2z19tuP9hs3GeaF1RBQqCl42SleCJ+073WXVu6BynmLaq8jUKUPZDNzsDGoZ1TkvmHJCJ7gbRQ1Y2ceJk/e45H15mlGjilfVu1Y2XxLFHle7xbryABHFCuLcK3SQT76YtSTbM7axToewvwqRyFcob/v71zj5Esu+/658yt7qqp7ppuzwy74/VY9EprsbFXYk1GDsggiB+JMQEnkQw2IrKcSAbJQUFCIl6sICuWRXgE8QcEWEwkS4kwK0Fkyza21yaRZYRJbHCMH2s8ZEd2e3Z2MjPbvdXdU9Vdtw9/nPO953erH1Pz6O7q2+crlep1H+feOvX7nt8bnvHef8o59z+AZ5xzvwB8H3gnQOz3+wzwbcKtf7/3fl8p4vCV0FxkhdMxAus2XeajALzJuco08hKLlQYhc1K9XW0I65Wp6ybnK01EjK9cDKHf7tXMU4O5UM58MEfo0NcK2kl/YYbe+hr9ufnajy7ntDUP7baykAC+HfueK/TXElFyWlsz06imRYWe88E0J+2lT4/TpDyUEFBwuxaeq7FajcWeW8UWS+r9TEYU9Mo+G0W3+m3U6Os0o8r/Mk4g+nwSIvG4mqnsoOCc+2fAXyXM54eccxI143P7w4SQiuedc2/y3l+627nNJnXNQCt+OcufI63IZUKyDneF79pQW+ureBKqKr162PciGWsCU76GorZETBdJJVrkB5CAXTTHg7oj2+Zd2MivZTNWq7WsmePY0GIR2NIMVcb4InDtRWg9HIX4ZSq/CM9D59GkAVkixIzbksQ1s520q475HfTbyLezHI+xIgKRc/8WgdS+Z064O45D9hosAAAgAElEQVRgXm8C/w94b9Swx7e7QmC/Ehh57y/dz3nvSCLe+28Ar9/l85vAm/fY5yPARyYdhKKztHIO0VIb2NIjkEworYoQkkknOJr7VcSSFYwyudgKvMFpH/4J6k1S+VoWOlVjpbJ1isHCNp1b4IYw295iszNDSUG33OBGca4ir/B8OxJiaidrNQhNJo0laU71YouWkFTG3gYSqFJx20SEKYpN0H2w/UGsfyNpZ1YbGlZ+DwUviHxeKhbH9g/kpLBqaYcytSkQYNJV1iFWO30WeMp7P3LO/RMA7/0vV+OIczv+2S5572+Y7+5qbuNJZh3r6BWZSMiNR2tBcjQrumjR7CetQlqKtpV/Qivxx0mC+xqBQKzz+3FCqXo8dFwa03OE7Hk5rWWOs0EAGqO0CEWASYORb2PZjFFmJ634RSD2ePNm/2sAvaBtcJFgOlLIbzRxieCsGe1Jkn9E0Vf23kjD0r18gkQ2iji7Zu55Fbh3DnieQCAxSmwiTeRI5vVTwC/vse2P23l9P5iKjHWPq1arJa1QAqUMbvPz5U2GzFaVdiXYVAJF+7Tj6t4mG9rQ1zKuyiWsFc4rAW97fLTKskqOq0JTC70/xagIQn9YtE3fkRK1xd1t5a0ILlubSmYiG0wQTqUcmHppdmuyUuizdWi3zH56r/PK9CWIuPTaPqtUvj2uEgwFEdD4Cms8s95+PgkOQ+333n/ee6/l41e4U8rx/WCGeh0pCU7lT5wn2d4VQXSeFAklzeMiSWuwq23Mftp+iXqYq4S5BKnMRNbPgUv+B0sCIggRj42CkgYiB7rGJ7+LCEf724guHU+mNEso14AbWvl7UoOqFqmzISFySwSr81k/iQixY84lf5PVRBTeOzLbimAq5368R9WG8omoCOP+aNy8NpgKEgGqZDuF4vaLHrMMGRVFJbwkiDZimCyktray669EU5fCfqWFzMdnqzVYoWrJZnawVWVZn17bCglzLaAIPpFWGRISu8ONytehCVAvXHibNqmCri0MKUKxjnX5dFTVdxTNdyoquUG3EtIyZQm7CW5LkDaXY0i7Ov9uWfAlLeOwr+eDjMa0LEs88pcMa+fS1UwSnXVXtuPzioCKj/fd8QS74+eB/7rnkODzzrmv3fPxHXXzDdRzNuQgl0BeIjnZJVD1/knqjm4JSmkIOoeNspKWY6OuRDoSvFqNr5ntdV4JU5GD5OUV6p0PrdFEhCUiE+EsUy84KUIbkQIKRgA3w3FauoEQVvwqNQKwlUhRorJljiNCE+RzkpYjDc+a6xThZYnSBhwAKTekT1WCpcpi3x2NnNcGUxMOMKKoopL69EL59VFJd32bjblT3GiHuHBbnkNmHgjRVyFiK62wAXpln34RfCMSkrvZJ2umJIX1zlGVNfetYM4SZgeBXEIy4jbMpdwNJRPquqyGISILJrRUZdfmvUjgijz79FCuh46jfBPr/LbCX6asarxR01KIsDAeGBBMWalPiPbRPRsPjU73r2VeJ+IRLOHtBY+j3J54NXZjP1uuc+4LJKONxQe995+I23yQICZ+e4/DvNF7f9U59xDwrHPuOe/9lyYdIBA0WJl0IJDEMmn1a+32LZJAl4+iRYjiUsSUzWG4RiAWG0FkhSPUiUT4JnWTlDU1WX+BNCat4mUakkBepG4ukjYi34I0IJvzomu9Evd9nGS2qsxJ58yYPamXOaTw2pnkt5EGIx+GvQ457XW/xq9FzvTxoATrs9HvsgZBK1oi+EJUBNKObycaOa8NpoJETrEdI67a9OizuL5KMQrJfW4A5UIyPEk4Jj/HZpWI16NflT5RwuFG0a1W8tpeNn8J99mxlXoxCgSy2TlFd30btxoJ5WwgmN7qFi+d7USzVxl9JGFs8g+IHMbLswBVkUVtJ9KQ+Sl8nhIWdyseqf3kx9iMBjIbritysPuOayz6vk+Ph3mRDbpV0IF8Taklru5/arsbsv931t0aT7aciES2HcPBgyn94L1/y37fO+feA/wU8Gbv/Y5cj3iMq/H5unPudwjZ6Xf3Z9smFUNcIVWMlaCSwJNgG8Vt5Nx9zLzuEAT0FZID20ZYaX8r6BXJBckfoZX4cjzWRVJJlefia63yJWQV7dUxx7aCV+dQYIDdX2ajRVKosY67bI4jbUpCewAM+qR8zlvpufVw0oysU13ahzVVPUGoQfYYidhs7oiNOrM+G0ji2vqcKg1kOb5WiZbd0ch5bTA1JPIILwDQXQ9LAOVlbM0R+3nUhXHwa3SrRlU2Ce4hXqzySIBk3ipLNoog5HsxfLgVzTQqR7K4vspmJ9hcT69tsTF3itkihPW2hzAahfwRCNpIMYL+ghIGUxl1EdZGjJeCesVdmYqKuNaHoLWIGMbNP+OrexEHpPLvs6TOkDur/g5rWtw4bLSa7eVugwP0bKPi7Njs+yISou79ROYs7yhH92cXngTOubcRHI5/0Xu/scc2c8Ap730/vv4J4Ffv+mQlSTBboSsfiExBsu1LyCunApJDHoITXFqJJQ5FXyk6y2oH8r08RxKgy2Z7nUP5GtbprUgtaUcSuvZ8VhjrWnQerejlj7hGItSWeSjwoOanAVpnCJ0N1ZL2LODqRCNSu0hdUxMh6VpbZswiGGvuk3Z4fuw42qYimS1SHa8XuVPtrEbOa4Op8InMsBW0i/U12kPorAbz0VYbbs/PsLi6FrsLpsZQQCXk5NAOPpA12sO06p9lSLsMAndUpLwICTU5nXurA7rlBmXrFLODrWiuCmG9M8NAam41aEdli9jB7xS352eCScus+FWI0BYkHNcQIJiApLmE94ogm62NX/vZay6jyaqeBR+oqc1mRa72XCkHJI13vIYYUOWdSOOw5FZS0KfHbbqV9jNeFl4EIg1m4sJyHspRMdHjPvGvCHGizzrnvu6c+7cAzrlHnHOfids8DHzZOfeHhFyQT3vvP3vXZ9okNYayhCHb/BVSS9sb5ruR2f4CKeppfNVuI7isCcduc416/oYEuIhC/hFLBBK4+l7fWd+ByEVKpkxEkBz8ei3/zkWz/Zp5tuShY1bH7pKc2T8EvpFISfcR6hFYCmXWNYtE5IfRPbZ+kBapsZVCqK1/ZQSh1/sMKULrInfqbNjIeW0wFZqIY5vesE/nB7D1Smh1YFSE/IyVhTbDhXYVTmqFks0LSbkMqecHUDVKSmaVkt76GjfnztIr+7xULAah3IJ+0eN0EfwyZavg9NoWwzZE7qHVCuY12iFqqxiVFTFJsEt4Wg1Dwlrl020XRZGDvuvTY54+a/SiHyP8RNY0pRBo6xTvxkrEtp7WeOvbgrLKwUn73a4dQ/4PmRZtSLI1aZ1mZ96KhTXf6dx3gveO0dbBr9i894/t8flV4O3x9R8Bf/qBnfQywX8hJ/UNUmVeObqVhGdX5dr2PHUHunUIt0i+BpnALpMivOZJob3L5lgaC9Qd9ZIKMulIoMucYxMJ10iahTQNEZnI4kp8foygSSmwYDy/RK9loqs0EpmLdKDX1E1q2m6Zeqn7K9RzZ3R/F0mVABS9prGskMjwBmP+EEgmtRlzgfsXYWz0vGZKNBEHlXN6s3MKNwqr/a05qnIkciT36VW5DAUlKyzWHOn9ohcEPCXdcoONdrcSfN1hiGjamOsEP0AROiS2h5vcnDsbBWyrVnF2szOTNJEBDBagGMXeIqPtSuuRqciakOTohp3lQCAJdznah9GsJlMYpKQ9XeNtuiyyUkVC6XzjIbRJKxihVrba3lb/1bY7TWxlpdXsZl6z761mo+NpfHcXvujYLlsTPY4NZkh5B1oxK4xUEU0KO1VSoEjiEklYd0glPSD5WWzoLiTBuEQyP12iruVA8iWIFFbi8a1jWUK0Y/aTsFUeijSacZeXPcZjJMFtfQ8aqzSfJeDGzXoLXiCYsHoEE5JKwpux2/Bh63K2Wp38QNbJL9+NtBRNq3FnvHW0c47UrvdMHNcdNJEmzmuDqSCRYuRpD2HrTDAfMQwO7Jlh6OEhM49MPyUt0xWwHkI7ZLbqiT4qCtplPYmvuz5AUVRyqK+0FytNJYwnOMxvz88wO9iqeoT7TvCLlK1AepudGWYH2yEMmKJKYpSpyY5JpKCVvcxIKlVvS4uMVyO2UKkW2+N83Fmerk8lXtRpcbOmrUhbUOCBPlP5eFs3C4gmrHqP9nAXw9n69FBHRhHocFJTFsRKdcVkj+MCRxLYNslQwlomFglam58hR7G0g/MkgSfyEMFYB7dW5+cJTmWrseg7EcjAbGt9M3ovjWWZtGK3kUvSoiCRpNbDl0mahR3fCrByMxFfy3zeOVcPa642UBfBhwk5Iz4RjbS3i9Q1JOszsj4n7bNkruMCibwwx9E9f4K6D4czpNa4d6id1cR5bTAVJOJPhVpVt+dn6C/MsHUmkIfqVcmMAkFgLvJS8KGU/UqQbdClV8YaUeWQ9jCs7k+vbSW/wmibYgTtchgq9JpVN1C1vxVmB1t01mOP8DKM0a0HTaRsneL02hb9hQ6352folhtVUIAlpNko7GXusR0BbQSXMtBtKO0Ki5XPQ50KlXtixy0SGDejjZ8PiP6StYpQ7Lbj4bzSWPbyaYicbb2zeofEcA/sMffFtoNBa7LHccEsiQAukla5MsVAnQhkbrpGWq0rOss60a3WMJ4PcoV62RKt/qFeRVg5FtJ8bMLfMkl7uka937r1vyjayUaCjcx2EtDSvvT5hXNJA9J1iNSEFoSckduEGphKOGwBrp5RLpJYpl56pWMe+g0gEZ71H+l+Wa1ExGRJpfrLzBDIbYl90cR5bTAVJOJiTcPe6lZoAEVwYA+j7EpdDzersNMutyt/hOpi9YsevWFwKBejEDY8sx6EYXu4abSIWNl3fbUye0GIArNE0h7G/JARId4/akvq0rfZOVWF+W4UXfpz85VvxGojAIvDlUq4yvQkTWJc8FsHtYhGVXiBWk5J2L5VaUEiE2kZNopqVAn35LsZT460xywpKu1D+6uEZCh82TU+m6jFUc8hsUUoJ8Jowsdxgkw+0jxsCKlMOVo5S9DaSKt56hFNUC9oKA1HxNGK+4kYoK4xSJvQs7SFJZJQtWYhrdq1bRXqGvFNkhDWMe35dd2WWCwB8vxObUbXwWlqLXFV9FBOcME6+yX4rU9D96Fltr1IivDSvnLWL5Hup7Q+kUmVxLgRD7J/2ROgmfM6YipIhBGsL5yivzBD2Sq4PT/DxtwpOregvxB+MQnm8eQ4RSbNMgwmlViivWyF3h+DBeiWGxSjbdwgFFUE6C/MUIxgcXWN2cEWvWGfdjmsCGZ2sM0w9hABYB0oArkVo23mVmOr3FFJ2Sp4+NYq7XIYmzFF3aksaZdDVlik367X5wr9UOZrK/TxLHDrEBdxghpqjbDRaiUtbNVcJTnasiyvKFeq49sESFUENiMHqIhJZj9bgysECGzU3m+azPREpq3a5/siNF5o1p9NK/GLpNWworK0irYmq3nzXsQwMMeQsLMObaj7MJbiZ4rmkoZhV+wStlDXMkak0vISrl+mHvWk+79CMPM8YY57gyDvn6AeTmsd5zuwlIjmhvl4EJMKOUcomKzmTw+nIALdv+eo+1lkvrOOcptYKPOcfEvS/GSi0z3QsZZI5FSR32mYn6HOZrugifPaYDpIpJVKrCu8tru+DQW84taA8eiic9ysBDYkM1SbTbrlBqfXtioh31mFMzdD0cTBQtAshu1ZZgdbzLwszSScf1i0a6Xe28M4Lk3Ul5M5a6udxlyMSjbmQk2tUREKM7bLIf2iRzEqK02pV/ZrYb9Ws2jHsGNFVUnwVtn58Rq75caO1T8QS7DUkxJtSRarYYkYFDUmApADXaVnrEPckraelZioc1r/jYhdn+9Wbn8HmvhnG+dOmW2kLVjzjcwxkEwwMrXYlbCEvgS2hJzCUiVgrclMK3xr3tJ3i8BgI0VWibikIS2RVu/L1IMAniMRj/w1uk4RqE3wU6SZxtqBKrpptLULycwA3yHVzDJ+iHE/iHwhctyrBMpS3EVjEZlarQ4SIS8TeEHXps+sBjnwgIM1T2KcPdDEeW0wHSQC0UkdHFRlK0RoDRbgxtn5mqO4EnijZOMvRiXq0zEs2mx2TuEjOfgODOaif+N68F111weVI18oRnD+1lrVkGpmPWzbWScWTAbaYZ9WGR5yrA/bs1VEV0krttWNmeDtXlXsUd39JJhtL/gNUzxS/VEg1b9SsIBMbtJyZHZaaS9WJjSVj9F4uuUG/aLHRtGtBH67HFZaiPWFKOlSjzDO2ZqvxvZqUbn9lMGeiMe28J0Iaqk9yeM4weZ9yHxihZuIRDZ7a6O34bwXSALYOtplPrJJdI9RLwcvgbhCIpgbZnu6aUGt1bjKhIxrSTYPQ0LcRp7J8a8wXCVOjkhRY5bIWjFZb3HGhNICzMT78GME8rhNIJOYq6EghCdIOSHyOeladb8U1nye+j2Uz2lEisYibqvoOGkei2a/qhjjdSbJE2nkvI6YChLZLuDM97foxAXJ3GpoR9u5Cr31tco2D6kZU3d9m97qFt1yo4rggrBin7u5jRvBzA9CgmAxChoJo6B5dK6Hz4WZ6CvZmDtVkQeEPBXfimQzBFox1HcUtlEEFwSSkmBvlcHE1Rv2K6e4VvyCmmGNl4bfLftbvcvbw83K5NZb3drhwwFin5IgvEMS5ChWHd6sfD9AFblm0aooKVU/FmyIsqLDUhXlFBo8XsFY1zSRX8TH+zzJ47igRVioylyi56+ThJrMVBKqEuzjznBpGBD8ENafIoK6QtJUZKISYUnLkJlHYcMSkjJZifREKhKk1iR1mWQqWjL7SDDfoN50S6YkaTgS+BoDW8bP8nwkD2+c3KeBH4k361b4zkaV6XoskVg/je6rDWEW2UjLEMkp10YkpONZsyDEsUzQ2bCJ89pgKkjEbQPrQVi3442cWQfm0jZapRcEQb0xdwq3HkxVM0O1fw0C8+WHZ1JD07lYQiXKtU5sMMUIZm5RNZ6C4AehFcmjE/JURkXc/0wYoyLGZl6O4cgEzUal44dFm9nBFi8V8oMEgV75M8p+5QAPWkAS3UKLerdACEJ8o92lGJUM27NVcyxhvLHVJm3KVkF7uEl3fRDMXdH3Mx5JpvEoys3WKLM+j5A3krSfgrLKjLflXjap91IJzxPo6k1U+7dJWoSEljUZSbhJAEKKvLphvpcm8DgpL0IraGuikjCXaUqajchBEVlLJOK4SDDPiGDkL4AkOCWc5QdZIiXkSQhfIRHRMvWoM/1m8tFI+5KJ64LKuwM8CoObgIvXp3Lw8QZ1zoXvBuY610jmQYUqW+KQNmLzVKzjXP4cEc4NgvYGKcdR2khNtQDif2BPNHFeG0xMIs65wjn3v51zn4rvzzrnnnXOfS8+v8Js+5Rz7rJz7rvOuZ+847FjmTCt8lkPpizfgY25TvAxRANLd7iRfBdm0d6nVwnM02tbsECwR9sfpiSwvSZwO5BCexgerTJtNyoCUcy8HI8xBM7G8UXTlhulUvGK2OoON9jszJDKvo+isz0I5VFRVBFTm8xGn0I9gklC3SYjCsVoO+XLtE7F7VM+iPwZ3TLcp7JVULZOVc/D9mz04YRlWm99rfKH2GONj8eWSEnFI0eVNqIe94oksw3FbMjzvjiCP9tBzuvqmiREZWKRXV1Cec18L0Enk4sihUQ+IgYbiWVX2tb0Mk99xa2foEMqya4xLbm03XnzucbSIWlOkISzxiVtSdrHBZJQt2QibUPjEfHpZskMJ9/HAILpaIMgrE8H/43ugXwh0ujkJxEZijghkbINib5C0qYgJUTqsxVSdWD5hHDQUhZpH6JJd09kEqnwSwQPl/AB4Ive+9cAX4zvcc69FngX8DrgbcBvxNa6e2MbWKVK6rv1aIf+3Dw3zs7HUuvBoa7w3d76WtAo1qlu/PnhzUpQQvp8ay5uJ81GkyjCjUIosVtPIcVeLXELAhGJfIjHMSurmfUUrRUIpWBYtEPOyjAI02FhHebtKsEQdjqsIYXgWp9Ft0zZ9ovD8K8YtpOfQv6QkhatsqzMVfIzCd31AcP2bLWNsveVnZ7MT0WlDWlMm7QZ0qZPj5dYrCKvRByh6GVZM8ulcikTJFIdzZ/t4OY1hGuykU1PkBzh1vk9Imkn1tku4S1hqbmrEFURyxJpdd0x2ymSS76NJXbmPIh4lOcBSQhDvVyKMs8V8aSHHOZXzDE0finMOo/MQ1Ywj4DzLpjpADoujaND3KAH3IZWNxGYjWzTdeq8CjLQuXRvpGmsmOPrfjxG8IXo/up6L7EzaZMrVGXp90MmEXDOXQT+CvBR8/E7gI/F1x8Dftp8/nHv/dB7/zzhJ3jDHUfxSPQ9lNBbHdAuh/TW1+iubwf/wuog1qtK0VPMAbeCOWrYnuXM1RDVNfNy2mRmSCAowvHXHzGXHMN9O7eAdvCbbLVjfaySMFmGpAibVRIRrYdtWQ9ajExip9e26A43ajW15HtQ4qDtNwKpp0c9CbFVM4EVoyCcFZG2UXSrHJV2Oaz5QICKfENByXDPylZBf26+KtWiaDBIVZEtLKEM42tbGdlqMKDQ43oY8MTFF+HQ/2wHPq8hLEQeox41tUaKcrJ5CuMCTaYmCWet/q2JSIL5IsHUpe+0nYS7nq3/xGo+yhNR1JcITStxkZnVpERggs1bkS9jxby2TmoR1pU4Dpn7bJguwChGYVXhvb30nQhShCE/BuwsrCjigkTYulcKm5Y56wqpbLzO9U0z5sr0uESwc99By84kAsC/BP4BQWcQHvbevwAQnx+Kn78K+IHZbpnUEKCCc+596t71xzfAL8QGUEVqBLUx1+Gls50QjjsMyYid68GPsdUmCPj16IBfHcAo5nVIg5Cz6hHgTNAa5lajGawFWw+ZAUX1Wv6Q6rOYZMiASlMycpPBI8H01VkN497sBKGtFT9QaQW2x4ZKocxS758evg85IEryU1SXOj3KaV8NvSgqZ76iwxQtNmzPsrIwz5mbW9U4RMQhXLgwBDaqxmfJIJisRpWWIdhijrfp1nrVW/+MNKuJcLh/tgc+r6E+t3npj5NT2TrI5d+wIbHXSKtdqJtbZCqSwJPG8jiphMmKOYeOrePb/BCrzYiIxs1MHeDPx+3lgxFByFEOiXSUpChntwhPJClIqIsoLo4dQ0K8+o1nSL3Mu+Fhhf8lgoCXb8VqHtJURGAiNTueJXP9lmieM+N5kpQLI/9KNb4z3DFjHRpNInc0VDvnfgq47r3/mnPuL01wTLfLZzuao3jvnwaeBvjR1wevSOd60BTmrm/TaW/xsi3Tvx5MT1sPBYd4S3+QqAHMRAGv2ltAmBBzBCJ4mco05eeC+WrmViAvt0rNie9W4+c3qVZ8v/LkU/tf9YJ5vdvi+37L4thjFuZ4Y89/h3/HaUbMDraCP2l9QNk6xcvngsqtJlqzgy36c2FJGRp6raFqAJD6ktQLShY1clH/FGkq7Wj+qn+uJMUJfCJyQh8CDmpeQ31uu9ddCtvIHPR16n4Am9sgpzmkirP2ttnIoJH5TgJS5KQwXmkgIh35Qi4QBO8o7PvU1j+648VPA/7xo7+aiHaJ5Px/nJ1mOWk2i1A1wdJ9sJqHzIzyg4hAHyOVz79gtpX2BLDoIvnuNjUMDnFeHwUm0UTeCPw159wV4OPAm5xzvwW86Jx7JUB8VjzUMvBqs/9F4Oq+Z3BRkLdg7uY2tEMG+5mboQyKzE1cT2ThBlRmKlpRM1kN5FCZoDpxm3V2ZtvG3A8XTVq8HCOvohnLiYAmkH3ThFk2aQ832ezMVGar7vp2pa0Im52Zytlva3uN+y6GsX6Wjr0bSkKnSauZ2NyS8bIue+Jw1f6Dn9cQ5IuEk8whcjrLX3CBVBdL801ahQSe9aso0kifSUBZzQP29nXI+TzP8Zrful6bCyK/ig2PthFYl6nn6egeS6PQfleo3z/d0w7wFeqlXebNMVkmmdr2wEk3Z3nvn/LeX/TeLxEci//Ne/+3gE8C74mbvQf4RHz9SeBdzrm2c+5RQr2C39/vHG4IzMHgbHCEry+cCi1q24lUfAc4F30PcwTDgtyhIoJzBNIoSAmCxM80OYbgrpMc5qRnmzvCOmFurJLI6hig6pIYfTLtYTDxyZGu6KxREU1fUVsYMstajDIREbRiSK80iXEfiEhD2skKr2BI6qyoPi8TOdXhUP9shzGvgbQalilKZhwRh3weS6SchGVSlJD8BEskYhFsdJPI4Rr1yC09r4y9vhLPI+f5cYAV7DKXyZQn85bVyrStnOTSOEbUizXKZ6JotQH15M3zBNOezFvSVtYgmLMmqeI74eMY4n7WIb8GPOOc+wXg+8A7Abz333LOPQN8m3Bb3u+937cj0dXTF/iVh95b/1DmG2smUsWDOYKf48fGDrQw9v7VnDjIad9bDf+4/kIwY4k02uWQYXu2ChNuUVY9WmSCslCLX4AQUVYnEpmvRhRVY6z03agyc403qdoV+rMdLR7YvAa4sHiV977xeJiLph4KtZWZyfqURiTH+Lx5fzFuLzOhCKZDvZ+8fE4r1LVB+bOkiUhjrAIKVGV4H0zHvD4w3BWJeO9/D/i9+Pom8OY9tvsI8JH7HFvGPWCW4DhXTbDQ5jcqnEUgk1Afa7MyYamIpTQN1cQKGkrSTmzxSxuBNYomMZGNukza7/cyhe3AEfzZ8rw+JpDpSUQgrUumOpGDQqdtwIA0kAsELW+R1GtFznybdQ8pWOAydfOhtJAOMOhyxxBfja2hmIqM9YwHC2XPK8qtGG1XmohySqQ9KEPd1s9qx8gEhfgqNFnZ7bb2F6SeJYKaUtXbGE+PT8Q59yHn3A9jH+qvO+fevsd2b4uJhZedcx+4v7Nm3DcGAD4FEEBdcxBkclI0mc2xkTNe5KPKADZjXfk2yoK3UJ6Ijjlp2ZMGz2vn/a4BJocK51wf+O5Rj+M+oViY44rDGP+f9PuXt4EAAAg1SURBVN7/ib2+dBcvef7uVyc70gfc17z3l+5lEM65DwFr3vt/vs82BfB/gbcSDCJ/ALzbe//tuzzXcZ/bx31ewxHP7SbOa4tpic347r3euGmBc+6rx/kapmL8nloOzhHjDcBl7/0fATjnPk5IOLzbP9uxnttTMS/uE0d+Dc2c1xWyOStjujC52n9eCX3x8b67PNMvOue+4Zz7TVsfy2Di5MKMjDuiwfN6WjSRjIy7jWK5sd/q0jn3BepNXoUPAv8G+HA844eBXwd+fvwQe4wwI+Pu0PB5PS0k8vRRD+AB4Lhfw9GP/wGGQnrv3zLJds65fw98apev7i25cCeO/r7eH477+OGor6GZ87rCVJBILBNxrHHcr2Eqxn9I5SGcc69UfSzgZ0i1Yy3+AHhNTCz8ISEh8W/e7bmm4r7eB477+GEKrqGB89piKkgkI6PC4cTT/1Pn3JOENeIV4G8DOOceAT7qvX+7937knPtF4HOEGgi/6b3/1qGMLqN5aPC8ziSSMT04pMxe7/3P7fH5VeDt5v1ngM8c/IgyGo2Gz+sjj846DgldzrlXO+d+1zn3Hefct5xzvxQ/f7Bd8A4YB97F735hu47e6XEMkOf24WGq53bD5vU4jpREYuLLvwb+MvBa4N0udJCbNoyAv++9/xHgzwLvj+N8sF3wDh4H28XvfqF4+kkeU448tw8d0zu3GzSvd8NRayJV4ov3fpNQkvsdRzymHfDev+C9/1/xdZ8wWV/Fg+6Cd4Bwh9HF737RrGqneW4fEqZ+bjdrXu/AUZPIsUvocs4tAa8H/icPoAveIeJAuvg9UHhCQdRJHtOPaZwD+yLP7QNCs+b1Dhw1iRyrhC7n3Dzwn4G/573frxPNVF2XM138Jt1ll88OfvzNUvunag7cCXluHyCaNa934Kijsx544stBwTk3Q/iT/bb3/r/Ej19UbLZ7EF3wDg7q4vd2Qp3SM8508Zua8Ter78K0zYE9kef2AY+/WfN6B45aE6kSX5xzswSH1yePeEw74JxzwH8AvuO9/xfmqwfbBe+AcGhd/O57oDTJdpzn9iHgWMztZs3rHThSTeQYJXS9Efg54P84574eP/uHPOAueEeA6Rq/QiEbgDy3jxzTM/4GzevdMBX9RDIyANyZS55LE/Zd+N1777uQkXGYaPq8PmqfSEZGwiHVGMrIOFQ0fF5nEsmYHjRc7c84oWj4vM4kkjE9mK4OcBkZDwYNn9eZRDKmC8c0QiUjY180eF5nEsmYHjQ8nj7jhKLh8zqTSMb0oOEOyIwTiobP60wiGdODhq/YMk4oGj6vM4lkTBca/GfLOMFo8LzOJJIxPWh4KGTGCUXD53UmkYzpwSGFQjrn/hPwp+LbRWDFe//kLttdAfpxVKPjlkmcMSVo+LzOJJIxPTi8XtR/Q6+dc78OrO6z+Y97728c/KgyGouGz+tMIhnTg20OtTFPrGD714E3Hd5ZM04cGj6vj7oUfEZGHYfbvOcvAC9677+3x/ce+Lxz7mvOufc9sLNmnDw0eF5nTSRjujB5UenzzjlbGvVp7/3TeuOc+wJwYZf9Pui9V2+JdwP/cZ9zvNF7f9U59xDwrHPuOe/9lyYeYUaG0OB5nUkk47jixn4OQe/9W/bb2TnXAn4W+NF9jnE1Pl93zv0O8AYgk0jGQeLYzetszso4qXgL8Jz3fnm3L51zc865nl4DPwF88xDHl5FxLzj0eZ1JJOOk4l2MqfzOuUecc5+Jbx8Gvuyc+0NC+9RPe+8/e8hjzMi4Wxz6vM6dDTOmBs79GQ//fcKtu8euA1zGyUTT53X2iWRMERqe2ptxQtHseZ1JJGOK0PBKdRknFM2e15lEMqYIzV6xZZxUNHteZxLJmCI0+8+WcVLR7HmdSSRjiuA51PoQGRmHgmbP60wiGVOEZtuOM04qmj2vM4lkTBGarfZnnFQ0e15nEsmYIjR7xZZxUtHseZ1JJGOK0OwVW8ZJRbPndSaRjClCs1dsGScVzZ7XmUQypgiH3L0nI+NQ0Ox5nUkkY4rQbLU/46Si2fM6k0jGlKG5an/GSUZz53UmkYwpQrNXbBknFc2e15lEMqYIzf6zZZxUNHteZxLJmCI0O4ol46Si2fM6k0jGFKHZUSwZJxXNnteZRDKmCM1W+zNOKpo9rzOJZEwRmq32Z5xUNHtenzrqAWRkJGjFNsnj3uGce6dz7lvOuW3n3KWx755yzl12zn3XOfeTe+x/1jn3rHPue/H5Ffc1oIyGo9nzOpNIxhRBK7ZJHveFbwI/C3zJfuicey3wLuB1wNuA33DOFbvs/wHgi9771wBfjO8zMvZAs+d1JpGMKYIckJM87h3e++9477+7y1fvAD7uvR96758HLgNv2GO7j8XXHwN++r4GlNFwNHteZ59IxhThhc/Bh85PuHHHOfdV8/5p7/3T9zmAVwFfMe+X42fjeNh7/wKA9/4F59xD93nejEaj2fM6k0jG1MB7/7YHdSzn3BeAC7t89UHv/Sf22m23YT2oMWWcTDR9XmcSyWgkvPdvuYfdloFXm/cXgau7bPeic+6VcbX2SuD6vYwxI+NuMY3zOvtEMjISPgm8yznXds49CrwG+P09tntPfP0eYK8VYEbGNOBA53UmkYwTB+fczzjnloE/B3zaOfc5AO/9t4BngG8DnwXe770v4z4fNWGTvwa81Tn3PeCt8X1GxpHiqOa18z6bfDMyMjIy7g1ZE8nIyMjIuGdkEsnIyMjIuGdkEsnIyMjIuGdkEsnIyMjIuGdkEsnIyMjIuGdkEsnIyMjIuGdkEsnIyMjIuGf8f4LYisAKKJJIAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.subplot(1,2,1)\n", + "plt.imshow(n['U'], vmin=-10, vmax=10, cmap='jet')\n", + "plt.colorbar()\n", + "plt.subplot(1,2,2)\n", + "plt.imshow(n['V'], vmin=-10, vmax=10, cmap='jet')\n", + "plt.colorbar()\n", + "plt.show()" ] }, { "cell_type": "code", - "execution_count": 23, + "execution_count": 31, "metadata": {}, "outputs": [ { @@ -58,19 +96,20 @@ } ], "source": [ - "n = Nansat(fn)" + "fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc'\n", + "w = Nansat(fn)" ] }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 35, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAARcAAAD8CAYAAABdJ+AhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAANfElEQVR4nO3df6jdd33H8efLtFaZDtv1toQkXTLIH0tlVrlkQsdw6tZMxfSfQgRH/ijknw6UDSSZsOEfBbc/xL/6R1BZwB8hoNJQZFuIigyG6Y1tZ9MYe7WuvSQ014mo+6Ou9b0/7ld2TG56T5L7zv2e+nzA5XzP53zPOe9L2me/50e/SVUhSevtdRs9gKTXJuMiqYVxkdTCuEhqYVwktTAuklq0xSXJniTnkiwmOdj1PJLGKR3fc0myCfg+8OfAEvA48KGqembdn0zSKHUduewGFqvqh1X1S+AosLfpuSSN0E1Nj7sFeGHi+hLwx1fa+fbbb6/t27c3jSKpy+nTp39cVXOr3dYVl6yy9huvv5IcAA4A3HXXXSwsLDSNIqlLkv+60m1dL4uWgG0T17cC5yd3qKrDVTVfVfNzc6uGT9IM64rL48DOJDuSvB7YBxxvei5JI9TysqiqXk7y18C/ApuAz1XVmY7nkjROXe+5UFVfA77W9fiSxs1v6EpqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIarFmXJJ8LsnFJE9PrN2W5ESSZ4fLWyduO5RkMcm5JPd1DS5p3KY5cvlnYM8laweBk1W1Ezg5XCfJLmAfcPdwn0eSbFq3aSXNjDXjUlXfAn5yyfJe4MiwfQS4f2L9aFW9VFXPAYvA7nWaVdIMudb3XO6sqgsAw+Udw/oW4IWJ/ZaGtcskOZBkIcnC8vLyNY4haazW+w3drLJWq+1YVYerar6q5ufm5tZ5DEkb7Vrj8mKSzQDD5cVhfQnYNrHfVuD8tY8naVZda1yOA/uH7f3AoxPr+5LckmQHsBM4dX0jSppFN621Q5IvAe8Cbk+yBPwD8EngWJIHgeeBBwCq6kySY8AzwMvAQ1X1StPskkZszbhU1YeucNN7rrD/w8DD1zOUpNnnN3QltTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIarFmXJJsS/KNJGeTnEnykWH9tiQnkjw7XN46cZ9DSRaTnEtyX+cvIGmcpjlyeRn426r6Q+CdwENJdgEHgZNVtRM4OVxnuG0fcDewB3gkyaaO4SWN15pxqaoLVfWdYfvnwFlgC7AXODLsdgS4f9jeCxytqpeq6jlgEdi93oNLGreres8lyXbg7cC3gTur6gKsBAi4Y9htC/DCxN2WhrVLH+tAkoUkC8vLy1c/uaRRmzouSd4EfBn4aFX97NV2XWWtLluoOlxV81U1Pzc3N+0YkmbEVHFJcjMrYflCVX1lWH4xyebh9s3AxWF9Cdg2cfetwPn1GVfSrJjm06IAnwXOVtWnJm46DuwftvcDj06s70tyS5IdwE7g1PqNLGkW3DTFPvcCfwV8N8mTw9rfAZ8EjiV5EHgeeACgqs4kOQY8w8onTQ9V1SvrPrmkUVszLlX176z+PgrAe65wn4eBh69jLkkzzm/oSmphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdRizbgkeUOSU0meSnImySeG9duSnEjy7HB568R9DiVZTHIuyX2dv4CkcZrmyOUl4N1V9TbgHmBPkncCB4GTVbUTODlcJ8kuYB9wN7AHeCTJpo7hJY3XmnGpFb8Yrt48/BSwFzgyrB8B7h+29wJHq+qlqnoOWAR2r+vUkkZvqvdckmxK8iRwEThRVd8G7qyqCwDD5R3D7luAFybuvjSsSfotMlVcquqVqroH2ArsTvLWV9k9qz3EZTslB5IsJFlYXl6eblpJM+OqPi2qqp8C32TlvZQXk2wGGC4vDrstAdsm7rYVOL/KYx2uqvmqmp+bm7uG0SWN2TSfFs0lecuw/UbgvcD3gOPA/mG3/cCjw/ZxYF+SW5LsAHYCp9Z7cEnjdtMU+2wGjgyf+LwOOFZVjyX5D+BYkgeB54EHAKrqTJJjwDPAy8BDVfVKz/iSxipVl70dcsPNz8/XwsLCRo8h6SolOV1V86vd5jd0JbUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1mDouSTYleSLJY8P125KcSPLscHnrxL6HkiwmOZfkvo7BJY3b1Ry5fAQ4O3H9IHCyqnYCJ4frJNkF7APuBvYAjyTZtD7jSpoVU8UlyVbg/cBnJpb3AkeG7SPA/RPrR6vqpap6DlgEdq/PuJJmxbRHLp8GPgb8amLtzqq6ADBc3jGsbwFemNhvaVj7DUkOJFlIsrC8vHzVg0satzXjkuQDwMWqOj3lY2aVtbpsoepwVc1X1fzc3NyUDy1pVtw0xT73Ah9M8j7gDcDvJvk88GKSzVV1Iclm4OKw/xKwbeL+W4Hz6zm0pPFb88ilqg5V1daq2s7KG7Vfr6oPA8eB/cNu+4FHh+3jwL4ktyTZAewETq375JJGbZojlyv5JHAsyYPA88ADAFV1Jskx4BngZeChqnrluieVNFNSddnbITfc/Px8LSwsbPQYkq5SktNVNb/abX5DV1IL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWU8UlyY+SfDfJk0kWhrXbkpxI8uxweevE/oeSLCY5l+S+ruEljdfVHLn8WVXdU1Xzw/WDwMmq2gmcHK6TZBewD7gb2AM8kmTTOs4saQZcz8uivcCRYfsIcP/E+tGqeqmqngMWgd3X8TySZtC0cSng35KcTnJgWLuzqi4ADJd3DOtbgBcm7rs0rEn6LXLTlPvdW1Xnk9wBnEjyvVfZN6us1WU7rUTqAMBdd9015RiSZsVURy5VdX64vAh8lZWXOS8m2QwwXF4cdl8Ctk3cfStwfpXHPFxV81U1Pzc3d+2/gaRRWjMuSX4nyZt/vQ38BfA0cBzYP+y2H3h02D4O7EtyS5IdwE7g1HoPLmncpnlZdCfw1SS/3v+LVfUvSR4HjiV5EHgeeACgqs4kOQY8A7wMPFRVr7RML2m01oxLVf0QeNsq6/8NvOcK93kYePi6p5M0s/yGrqQWxkVSC+MiqYVxkdTCuEhqkarLvjx744dIloH/AX680bOs4Xaccb3MwpyzMCNs7Jy/X1Wrfgt2FHEBSLIw8X9cj5Izrp9ZmHMWZoTxzunLIkktjIukFmOKy+GNHmAKzrh+ZmHOWZgRRjrnaN5zkfTaMqYjF0mvIRselyR7hhN5LyY5uMGzfC7JxSRPT6yN6kTkSbYl+UaSs0nOJPnI2OZM8oYkp5I8Ncz4ibHNOPG8m5I8keSxEc84myfIr6oN+wE2AT8A/gB4PfAUsGsD5/lT4B3A0xNr/wQcHLYPAv84bO8a5r0F2DH8HptuwIybgXcM228Gvj/MMpo5WTkb4ZuG7ZuBbwPvHNOME7P+DfBF4LEx/nkPz/0j4PZL1kY356U/G33kshtYrKofVtUvgaOsnOB7Q1TVt4CfXLI8qhORV9WFqvrOsP1z4Cwr5ygezZy14hfD1ZuHnxrTjABJtgLvBz4zsTyqGV/F6Ofc6LjMwsm8R3si8iTbgbezcmQwqjmHlxtPsnL60xNVNboZgU8DHwN+NbE2thlhRk+QP+0JurtMdTLvkdrQ2ZO8Cfgy8NGq+tlwpsBVd11lrX3OWjn74D1J3sLKmQzf+iq73/AZk3wAuFhVp5O8a5q7rLJ2o/681/0E+TfCRh+5THUy7w12XSci75DkZlbC8oWq+spY5wSoqp8C32TlL8gb04z3Ah9M8iNWXo6/O8nnRzYj0HOC/Btho+PyOLAzyY4kr2flb2o8vsEzXWpUJyLPyiHKZ4GzVfWpMc6ZZG44YiHJG4H3At8b04xVdaiqtlbVdlb+uft6VX14TDPCjJ8gfyPeRb7kXe/3sfKJxw+Aj2/wLF8CLgD/y8p/AR4Efo+Vv6722eHyton9Pz7MfQ74yxs045+wcpj7n8CTw8/7xjQn8EfAE8OMTwN/P6yPZsZL5n0X//9p0ahmZOWT1KeGnzO//ndkbHOu9uM3dCW12OiXRZJeo4yLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWvwfldf/t/Dmb/QAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAADxCAYAAAAKooq5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAaNUlEQVR4nO3dfaxcd33n8fcnzgMlTZdEzoNjJxCpFtRZbSm9MqBoJSgJOG42pg/sOiuoRdkNdJMVaFtRp5G6aBFS9qFsVUHI3qWobguklmgaL3UJxtsqLWqa2CFAHMcbN1BibNkyu02CQhKu/d0/5kwY38zce+7MmZnffM/nJR3dOQ9zzu84n5PvOWfO/EYRgZmZ2TDOmnYDzMxsdrmImJnZ0FxEzMxsaC4iZmY2NBcRMzMbmouImZkNzUXEUpL0GUknJD3aM+0iSXskPVH9vXDAezdJOiTpsKTtk2u12dJKzLWLiGX1B8CmRdO2A3sjYj2wtxo/g6RVwCeB64ENwE2SNoy3qWa1/QGF5dpFxFKKiPuB/7to8hZgR/V6B/DOPm/dCByOiCcj4kXg7up9ZlNXYq7PbmIlZk34SSmeq7nsMTgAPN8zaT4i5pd526URcQwgIo5JuqTPMmuBp3rGjwBvrNkss5fJnmsXESvGc8D7ay77EXg+IubG0Az1mea+gWxo2XPtImLFEGMP5HFJa6qztTXAiT7LHAGu6BlfBxwdb7Mss+y59mciVoyzgB+rOQxpF7Cter0NuLfPMg8B6yVdJelcYGv1PrOhZM+1i4gVQ8A5NYdl1yV9Hvhb4LWSjkh6H3AHcJ2kJ4DrqnEkXS5pN0BELAC3AvcBB4GdEXGgsZ201smea7kreCvFq6Wo+/D6v4P9Y7p3bNao7Ln2ZyJWjO4Zm1km2XPtImLFmMAHkGYTlz3XmffNZkz2MzZrp+y5dhGxYnSfYjHLJHuuXUSsGNnP2KydsufaRcSK4kBaRplznXnfbMZkP2OzdsqeaxcRK0b2p1isnbLnOvO+2YzJ/gGktVP2XLuIWDGyX/ZbO2XPtYuIFSP7Zb+1U/ZcZ943mzHZz9isnbLn2kXEipH9jM3aKXuuM++bzZjsZ2zWTtlz7SJixRC5n2KxdsqeaxcRK4aAc+omcmGcLTFrTvZcu4hYMSQ4O/HBZu2UPdcuIlYMCc5ZNe1WmDUre679G+tWjO4ZW51h6fXotZIe6RmekfShRcu8RdLTPcv89jj3zdqrqVx31lVetn0lYsWQ4JzzRl9PRBwCXt9Zp1YB3wXu6bPoX0fEDaNv0WywpnINZWbbRcTKMZ4H6t8G/H1E/EPjazarY3xfFCki276dZeXoHmx1hvq2Ap8fMO/Nkr4u6S8kXT1kq82WNp5cQyHZ9pWIlaV+IldL2tczPh8R870LSDoXuBG4rc/7HwZeHRHfl7QZ+DNg/cobbFZDg7mGsrLtImLlEFD/KZaTETG3zDLXAw9HxPHFMyLimZ7XuyXdKWl1RJys3QKzOprPNRSUbRcRK0fz945vYsDlvqTLgOMREZI20rm1+71Gt24G4/pMpJhsu4hYOQQ09BSLpFcC1wHv75n2AYCIuAv4ZeDXJC0APwC2RkQ0s3WzHg3mGsrLtnzcWCnmzlfs21BvWe1jf83LfrOpyp5rX4lYObL3mW3tlDzXiXfNZlLi7iGsxRLn2kXEypH8jM1aKnmuE++azZzkB5u1VPJcJ941mzkNP8ViVoTkuXYRsXIkP2Ozlkqe68S7ZjMn+cFmLZU814l3zWbOyrqHMJsNyXPtImLlSH7GZi2VPNeJd81mTvIPIK2lkufaRcTKkfyMzVoqea4T75rNnOQHm7VU8lwn3jWbSU6kZZQ414l3zWZO8qdYrKWS59pFxMqR/LLfWip5rhPvms2c5E+xWEslz7WLiJUj+RmbtVTyXCfeNZs5yQ82a6nkuT5rXCuWtEnSIUmHJW0f13Yske7BVmeYEufaVmwGcj2KsRQRSauATwLXAxuAmyTV/JVha7VVNYdlSPq2pG9KekTSvj7zJen3qmLwDUlvqLFO59qG01CuYTzZHsW4at9G4HBEPAkg6W5gC/DYmLZnGTR/2f/WiDg5YN71wPpqeCPwqervUpxrW7nx3M5qOttDG1cRWQs81TN+hEU7Ielm4GaA888//2df97rXjakpVor9+/efjIiLBy4w2adYtgB/GBEBPCDpVZLWRMSxJd6zbK7B2W6jJbM9+aezhsn20MZVRNRnWpwxEjEPzAPMzc3Fvn0vuyqzZCT9w9ILsJJErl50KT9fZaorgC9LCuB/LJoH/QvCWmCpA23ZXIOz3UZLZrvZXMN4sj20cRWRI8AVPePrgKNj2pZlsbKD7WREzC0x/5qIOCrpEmCPpMcj4v5FW1vsZQVhEefaVq7ZXMN4sj20cT2d9RCwXtJVks4FtgK7xrQty6LbPUQDH0BGxNHq7wngHjqfZ/QapiA417ZyDeYaxpbtoY2liETEAnArcB9wENgZEQfGsS1LpKFHISWdL+mC7mvg7cCjixbbBfxK9STLm4Cnl7tn7FzbUBp8xHdc2R7F2J5MjojdwO5xrd8SEvCKRtZ0KXCPJOhk/HMR8SVJHwCIiLvoZHMzcBh4DnhvnRU717ZizeUaxpjtYc3o11sspYZ6O60ewf3pPtPv6nkdwC2jb81sGQ324ltitl1ErBzJu4ewlkqe68S7ZjPJibSMEuc68a7ZzEn+4z3WUslz7SJi5Uh+2W8tlTzXiXfNZk7yH++xlkqeaxcRK0fyMzZrqeS5TrxrNnOSH2zWUslznXjXbOYkP9ispZLnOvGu2UxK/BSLtVjiXLuIWDmSn7FZSyXPdeJds5mT/CkWa6nkuXYRsXIkP2Ozlkqe68S7ZjMn+cFmLZU814l3zWZO8oPNWip5rhPvms2iSPwUi7VX5ly7iFgx4ix4sbkf7zErQvZcu4hYMUKwsKruLzafHmtbzJqSPdcuIlaMkDh1dt1IvjjWtpg1JXuu65ZHs4k4tWpVrWEpkq6Q9JeSDko6IOmDfZZ5i6SnJT1SDb89tp2y1msi11Bmtn0lYsUIxKlm+odYAH49Ih6WdAGwX9KeiHhs0XJ/HRE3NLFBs0EazDUUmG0XEStGIBYaONgi4hhwrHr9rKSDwFpg8YFmNnZN5RrKzLaLiBUjEC/W7x9itaR9PePzETG/eCFJrwF+Bvi7Put4s6SvA0eB34iIAytrsdnyxpFrKCfbLiJWjBVe9p+MiLmlFpD048AXgA9FxDOLZj8MvDoivi9pM/BnwPqVttlsOU3nGsrKtj9Yt6KcYlWtYTmSzqFzkH02Iv508fyIeCYivl+93g2cI2l10/tjBs3lGsrLtq9ErBhN3TuWJOD3gYMR8fEBy1wGHI+IkLSRzgnV90beuNkiTX4mUmK2XUSsGJ3L/kYieQ3wHuCbkh6ppv0WcCVARNwF/DLwa5IWgB8AWyMimti4Wa8Gcw0FZttFxIrR+QDy3NHXE/E3dLq9W2qZTwCfGHljZstoKtdQZrZdRKwYAY1d9puVInuuXUSsII1e9psVIneu8+6ZzZyGv9lrVoTsuXYRsaJkPtisvTLn2kXEipH9jM3aKXuul/2y4aBeIyVdJGmPpCeqvxf2vOc2SYclHZL0jnHugOURiBc4r9bQBGfbJmHSuZ60Ot9Y7/Ya+VPAm4BbJG0AtgN7I2I9sLcap5q3Fbga2ATcKSlvGbbGdM/Ymvpmbw3Oto3dFHI9UcsWkYg4FhEPV6+fBbq9Rm4BdlSL7QDeWb3eAtwdES9ExLeAw8DGphtu+Uz6YHO2bRKyF5EVfSayqNfIS6tuiYmIY5IuqRZbCzzQ87Yj1bTF67oZuBngyiuvXGm7LalpPU/vbNs4Zf6eSO0OGJfpNfKMRftMe9lX7iNiPiLmImLu4osvrtsMS6zbPUSdoUnOto3TtHI9KbVaPaDXyOOS1lRnamuAE9X0I8AVPW9fR6dPe7MlTeMpFmfbxs1PZw3uNXIXsK16vQ24t2f6VknnSbqKTj/2DzbXZMuq8xTLubWGJjjbNgmTzvWk1bkSGdRr5B3ATknvA74DvAsgIg5I2knn5xoXgFsi4lTjLbd0Gu7ttA5n28ZuCrmeqGX3bJleI9824D0fAz42QruspSZ52e9s26Rkvp2VtzzazMl+79jaKXuuXUSsGNkPNmun7Ll2EbFidLuHMMske65dRKwY2c/YrJ2y57r2lw3NJqGp7iEkbao6STwsaXuf+ZL0e9X8b0h6w1h2yIzmcg3lZdtXIlaMQI10D1F1ivhJ4Do6XxB8SNKuiHisZ7Hr6XzPYz3wRuBT1V+zRjWVaygz274SsWI02D3ERuBwRDwZES8Cd9PpPLHXFuAPo+MB4FXVt9PNGtVwtyfFZdtXIlaUFdw7Xi1pX8/4fETMV6/XAk/1zDvCy8/E+i2zFjhWv7Vm9TSUaygw2y4iVoxAvFi/64eTETE3YF6djhJrdaZoNqoGcw0FZttFxIrR4L3jOh0lujNFm4gmPxOhwGz7MxErRoP3jh8C1ku6StK5dH6NcNeiZXYBv1I9yfIm4Onub4iYNanhz0SKy7avRKwoTTxPHxELkm4F7gNWAZ+pOk/8QDX/LmA3sJnOrxM+B7x35A2bDdDgr3EWl20XEStGk1/KiojddA6m3ml39bwO4JZGNma2hKa/bFhatl1ErBgN3zs2K0L2XLuIWDE6T7Hk7WPI2il7rl1ErBjZ+xiydsqeaxcRK0rmg83aK3OuXUSsGNnvHVs7Zc+1i4gVI/tvUVs7Zc913j2zmbPC7iHMZkL2XLuIWDGyX/ZbO2XPtYuIFSXzZb+1V+Zc590zmznZH4W0dsqeaxcRK0b2g83aKXuuXUSsKJnvHVt7Zc61i4gV4zRnpe4ewtope65dRKwomS/7rb0y59pFxIqR/d6xtVP2XLuIWDGC3PeOrZ2y59pFxAqSu3sIa6vcuc67ZzZzsl/2Wztlz7WLiBUjEC9MoI8hSf8V+BfAi8DfA++NiH/ss9y3gWeBU8BCRMyNvXGWTvZcnzXKm82a1O3ttM4woj3AP42Ifwb8H+C2JZZ9a0S83gXEhpU917WLiKRVkr4m6YvV+EWS9kh6ovp7Yc+yt0k6LOmQpHeM2khrj1OsqjWMIiK+HBEL1eiDwL/tl2vgUuCfdN/nXNuwppDrB4B1Ize8hpVciXwQONgzvh3YGxHrgb3VOJI2AFuBq4FNwJ2S8t4QtMZ07x3XPNhWS9rXM9w85Gb/E/B4z3hvrn8AfFXSfkn/EefahjClXP8q8BcDmwRfrnI97PpfUuv6SdI64OeBjwH/oZq8BXhL9XoH8FfAb1bT746IF4BvSToMbAT+dtTGWm6BOHW69v+XTy51KS7pK8BlfWbdHhH3Vsv8Z+Bi4N/TP9c/B/wJcD3wdWCXc20rNYVc3w4sAJ8dsJprIuKopEuAPZIej4j76zZwsbo34X4X+DBwQc+0SyPiGEBEHKsaBLCWzqVU15Fq2hmqCngzwJVXXrnCZltGcVq88Hwz3UNExLVLzZe0Dfg3wI3AOT2zenP9iKRLIuKEpKc4M/99c12t29m2l0wh1zcAb4uIGLCOo9XfE5LuoXMyNHQRWfZ2lqQbgBMRsb/mOtVn2st2JiLmI2IuIuYuvvjimqu2zCLEqYVVtYZRSNoEfBS4NyK+OmCZ8yVd0H0NXEGncJzR5P774Wzbj0w4178J3BgRzw1YZnGu3w48Osp261yJXAPcKGkz8ArgJyT9MXBc0prqKmQNcKJa/gidA65rHXB0lEZaSwQjH0g1fQJ4FfBuSe+mc+mPpC8AZ1d5/jHgf1V/HwS+Afy/nnU411bPZHN9Hp1bVAAPRMQHJF0OfDoiNtN5WOSeav7ZwOci4kujbHTZIhIRt1E9KibpLcBvRMS7q2eStwF3VH/vrd6yC/icpI8DlwPr6RyEZkuKEAs/HP/BFhE/2Tvek+tf6uY6Iu6Q9EfARRHxYUlX41zbEKaV657pR4HN1esngZ9ucrujPJh8B7BT0vuA7wDvAoiIA5J2Ao/ROcO7JSJOjdxSawFx+tTUv//qXFvDisj12KxozyLir+g8hUVEfA9424DlPkbnSS6z+gKYzGX/mZt1rm2cppTrSclbHm32nBY870haMslznXfPbDYtLL+I2cxJnGsXEStH54cXzHJJnmsXEStH8oPNWip5rl1ErBwB/HDajTBrWPJcu4hYOQJ4YdqNMGtY8ly7iFg5kl/2W0slz7WLiJUj+cFmLZU81y4iVo7kB5u1VPJcu4hYOZIfbNZSyXPtImJlSXywWYslzrWLiJXjNPD8tBth1rDkuXYRsXIkv+y3lkqeaxcRK0fyg81aKnmuXUSsHMkPNmup5Ll2EbGyJD7YrMUS5/qsaTfA7CXdM7Y6wwgkfUTSdyU9Ug2bByy3SdIhSYclbR9tq9ZayXPtKxErx2ngBxPb2n+PiP82aKakVcAngeuAI8BDknZFxGOTaqAlkTzXvhKxcgRwquYwfhuBwxHxZES8CNwNbJnIli2X5Ll2EbGy1L/sXy1pX89w8wq3dKukb0j6jKQL+8xfCzzVM36kmma2colz7dtZVo6VPcVyMiLmBs2U9BXgsj6zbgc+BXy02uJHgd8BfnXxKga00Gxlkue6iCKyf//+70s6NO121LAaODntRtRUYltfveTcBh+FjIhr6ywn6X8CX+wz6whwRc/4OuDoStvhbDeu1HYOznbCXPcqoogAh5aqvqWQtG8W2gmz1daXTKh7CElrIuJYNfoLwKN9FnsIWC/pKuC7wFbgXw+xOWe7QbPSzjPkzPVLSikiZh2TeZ7+v0h6PZ1zxG8D7weQdDnw6YjYHBELkm4F7gNWAZ+JiAMTaZ3lkzjXLiJWjgl9szci3jNg+lFgc8/4bmD3+FtkqSXPdSlFZH7aDahpVtoJs9XWjgB+OO1GNG5W/ju4neOSM9cvKaKIRMRMBGNW2gmz1daXdJ+nT2RW/ju4nWOUMNe9iigiZkD6juqspZLneupfNiypfyJJV0j6S0kHJR2Q9MFq+kWS9kh6ovp7Yc97bqvafkjSOybY1lWSvibpi6W2ccWCTvcQdYbCOdcjtTdXthPlup+pFpGeflyuBzYAN0naMMUmLQC/HhE/BbwJuKVqz3Zgb0SsB/ZW41TztgJXA5uAO6t9moQPAgd7xkts48qU1T3E0JzrkeXKdpJcDzLtK5Gi+ieKiGMR8XD1+lk6QV5btWlHtdgO4J3V6y3A3RHxQkR8CzhMZ5/GStI64OeBT/dMLqqNQ5lQb6cT4FwPKWW28+S6r2kXkWL7J5L0GuBngL8DLu1+iaf6e0m12LTa/7vAh+l8jamrtDauXJ6Drdh/88JzDRmznSfXfU27iBTZP5GkHwe+AHwoIp5ZatE+08bafkk3ACciYn/dt/SZNvV/4766j0LWGcpW5L95ybmGxNnOk+u+pv10VuP9uIxK0jl0DrTPRsSfVpOPd7sUkLQGOFFNn0b7rwFurH5w5hXAT0j648LaOLwZvS+8SHH/5jOQa8ic7Ry57mvaVyIv9eMi6Vw6H5LtmlZjJAn4feBgRHy8Z9YuYFv1ehtwb8/0rZLOq/qiWQ88OM42RsRtEbEuIl5D59/rf0fEu0tq49C6fQzVGcrmXA8hbbbz5LqvqV6JFNg/0TXAe4BvSnqkmvZbwB3ATknvA74DvAsgIg5I2gk8RueO5i0RMa1zjllo49KSfLPXuW7crLSzvyS5HkQR5d1CtHbSK+eC1+2rt/DXtH/menO1Vsqe62l/JmJ2phl9QsVsSYlz7SJi5UjePYS1VPJcu4hYOSb04z1mE5U81y4iVo7kZ2zWUslz7SJiZUl8sFmLJc61i4iVI/mjkNZSyXPtImLlmNCP90j6E+C11eirgH+MiNf3We7bwLNVqxZm7dFLK0TyXLuIWDkm91vU/6r7WtLvAE8vsfhbI+Lk+FtlaSXPtYuIleM0E/1hnqo7kH8J/NzktmqtkzzX0+47y+xMk/3xnn8OHI+IJwbMD+DLkvZLurmxrVr7JM61r0SsLPV74VktqbcvifmImO+OSPoKcFmf990eEd0O/G4CPr/ENq6JiKOSLgH2SHo8Iu6v3UKzrsS5dhGxWXVyqQ8EI+Lapd4s6WzgF4GfXWIdR6u/JyTdQ+dX81xEbJxmLte+nWVtdS3weEQc6TdT0vmSLui+Bt4OPDrB9pkNY+K5dhGxttrKokt+SZdL2l2NXgr8jaSv0/mNij+PiC9NuI1mKzXxXLsreCuG9IaAr9Zc+pUz12W2tVP2XPszEStI8q/2WkvlzrWLiBUkeU911lK5c+0iYgXJfcZmbZU71y4iVpDcB5u1Ve5cu4hYQYKJ9g9hNhG5c+0iYgXJfe/Y2ip3rl1ErCC5L/utrXLn2kXECpL7jM3aKneuXUSsILnP2KytcufaRcQKkvuMzdoqd65dRKwgE/71HrOJyJ1rFxErSO7Lfmur3Ll2EbHC5L3stzbLm2sXEStI7jM2a6vcuXYRsYLkPtisrXLn2kXECpL7KRZrq9y5dhGxguR+isXaKneuXUSsILkv+62tcufaRcQKkvuy39oqd67PmnYDzH6ke8ZWZxiepHdJOiDptKS5RfNuk3RY0iFJ7xjw/osk7ZH0RPX3wpEaZMnlzrWLiBWke8ZWZxjJo8AvAvf3TpS0AdgKXA1sAu6UtKrP+7cDeyNiPbC3GjcbIHeuXUSsIN0PIOsMw4uIgxFxqM+sLcDdEfFCRHwLOAxsHLDcjur1DuCdIzXIksuda38mYgU5dh98ZHXNhV8haV/P+HxEzI/YgLXAAz3jR6ppi10aEccAIuKYpEtG3K6lljvXLiJWjIjY1NS6JH0FuKzPrNsj4t5Bb+vXrKbaZO2UPdcuIpZSRFw7xNuOAFf0jK8DjvZZ7rikNdXZ2hrgxDBtNFupEnPtz0TMfmQXsFXSeZKuAtYDDw5Yblv1ehsw6AzQrARjzbWLiLWOpF+QdAR4M/Dnku4DiIgDwE7gMeBLwC0Rcap6z6d7Hpu8A7hO0hPAddW42VRNK9eK8C1fMzMbjq9EzMxsaC4iZmY2NBcRMzMbmouImZkNzUXEzMyG5iJiZmZDcxExM7Oh/X/l5qoSf/dyzwAAAABJRU5ErkJggg==\n", "text/plain": [ - "
" + "
" ] }, "metadata": { @@ -80,10 +119,125 @@ } ], "source": [ - "plt.imshow(n['U'], vmin=-10, vmax=10, cmap='jet')\n", + "plt.subplot(1,2,1)\n", + "plt.imshow(w['U'], vmin=-10, vmax=10, cmap='jet')\n", + "plt.colorbar()\n", + "plt.subplot(1,2,2)\n", + "plt.imshow(w['V'], vmin=-10, vmax=10, cmap='jet')\n", + "plt.colorbar()\n", "plt.show()" ] }, + { + "cell_type": "code", + "execution_count": 38, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "----------------------------------------\n", + "/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc----------------------------------------\n", + "Mapper: netcdf_cfBand : 1 model_windspeed\n", + " colormap: jet\n", + " dataType: 7\n", + " grid_mapping: stereographic\n", + " long_name: Wind speed\n", + " minmax: 0 18\n", + " name: model_windspeed\n", + " NETCDF_DIM_time: 43653.70576388889\n", + " NETCDF_VARNAME: model_windspeed\n", + " short_name: W\n", + " SourceBand: 1\n", + " SourceFilename: NETCDF:\"/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc\":model_windspeed\n", + " standard_name: wind_speed\n", + " time_iso_8601: 2019-07-09T16:56:18.000000\n", + " units: m s-1\n", + " wkv: wind_speed\n", + " _FillValue: 9.969209968386869e+36\n", + "Band : 2 U\n", + " colormap: jet\n", + " dataType: 7\n", + " grid_mapping: stereographic\n", + " long_name: Eastward wind velocity\n", + " minmax: -10 10\n", + " name: U\n", + " NETCDF_DIM_time: 43653.70576388889\n", + " NETCDF_VARNAME: U\n", + " short_name: U\n", + " SourceBand: 1\n", + " SourceFilename: NETCDF:\"/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc\":U\n", + " standard_name: eastward_wind\n", + " time_iso_8601: 2019-07-09T16:56:18.000000\n", + " units: m s-1\n", + " wkv: eastward_wind\n", + " _FillValue: 9.969209968386869e+36\n", + "Band : 3 V\n", + " colormap: jet\n", + " dataType: 7\n", + " grid_mapping: stereographic\n", + " long_name: Northward wind velocity\n", + " minmax: -10 10\n", + " name: V\n", + " NETCDF_DIM_time: 43653.70576388889\n", + " NETCDF_VARNAME: V\n", + " short_name: V\n", + " SourceBand: 1\n", + " SourceFilename: NETCDF:\"/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc\":V\n", + " standard_name: northward_wind\n", + " time_iso_8601: 2019-07-09T16:56:18.000000\n", + " units: m s-1\n", + " wkv: northward_wind\n", + " _FillValue: 9.969209968386869e+36\n", + "Band : 4 swathmask\n", + " colormap: jet\n", + " dataType: 1\n", + " grid_mapping: stereographic\n", + " long_name: Swath mask\n", + " minmax: 0 1\n", + " name: swathmask\n", + " NETCDF_DIM_time: 43653.70576388889\n", + " NETCDF_VARNAME: swathmask\n", + " short_name: swathmask\n", + " SourceBand: 1\n", + " SourceFilename: NETCDF:\"/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc\":swathmask\n", + " standard_name: swath_binary_mask\n", + " time_iso_8601: 2019-07-09T16:56:18.000000\n", + " units: na\n", + " valid_range: {0,255}\n", + " wkv: swath_binary_mask\n", + " _FillValue: 0\n", + "----------------------------------------\n", + "Domain:[579 x 528]\n", + "----------------------------------------\n", + "Projection(dataset):\n", + "PROJCS[\"unnamed\",\n", + " GEOGCS[\"WGS 84\",\n", + " DATUM[\"WGS_1984\",\n", + " SPHEROID[\"WGS 84\",6378137,298.257223563]],\n", + " PRIMEM[\"Greenwich\",0],\n", + " UNIT[\"degree\",0.0174532925199433]],\n", + " PROJECTION[\"Stereographic\"],\n", + " PARAMETER[\"latitude_of_origin\",67.27],\n", + " PARAMETER[\"central_meridian\",8.04],\n", + " PARAMETER[\"scale_factor\",1],\n", + " PARAMETER[\"false_easting\",0],\n", + " PARAMETER[\"false_northing\",0]]\n", + "----------------------------------------\n", + "Corners (lon, lat):\n", + "\t ( 4.41, 68.40) ( 11.45, 68.41)\n", + "\t ( 4.75, 66.04) ( 11.13, 66.04)" + ] + }, + "execution_count": 38, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "w" + ] + }, { "cell_type": "code", "execution_count": 6, diff --git a/openwind/sar_wind.py b/openwind/sar_wind.py index 3802492..cc6bd4d 100755 --- a/openwind/sar_wind.py +++ b/openwind/sar_wind.py @@ -208,8 +208,7 @@ def _get_aux_wind_from_str(self, aux_wind_source, *args, **kwargs): wind_u_bandNo = aux.get_band_number({ 'standard_name': 'x_wind_10m', }) - self.set_metadata('WIND_DIRECTION_SOURCE', - aux.get_metadata(band_id=wind_u_bandNo)['SourceFilename']) + self.set_metadata('WIND_DIRECTION_SOURCE', aux_wind_source) wdir, wdir_time, wspeed = self._get_wind_direction_array(aux, *args, **kwargs)