xarray_tutorial/grid_soiltype.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/manmeet/anaconda3/envs/py35/lib/python3.7/site-packages/xarray/core/merge.py:16: FutureWarning: The Panel class is removed from pandas. Accessing it from the top-level namespace will also be removed in the next version\n",
      "  PANDAS_TYPES = (pd.Series, pd.DataFrame, pd.Panel)\n"
     ]
    }
   ],
   "source": [
    "import xarray as xr\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import random\n",
    "from scipy import stats\n",
    "from datetime import datetime, timedelta\n",
    "from matplotlib.axes import Axes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_dir = 'soiltype.nc'  # daily\n",
    "ds_soil = xr.open_dataset(data_dir)\n",
    "\n",
    "data_dir = 'grid_spec.nc'  # daily\n",
    "ds_grid_spec = xr.open_dataset(data_dir)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x7f1c4a39d290>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAEKCAYAAADNSVhkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXd4VFX6+D/vTDLpoYYAAQwlVAULCvaCQKKu6NqQtbvLD3tX1HXVtawVLKuyrLqrrq5l1RW/GgL2tVAtqCAQIfReQhLSZub8/jh3kkkyM5nJtJTzeZ55cu899U7u3PeUt4hSCoPBYDAYYokt3h0wGAwGQ8fDCB+DwWAwxBwjfAwGg8EQc4zwMRgMBkPMMcLHYDAYDDHHCB+DwWAwxBwjfAwGg8EQc4zwMRgMBkPMMcLHYDAYDDEnId4diATdu3dXubm5zeZbtXQNAIMPGxDlHhnaEua56FgsXbp0p1IqK5w6Jp6YpnbtdgXX3rLqIqVUfjjttUfahfDJzc1lyZIlzeYbbzsHgPlL3op2lwxtCPNcdCxEZF24deza7WJRUb+g8tp7re4ebnvtkXYhfAwGgyGWKMCNO97daNN0KOEz321GtoammOfCECoKRa0KbtnN4JsOJXwMBoMhUpiZT3gY4WMwGAwholC4TDiasDDCx2AwGFqAGyN8wsEIH4OhA+DR6LMlJesLNkHsdgAkI71pAbe1pJSQoI+dTgCUywUunaZqalA1NfX1OBw6HcCtkMQEXR6QJEfdMUqBzQZ2y8zQbgcRq182sEnduRJBqmt0mssFVdW495XVdbOo4uWWfylhoACXET5hEVcjUxG5QUR+FpGfROTfIpIsIl1FZL6IrLb+dolnHw0Gg8EXblRQH4Nv4jbzEZEc4FpguFKqUkTeBCYDw4GPlVIPich0YDpwWzT6MMExBVtaKmCN4mqdKGdtXXpLtKDG28/DlphQPwIE5tW+rtNs59SNEhuj3Pohne96I+Q2m/TBGuUm9OoJjkTc23fWpcVrpGiIDePt5wEgdjv2fjn1M4p2SsGAm/VBdTWqrByAufv+EfV2FVBr9nzCIt7LbglAiojUAqnAZuB24AQr/SXgM6IkfAyGjoJ9xGBUUiIqoX6xI1KvTtXob6TrJyOpwamtsqs+cLqhsjpSrYSEQplltzCJm/BRSm0SkceA9UAlME8pNU9EspVSW6w8W0SkR7htTUy5UB+MGIhU1+JOSwpcoB1T0O96faAUqrwCV2lpXZqxd2lbjLedQ0J2D1SVfgG7vPZCDFFGgcvInrCI57JbF2AS0B/YC7wlIheEUH4qMBWgX7/g3Fw0xj12BNTo5TH7nv1QXoG7dF+L6vIQaNksVi93SUgEQFVXIzYbtmxLftfWBihlaA/YEj0b/EngVri6ZwCgEvWMR9mtjXy7/+1ecbkDpscbcWqFBynZhLOsXuDGcvCkPRwYwiGeT9jJwFql1A6lVC3wDnAUsE1EegFYf7f7KqyUmq2UGq2UGp2VFZaPQIPBYAgRwRXkp9maRPJFZKWIFFv73I3TRUSestKXicihzZUNpLglIrdb+VeKyESv6+eLyI9WG3NFJKo+6eK557MeGCsiqehlt3HAEqACuBh4yPr7Xtx6GAMmOKYAYMtMR5KScG3fUZfmUVSIFfld/wCAu7xCt1/zWth11jntNEt6AZngmIK9a2d9kpiIsmbg7qpqxJGIqqmftUbi/2IID61wEL4yh4jYgWeA8cBGYLGIzFFKLffKVgDkWZ8xwHPAmGbKTseH4paIDEcrdo0AegMfichgQIAn0QpgO0XkEeBq4J6wb9IP8dzzWSgi/wG+BZzAd8BsIB14U0QuRwuoc8Jtq6jylbrjgj7XIskOAOzfrqrP1LMH7t17m+Rvi/h7OU1InIy9V08AXH26g+pGwtY9OrGisoGGnqENIDbce7WQKqp8hfH283DXanuchKzuUFaObddunTUzA5IcuDtZ2p2JdtxJ+ufvdthI3F2JrN0MgGvvXkT0oog9MwM1qC9VWSm6HgViaWYqmxCB92/IuJJ0X1L2ZZGQnUXhqodj3gdt5xORmz8CKFZKrQEQkdfR2xHewmcS8LJSSgELRKSztSqUG6DsJHwrbk0CXldKVQNrRaTY6sMStABKE5FdQCZQHIkb9Edctd2UUncDdze6XI2eBbUb6lSf++Tg7NMdW7Uexdp2R2aDOL/rH5CsbgA4u6WjEvSPoqaLFrK2Gv2ycIRY78T0iykqf6nF/ZqYfnF9H0f+Eamqxb15a901W+dOuHfuqjsvqnq1xW11JPJH3NnkWkH2FXHoScfGHRnJmwNs8DrfiJ7dNJcnp5my/hS3coAFjetSSn0jIlcAP6JXn1YDV7X0poIh3qrWMadw41Px7kIDbJnaulx5KTtEa7lt/+mjm45U+6Q2yZe8s4aEJSsj0qbtkOFIrZlRNYfYBJdnlmK3Y+ubA4AtORFn1zQSN+tZOTZbnYJvQu+eqIr92pMAYB+Rh1TrmQ+VVSiXC8mxZrrpydhqnNqDgIW9XGvJJe4oxblxE+B/edSzPAzU2cKJ3Y4kJWHrlQ1A+Ygs3AnRnwo5yvTz5D2QiTUhzny6i4h3wLHZSqnZ1rGvShrr0fnLE0zZxvgsIyKJwBXAIcAa4Gm02cv9zdTXYjqc8Ik13qP/YCnImgbdu+Au2Vh3zV1dBdZSiBw+om7JRNmExMh01Sfjj3kAcbqw79b7QJSVo7K74UrX6uq2ypq65Zq5u/8exZ60PyYkTgag5uRDaS3K/wVDpuNaUwLEfs+xLaEQXMHra+1USo32k7YR6Ot13gdt7xhMHkeAsttEpJc16/FW3PJX18EASqlfASyj/ybKD5HECJ8YoI49BIAagcRvVwP+rbALsqaFXP/+YwdTnenbc0JLqO7qwH7QoLDVbYvKX9IeHxx6wa/2yOHY19YHkXRt078H85KD0rMOxW65MEveUY292jIB2LGPxJV7UJVVAEh6GnTN1BntNtzlFditJdea7HTEqQe+4srAmdoLx249u7H9shbXQQNxO/RzklBegztF/1/cfbpDn+4k7Cr3279glBwKBt3C7iN7hXjnoZG8y0nycuv92qUzhRuejGp7gYjQsttiIE9E+gOb0MoAUxrlmQNcbe3pjAFKLaGyI0DZOfhW3JoDvCYiM9AKB3nAIiAbGC4iWUqpHWglhhWRuEF/GOFjMBgMIaIQalT4Az6llFNErgaKADvwolLqZxGZZqXPAj4ETkErAOwHLg1U1qr6IXwobll1v4lWSnACVymlXMBmEbkX+MLyOLMOuCTsGwyAET4WExIno9wKe2ZG3TVVWQmAu0YPSUNRF/b42LKlJAddpmD4HUHnbY2cfPyDuJLsuBP1iNDtsBHM3Y874UEcJTugQn/f2G1gt6Oq612nzN05209p/xTkXNPgvHDT00GXzc+8VGuIgd5XSdYLY06PnzxVb2LYXtXITz7uARK3laG26hnq3NIX49yj1oM2Mo2MmaRS6kO0gPG+NsvrWOFn899XWev6LvwobimlHgAe8HF9FjCraYnoYIRPlFHDB1LRW7+4On+5HjplBlwq2PmbwQC47ULSSG3j5UwWlA08Ay0VZWeRzmQ7ld0au9nP8JkXUsnokUbCfqfP1Iqzj8CZovvrShQYWK/Ik1ClSN9QFYEet1GsPbzk9aUk7XQga/WmP717sG+otglUuam4EwVnUtP/edp2J8mbt+G0hENSSnKdI9HaXp1ILHfCd1pjV3pksb9XMlVdPKP1+mGBKEVCFdSOSAPAPrYHmWsrQ74d95ZtKLtednMlRu4ZtdfqpcT0TTUk7qnCnd2FoiX3RKz+lhIhVesOS4cTPgVZ08CajZQf1qdO+6upzpePspZfNNcWrWXjvVfhUacGsKUGrq0g9wYAag/IwlapX9r2zTsh00dcFQMA44/2UroRQercIpWzP08LaUepnqHaLI2vYMel423nYHM4kCEDANh7UBc6h9C343/zaN2xY3c1CTv13onasQvprp1gutOSUckJuFL0T66qayJpIbTRGsjvPhVo2Sy0vaGU4FKt1wVRW6DDCR9/7Lp0DK5GKkdirax0XV5F4s5yKK0Iuj5bSjJF5S8xIXEy6Uv1yM2p3NrVfY6fTVml8DzPrmTYn9w2Hu69gxKhBTp3rmShurM2GOSgpv75kkr195ZZvL/FffP2GBEKWy8a6VMp1V41CHuNQiyF1sRyd93IPGQOHgqA2wa1nZJIqtVq0bK3jPQN1sPoUpT3T6M2VXfGUeYmebcWvDUZdqrHH0zKl78A4Fy7HjniQN2vrfvYN6oHcrpWsspYuomM+cvJsEJ6lI0bimOfFtLJP23EtX0HtmGDACjP68zOUZ4BVCrizib7jeb3nosqXmbCEX8GYPsYfzPl0LA5AUsRI2n5Jlw7dgXMH0vcZuYTFh1K+ExM/h22jMj8KACOnPI4nT8v8Zs++vIZdI1Ya4ZwmZh+cd3ejWu39uzQXvdrosnE0fcAIGv0MmFHVLHXCgcd6vUZcTrct7e7YDD7e4Q2Ytl+WDJ6jVwv74j7ANK2NvRp61Gnti/8GfeoPPYM1SPH7dPGkFiuR8Y1nQS3HRKs5XRng9W5DBIqwWW5IYiAIk2bx/N/2t/D/wKVuDOoX/1oPPvqhKNUL6Wl7HKRNu9H33XY7dgO6EvZoE4A1GQILj+aEnp27PX89Gj8j0qAugW1bL/9BtiX22ggdJDHB0V98N5Oa104Sp049uk2bTVuyg6wPFc4IbECqo7SMyiUYm+eTsv6zkZldxu1VldK+/ejc3H9vty+A2ykb9I//0Rrdig7tEBOX7eZDEvZYtM5uSgbbD9nGAApO1ykl/hXybatLAHAddxBAe/dHxkb3HRatAnK9WzXo3QSiwBxoRBJhYOOSocTPgaDwRAJXPFwbNeO6BDCx/i9MhjaPuNt51B1+hGkbNMakhU5KdSm69lHQqUis/AnIFZhtEPycGDwQYcQPh52XjSaykZxUcXyVJ+9pBp3orBjlF62UM3sn5cO0A+eM7k/brtQo1dscB3b1ItGTaeGIyRXiu86/S31GELE2v/PWuYkoUz/gyt7JLHt0oOxW5rdWV/vrMtuGzaI7WO7kbpDb+TbahWu5NYxqt09zI62H/RNWT/fD+qmE5pqXO7JS6DvB3qJLXl7KvZfLG8TGRm4SvdBul6jk4QEtufnAvVLwyk79ZeasqOG0iF6Sa7L7k5N2nBbnhiS9ir29wzwHbqh50K9pGavdGIv1WvRrl+KcRI4KGNrwW203cKi3Qofj5sa5+A+7fcmDQZDXNCORY3wCYcO8V5eO00bbjZWTkndAj0Wam/BUuNkY353nCGa2uzr3zpGyIZ6krRzaBwfakfCnlF0fualuKv0aNtNvb8ytWotPXaVQpoe5rtHZlF2QPvQ+PDM9HouqiZp9TZwaUWZ/cO7kl7dGwDbJksdvaxekaDHF9sA2H58NhW9oXSgZSjsSKHbPxbqY5r65fN8pxMcU3D+v8ObrDR4yFyrKDvAweIXbwz6Xua736Ig+wr2j+6v208SOq3WigkJv27GWe5fESLSKIRaoxUUFu1S+Iy9YIaXvpDBYDBEFqUwRqZhElfhIyKdgeeBA9Ez2cuAlcAb6Ch9JcC5Sqk9odb96/VDAHD7WBJPtcKAZP19IVhGd8V/9ufx3BBrxAW9FjjrYsPsHmYne7G2NNwwrmlIPM9MJ+v7GlJ+3QVOrVLsdjTM628juqjqVY498zG2HBnESNYTRcXbrrSVTn5zPq8laf73deeFjTxTT0y5ENCzQFtyElt/mweAo1whLn2DPd5bBW6F+wBtAMsvayEpqdkgg/NqXiO/02VsuVirXO/vVW8+4NgLXV5ZpE9CdBdXuO250ApEDTFGpmESb9H9JDBXKTUUGIV24e2JPZ4HfEyUY0oYDAZDqCj0zCeYj8E3cftmRCQTOA54AUApVaOU2ouOMe4ZVr0EnBGfHho6CvmZl8a7C22O/IPuJP+gpuG8OxIubEF9DL6J57LbAGAH8A8RGQUsBa7Df+zxBojIVGAqQL9+Df2CJe9y4szwv4RS1U0/EGrMQaiv9bJE7geVbDwhhZpubr/lDOGRvM2Gvdp3Wm0G2C3PD32f137EPA4sxx99f53HiNwPKkncUb+x7E7VS2vhejlO//pXespAQC/zVWU3fA4SS/Uzc8D7ZewemU63paUA2HbswdVbe77Yb4Uk3z1cP3tVPeL3LCXutZE071vAf7C+ospXGpwPu2smAKVD69cUt48ZRL9CrYL+WeGtFPS9DtU5OK0c7xAM423n1HkBSdxSCv0PoHDVw0HeTetDIZEKJtdhiafwSQAOBa5RSi0UkScJYYnNioE+G2D06NEt9OxoaM2MmK5fhr3j3I+2zsTk3wF6b8sf+d2nwhXDAtZz0E0z6dO43Mg/MnfZ/T7ztwY8jk7LBqWR+c63EatXAbXGt1tYxPPb2whsVEottM7/gxY+/mKPB03Kii30SssBIPOH7bjTtVXnqsszUQmKnM8tN7nfLKsrk7hmGxkDD2BHLyPHIomt0kbfT/TI2fHBAsC3M09P8D3QKry2Q+pfhKsvSGborHqdk8LlD0a0j40VESYkTgaxsXeKVkKp6ipUW+qTtuINdC8O7EzTcy/2g4aw8nJtiKkSY/tc1XZ2s33aGLL/viToMgmW0/bKnIZ93TAugdQtepS/6roDcCfp9CF/D14PqP05cJWIxfMRkXz0/rcdeF4p9VCjdLHST0FHMr1EKfVtoLIi0hU/ilsicjtwOfqndq1SqkhEMoD/eTXbB/iXUur6iNykD+K2IKmU2gpsEJEh1qVx6NCuntjj0DD2uMEQVSY4pjDBMSXe3YgaYy+YQUHONXWfSOGZWbU28jtdFrW6FdrDQTCfQIiIHXgGKACGA+eLyPBG2QqAPOszFXguiLI+Fbes9MnACCAfeFZE7EqpMqXUwZ4POoz2O2F8Rc0S73njNcCrIuIA1qBjk9vwEXs8FApLZtYdD79zJn0f16O/5O2j6bGkluQfSnRi964oK9jX8hs6Iy43OMyeT9jU6hFh5x8T6PWeduFSuP6JgEW83amcfNwD/HKRA7Fr1ziJWxNRv66PUmebYs/th7tTKt2+1oaWzuI12NPr9zma8x3mfS8FvXT043WXD6Iy241Kd0Whx01xJ8HeQ13snaH3WZJ22emy0v/sa+7O2Rzz28cAKDuo4W/AleSmrFvTMmvP6UbuA8HPrOJBTTfts2rL8YpOy/WeHj9Epu4IzXyOAIqVUmsAROR1tNLVcq88k4CXrXDaC0Sks7UqlBug7CTgBKv8S8BnwG3W9deVUtXAWhEptvrwjacxEckDetBwJhRx4ip8lFLfA74MbHzGHvfH6p82UpB3K4WrHwH0CFYSrVu7/pAwe2mIBAU9r8Rdrtd1mrMRiTX5XX4P7o436Dj67Mfqjl3JtrBeBsed/iipn+n3ZWsIfzD+6PsDeMQLH6UkUr7dcoANXucbgTFB5Mlppqw/xa0cYIGPurw5H3jDEnZRI94zn6jjKIWqk0YB0H2ZE0dpDTt+o43patOFfQP1S8eWUhO3PrY3VI227M36tgJnn+4krNkcWgUi2FKcJK7To9Z+RZW4DhvCR5/fEemuNsDeRe/PONfo2VpRBJxbFm55BoCCQbeAW7HiT1kASGKMhJ3DinraqZZtA2DrRP2ekb2J9P684btl4xnaOFfswb1zaga4WPnEIfT7sHXuk66YMZRhM7RW4uAXkig+39q4i8DMRyscBC3euouI9xRxtqUwBb5NlBt/of7yBFO2McGUmQxc2Ew9YdPuhY/BYDBEHgnFgHSnUsqfC5WNQF+v8z5A49GavzyOAGX9KW4FbM8ye0lQSi1t9q7CpF1ZQBXk3ao1lQytmvG2cxhvC3krz2BolonJv4uJAoRWOJCgPs2wGMgTkf7W3vdktNKVN3OAi0QzFii1ltQClfWnuDUHmCwiSSLSH63EsMirrfOBf4f4dbSIdjHzqerZ8DZ2Xno4u8fqzWpHRjn7T2hcoqLuqJmwPYaWkKSXb9bdoE/ty/QyZ8o26Pb3r5stbl+0AvdFBzHg0Z/rrs3d83zk+xmDuj24N27GlpHBsL/oZ3XVH3pgz61oplQUSall61gdv2fwzLWUHp2LI7U29HqSnGyaopefUoaMjGQPw2J9fhpdvgM2bgGg5NaDsA2MrNfrSHgvUEo5ReRqoAitLv2iUupnEZlmpc8CPkSrWRejVa0vDVTWqvohfChuWXW/iVZKcAJXKaW8tWDOtdqKOu1C+BgMBkMsiaSHA6XUh2gB431tltexAq4Ktqx1fRd+FLeUUg8AD/hJGxB0x8OkXQif5K0uKk/QuqBrb8siPauUtDj3yVBPZrHeXK9Js5GQ7SfAixd1bl8uj2avYktj7wJHTnmcfXmxUbv2R2LePgBWPJRDcloFiXb/ShAVu7XroO5fJSBe2RIrFJ2Wajfxe8bGPxRv4+85/7U/ANDtZzdJ/9OumFZHqC13+9q1iDntQvgYDIbWw/D/3gPA8jPuiWs/oolSUOs2wicc2oXwqe0LZdfpUVyPZjUNDbFkd2k6/ZZakTJdLpzbQvaW1C5Jf30Bnb8bBMC6vySRlhQ/Vf/U7s23XWHXs5rKbkKfJ720hhMSWPmANmXoMngXldWtaxfV4wbpxAkPs+0IK75TUfj16mU3I3zCoV0IH4PBYIg1kfLt1lExottgMESFgtwb4t2FqBFBVesOS7uY+dhtijRH5Jct1m/UcVqSSxwkHbabTilVEW+jPVPjTKD/s4ra7Aw++uJOJh52N7ZDhlO09N54dy3ueHt5Hnz/TFKODNELRIzpl7MLgPXu7th6aqUR1Skd2byDwfdotzqS5EBla8UfKY2sWnO4OL5YRpfMSLraMstu4WK+PUNMmJhS760jv+sfyO/6hzj2xmAIHzcS1Mfgm3Yx87GJIi2x4czn1+161lKzK4W8IZsA2P9MH7aNtpF3ZEmzdZbs6sqw29bVnZcePxCnpctQ0ctG5bF6ZDewx84I3EH7Yt1u7UOr35+cQCVqtf4eN47vTN/ZywOU7KAITZ7f1oLTbaN4Uw8yF2mFg4wE2DxJRw7uPX8HZHVlboTjK0WDoqpXtQPZCKG13aLpurT90y6Ej8FgMHgz7oQH2XBNQzuqSFpPmjDa4dMuhI8dRXqCHjnuub0v4nTT/wdtSmbrl0PJ2XqkVnNWJYPvLaP6o14ArJli4/Aha33WmfRJJlvOyeT7p/WmaX7mpbgq9gOQBiS8pkf3Gy4ZwrAzV0Xt3toC2yrTKX1fB7tO2eHmgMX16tSFqx5m5A06vlLVIfupOG5wXPrYmpFhZXXPb2tj0XL9yv7hybavPLD1X/p3z2mRqc8sqYVH3Pd8RMQuIt+JyP9Z511FZL6IrLb+dol3Hw0Gg8Ebo+0WPnEXPsB1wAqvc5/hX+PFqGtnkjv70Xh2wWAwtEIiEUa7IxPXZTcR6QOcinZyd6N12V/4V7/sr01kx339AUhaqCNFeXw8HX7ZDNI3Kha/eCMTD7mbFbd2IftTfdvD/rSO9He0+vQPO3LYvS2T1GJtoX3Af0vYfGZuXRtd5iXy3eYRAKifM0jQK3BUZSnSE9u/CvaXJToEsf2nNI48bVmDtE1PDiINN1+/cZPPss7jdUCv/jMSgRo++uLOqPa1reFck05639b1DP1633AAhq+wvFNEbq8+JjhWb+Hofr7VvSMRRVspwWkES1jE+9t7ArgV8PZo2CD8KzqWeBNEZKqILBGRJa59cXRNbwhILGKrGAzxwCy7hUfcZj4ichqwXSm1VEROCLW8FYZ2NkBG5z7K7RC+eO+WJvm6rKgXTEXf3ct42znIGB13xLlpM9+9fBQApWOqGHpLMWLX6pOVr6VxXJclDeoan7tSH+SG2tu2zaInDiP3X98A4DrxMLY8ZeeXJ7VfMioSyPvPogClwenU32nC8nUB83VUen/pIv2k+CscbK7SYcS3PT6AlP9bCOiAL94GsW2FXeP7M8gevWCcnj0fQ8uJ57Lb0cDpInIKkAxkisi/8B/+tVnyR90FgLLZqOyXDlbFAAUDbqa2b1ejnxJD8kf+Md5dMBiihhE+4RE34aOUuh24HcCa+dyslLpARB5Fh319iIbhX/1iK68mtWSfz7SdB6eTXbSp7rzi7LEAfP2m3p+Y4JgCQGL5aNz7ylEubRuw46PBZJ6/pkX31t4YevXPfNtDzxCzF1fiHDmAYbfoWcyh87cx97Jj/ZZd/fBh2H7VP9LiW4ZRfFvbV9mNNMnvL+GD48dw/rgv49qP91dq79Sc5qbk7bY32/Gm82tL+CLzcCqzrQsC3X7yrO6/Fnb9kbTzEZF84El0NNLnlVIPNUoXK/0UdCTTS5RS3wYqKyJdgTfQ6zQlwLlKqT1W2u3oaFku4FqlVJF13QH8Fb3n7gbuVEq9HZGb9EG893x88RAwXkRWA+Otc4PBYGhVRMK9jojYgWeAAmA4cL6IDG+UrQDIsz5TgeeCKOtTa9hKnwyMAPKBZ616AO5Eb4UMtur7vGXfTHC0CiNTpdRnaK22gOFfDa2bEwoeoeQ39eOZoXHsi8EQTZTSrociwBFAsVJqDYCIvI7W+PX2QzUJeNkKp71ARDpbWxK5Acr60xqeBLyulKoG1opIsdWHb4DLsH62Sik3EFXfYa1C+IRL1QEOjv+3783F/a4kvl5xeN25Z7nNw7ya+il4wec3496yDYDeD3/Nt69qzwgqI43tx3bHXq2du/32lo8i2v/WztD0LazYodXMP/r8DgAmpl0EwFurD6HnNie22qZB/Ab8s4QNb3aj74e76i8GVJrvmMx3vUFB/xt5+4CDAbhwyCJs4j+kdTSodifSdZFXILhLY9p8xPH8ro89Q9voVXW140yO7B5NhJbdcoANXucbgTFB5MlppmwDrWER8WgN5wALGtclIp2t8/usbZBfgauVUttaclPB0C6ETyh4XppFFS/r8/SLAZDERKRLZ7/lDJo61Wm7capoaL3kd/0DKx4ZTL8o1R/ink93EfFWnZ1taesCPtflGo/k/OUJpmxj/JVJAPoAXymlbhSRG4GsQsOvAAAgAElEQVTHgAt95I8I7UL4ZKeUkWrzraq6tboTe4am6HxvryTQeFLt2UtR5St+0/v/9XEAv221Z0onaYO9Lm/oc4/wBii47RY2zdTfsW1NZzqt0s/30m8OJR3F3GX3x7azbZDCtTPqjk8oeIThD2hD3iGpW6Pe9s/lOaz643AyqeGzue1ravrbv8wHaDCTXPKPyNStghc+O5VSo/2kbQT6ep33ARoHd/KXxxGgrD+tYX917UIrM7xrXX8LrZQQNdqF8GkJJ054mMRPv0OSkhpcz/uLdoLZ59MaXCl6Tbeyqx17jYIjY95Ng8HQSomQY9HFQJ6I9Ac2oZUBpjTKMwe42trTGQOUWkJlR4Cyc/CtNTwHeE1EZgC90UoMi5RSSkTeR+8TfYLed49q/JN2IXxsuMmw+XZPMiptA6NutJZFb4R3tx1C+X196tJ/eVLvZUwd+wVvrj0EfvTfzlXj5gH4bcuDG8HW7Oy3bXHbgfreXz3q1CZphcWPcuLEhwHoRi2fFrWv0XOs+azwVgY9qmdC85Qw5vjljOuyoplSgXEj1KgE7Nbc3y7uumd0wxmdSGEThRufCq/jrYjVdw7j7uPeiVr9SkVmz0cp5RSRq4EitLr0i0qpn0VkmpU+C/gQrWZdjJ6dXBqorFX1Q8CbInI5sB44xyrzs4i8iRYsTuAqpZQn9sRtwCsi8gSwgyjv/LUL4WMwGAyxRXBFRtsNpdSHaAHjfW2W17ECrgq2rHXdr9awUuoBtD/NxtfXAceF0vdwMMLHYDC0G+rCtd9/SNTbCmHPx+CDdiF87KJIs1UHlfeCXgvgWe8rv9Yd3TF0br1xyjkt78/frjiLktMcAPzl1H+3vKJWSGV3h8/rZqktsuQ9uxGA5Xf2Ys2TQzn1z1oBwR5QZaaeO5acSd+XEqjqpn/iSbudACR/qteVbf1y2P6E1ljsnt0uXgMNSM0rDfqd0BKMb7fwaY0eDgwGgyFkDrppZuwaU3rfJ5iPwTftYshjQ5EaxVEOwDXzL0Jq9EjnqVNfCpi35DI3eZdqo9c7XOcz4wytlvyfHaM5NyuwB+jWyoyS8QB0+mZ9nHvSQajVM5UDB29g2ze5PPboZABuu+W1gM/6jW/oPeKBf9JeqT+rfd1v3hHT9ct69YXt5w05+Cs9M7y687vN5AwfE0Y7PMzMx2AwtHkmHnZ3TNtTlsJBMB+Db9rHzEfcUV3fBTj+kBVsP1W7H7l74Ok8Nty/19+/jX2FqY95Qj8qZp10IgAVo3pzx8ARlI3SRqpjBq/lu805deXuGvkBOYl7onMDYbCmpge2P3cDQFWsjXNvOgbeas/HPPEYKdu0ev9t/zuHF098wW85jy1lyT1jWHVnYA/iOTMWA1BzkuXN+vowOtwKuPfKl+ls3x+z9sySWni0C+FjMBg6Jp6QKHJQXszbNtpu4dEuhI8NRZpE1+XNldmfcOWZ1wKQddde3n9Bq3JO7rKwaWaB2ac+D8AfCi+n8+s6mupBKYuYP/souixy8N0zTUelR52Xyol3fgXA6Z2+i8ZthMSf1/0GgJo7slEebw9jY/8j7+h8+c7NHHvmYwAMueZHrrhtGgAvX9TUKLTXAmfQ9dr76Vl34vxvI9DL+GKbuZechL0xa08rExjhEw7tQvgYDAZDrDGq1uERN+EjIn2Bl4Ge6Kh5s5VSTwaKwNeemHj4vQAULY7tRqmhY5Lf5feUjRtKeqG2F5KB/Zroanlv2hctvTeGvWubmD2f8IjnzMcJ3KSU+lZEMoClIjIfuAQdge8hEZmOjsAX0ILRhiJZgl9uaCmX3zQHgPfPGMu7X2sntZec9pXPvMl27QX6vdOebHB98HVbeet3J/ss8/UbN3HYVL285braRoZdbzKflvFDwH59UjGMvxbmA5BeItSm6+v/njYjQCnN0qoDAHju1+Po10kvW/z4aR6575XBd9qvoLCdpK5ddL92t7txQJvgf+/eXHc83qYtoKcecgEvj/xng3z7+gX/k15zcW8AxNUbWy30nRu7ZatIsfk6Havr1dzmn/VIohDcRpMtLOImfKxAR55gR2UisgId6MhfBD6DIS4UDNaPX+Gqh+Pck9CYmHJhwBAhgcjP1PZCu888iK7vaq8Ic/dFKBZBO8FMfMKjVez5iEgucAiwEP8R+PyXB5LF1Vy2sBmXuhKAf408jWEzdYC/T04cxinpPwVdx8u3nA694Yv3bvGZvnS2RxHhBnL/pjeZzz7F/4bw++UH8ebDE+kGLPrnjUB9wLxbJ57FUwPe9Fv2f5UDeXvSMQB0TU5k3WGDAMhEMW/Rn4K+J0Nsme/Wav65zz1G8qiGz/25V+n4NZ8d1q3B9aT3k1Gujag5eoVA/bKGlXfVx1ma4JhC+aTDAEhbV9FmzCf7zNWRnpOvjf7vvwFG4SBs4i58RCQdeBu4Xim1TyS4f6iITAWmAvTOMVE1o43HdYnzmH0ArDjT7FW1BYa9q/duBmfFuSPtETP1CYu4Ch8RSUQLnleVUp7gG/4i8DXACkM7G2DkyESVJLF7Eibc9QXfHJ4JwBv3T+SoB1eTbfcf42eDMw2A6x6+im7vfx10O0nbtVD1dW+f7B8IwNxpx5NJFR99fkddmifKaEHerSR9qqi1xrGJjX4tSVKLs1t63XnWRxswNKW1LrflXbWQS764icTf62in/xzyL87P1Cr6W78+gfeWHUwK+rmcmfs2599+C5zXj4Wv3Nikrnk1r5HffWrd+YYrRgJQMaKail4jo30rLWedDtwZy9+/h4488xGRi4HrgCHWpRXAU0qpl/2XakjcdsxET3FeAFYopbx3Cz0R+KBhBD6DwWBoFSjA7ZagPu0NEbkI7Q/jJnQ01BzgVuA6Ky0o4jnzORq4EPhRRL63rt2Bnwh8BoMhMoy5aAZdvt8NwNxl9zeT2+ATBXTcmc+VwJlKqRKva5+IyFnA62gTmmaJp7bbl+B3X9NnBD5/iAiOIPeKIsHULot5Z1q9wsCNq87lreENtYpc1hLXaY/eSq9/6RDIPVjBXLd/n3CN6f2l9trguKzhvR318bUMfbgMgAQqmPtzk6CEAKg9exn3xi0k79Tl51/9aIP0canr2f38pwAUnjMGlZHK3B9917W/l76fyQMDq30bYsf8Rs/SxMNSWTlNL6PaM2qxOVx8fMTfrFRhx2GQFcCZwdyds6PU0+hRMW4EAA4pinnbkbLzEZF84El0KOznlVIPNUoXK/0UdBjtS5RS3wYqG8heUkRuBy4HXMC1Sqki6/pnQC+g0mp6glLK17ZHZiPBA4BSqkREMoO9b6OobmgRJ05onfsgBkPMUEF+AiAiduAZoAAYDpwvIsMbZSsA8qzPVOC5IMpOR9tL5gEfW+dY6ZOBEUA+8KxVj4ffKaUOtj4+99upF06hpjUg7tpukUCAJImtHP3g5kcAOO7dm6n5sic261/+YcUAHnj/twx+ZDUAvVlFYQtHlAkf6WHqTRtO46/96sO0D36uFndmMvO+uStgec9I1uNNYccVdva6UwA40KG9/16Q+QsAM28eH7Cuziv17GnopM28eNWZod6KIQY09kow3nYOJ92lZ+hHnraMwrMeo+u5+v9Y4ky0crXtZbedB3qUcmI9jpZIKRwcARQrpdYAiMjraFvH5V55JgEvK6UUsEBEOlvKWLkByvqzl5wEvK6UqgbWikix1YdvQujzMBFZ5uO6AAOCraRdCJ/WTn6ny+qO55a+GPX2PHHslUvbPsioIYGyt5gJR94H0KwQNLQ9jrhE6wC5rTFxVlFJXZp3uIcOTfDLbt1FZInX+WxLWxf0Zr23mulGYEyj8r7y5DRT1p+9ZA6wwEddHv4hIi60FvL9lsBrzDAf10KmXQgfARJjvILYyRppLfrtDE6/6UZ+d8E1dWm/fnojBY9eGXYb811vAHDY1JnY7ylkt1sbCNr3hhazxOM/7obvV7D4z9ot0DH3LOCO7vVRVfv39TfD1vT4eBMA9w46j7y1W3B1zwipD4bY02RPKOVCqv8vG4DTemuvBYfHvFeRpedCvS+6baqijz2G9n4KVPCabDuVUqP9pPmqpPEL31+eYMqG0t7vlFKbLHdnb6MVwpooDyil1jXTRlC0C+FjMBgMsSciy24bgb5e532AzUHmcQQo689e0m97SqlN1t8yEXkNvRzXRPiIyFoaCjnxOldKqYF+79YLI3yiwHXfnQ90iXc3DIZmKci+AoDCbc/VXTvgxUfIDvBqGPHePQDUOvVMY9VZHXTZNTLabouBPBHpD2xCKwNMaZRnDnC1taczBii1hMqOAGU99pIP0dBecg7wmojMQNvo5AGLRCQB6KyU2mkZ/58GfOSnz41ncTbgXOBmIOhAZO1C+AiCPYaq1t6kksDoW75l3v/p/0evo/XyVPHTvevyrD7nj2G1sXT2DTy/6h3e367DHTv+XsZ7x/zVb35/ziTtohh792IeG/UG+aPu4uYX0pi3Qi/fDvuLdpvTYJvTi8I12s/c6MtnsOaiXvSf6SejodVSVPkKBbnad+Dif/dn3xn+fzPzHtTulBKt5eXJl5wOwJrd3fyWiQefztNOX48/VfjzM88zOqk6do1HQPgopZwicjVQhFaXflEp9bOITLPSZwEfotWsi9Gq1pcGKmtV7dNe0qr7TfQv3QlcpZRyiUgaUGQJHjta8PzdT593AYiIDb00dwvwPXCqUiroF0O7ED4dkUGPzKT41qbRUINhYsqFyODcyHbI0KYZ9MhMbpn0X+vs1JDKjrefV3fs2ads90TQyFQp9SFawHhfm+V1rICrgi1rXd+FH3tJpdQDwAONrlUAhwXTX0tAXQbcAHwJTFJK/RpMWW/ahfARIIH4ORd9rNc38Id6TcVyVcO39x5anyECPhqefua39PxCx1v59fzOLapj5b4e1EzQdahVJWw4syt5m5YCegjUeIPaF1lz11I1ZQArnwpao9LQiigs0TOagp5XsuJu/T+UTrXYN9Xn+XLiTFIkpUG5/wz6oEld5aqGkx65KXqdDYHPP7iVwfcn8uPlT8eszQ4cTG4t+pXxBHpWNUpERnkSvfx0BqRdCB9DdDhukvaIkLZ8BwCFqx+JZ3cMbQBv56Rt0WNCSLRDv21B8hF67jfK+nijgI4jfCoU/HNfvQLHpPTVdLenxq0/p/14AVy5i28mPNR85iD54cn6Jba8t1pmGFj5x15suF2PdnMf/g737j1BzXa8Kdz0NOOPvp/VeUkt6oOhdVC49dm644nJvwMg9Uy9X9LLnhZUHZ0kmaXTn2HQkGkADL3t52ZKRJfcexdyWKk2efjo+kej/g6IgyPtVoFS6hJ/aZZ/t6Aw7nUMBoMhVIJ1rdPxBNTMYDO2i5mPwT/jTvoLv56nXakMCqHcIVfOxJiRGsLBo4jQPpUQpCN7tQ5E0F9KuxA+W1Z04r/HDa+zVej/9OPMPu15AA5O2kc3W3DLCJHgl9oKMu622psQnTZaqrq95sxkPjpb79uclHt9s/lvmvkqAHf8eEaD6/O/Ck913NC6KKrS/+eZK/QDa2N3SOXXnGnt7ZwJW13lABzz1s0MvGlh5DoZBPNqX687Pvn4Soovt9HtS708nLzbRWWWRympaTC9FtHxZjXBEPS30i6EjyEyTHBo+7TVjx9G5yBXZD1l5tW8FrV+GaJPQb/rmRD7qARtG3e8OxAfRORHfAsZAbKDrafVCp/mYlx4kzeqH4VL6i20115zEzrIHoz/7AaKhjZVE400n1Xp2eZdN98AOfDVf26OepvB8PEntzPgtQcB+PH0v5EiejHt2CGrA5b7+q6nSBK9XHfb9530xUnR66ch/gxK2hp2HT3tOp5Q8eRZnPzBZc3kjh51YeV9xNWUWRGY+XTsYHKnRaISv8JHRN4nwBRKKXV6JDrgp21PnIrxaF9Ei0VkTijWs4bQOGncQ82ORI7/zaOsO0VI2WIn95UNDdImOKaY2U8bpaDXVZCY2HxGQwM6sLZbUI5FReQbpdSR/tIDvW8eC7lXkSOYGBethmkbj2R9vl5bzmBVq7NvWDNFjwIHvg6rT/gnADP6FjbJd+L32gB1bte/AfUvI8de6FLs9Fn3tisOJ/u5xZHtsCHmVPwjkVNTqyJa5wkzv45ofa2ODip8QiA5UKLfhX2l1OdKqc+Bgz3H3tci3ctG+ItfYTAYDIa2QUDxHMyez8XovRdvLvFxLZI0G6dCRKaiQ8rSr1+/KHalfTPeZvn+OTEot04BOfm4B1h/rZuBN+6qu+bevQeAovKXwq7fEDoFg2+rOy5cFZ/Q5+d9M43vPhnCgJkr6q61ttWBltBRl90iRaA9n/PR7rn7i8gcr6QMYJfvUhGj2RgXViTA2QCjR4+O22Nw1/YDWT/BASjm7vbpBLbVMPjW7bz1iVYeKEit4rzjJ6PdNEHSshJOTvesajoalMteWMHAp1cC8GnOQkY9cgWqtAyA3oVbIbcvtT0zY3IPhsiz/uJB/DjiueYzhsj1Xb/nD+tOiXi9rQJFR3avEywBv6BAM5+vgS3AUcA/oE75vwzwFb87kgQT48LQhsh9+SFy3k8k9Z366Knt0/jQEAxD755Jj+/r9xG/eO+WOPamhZiZTwNE5GhgilLK44H7wkD5/QofS6NhnYjMAh4FvgVeBL7wE9c7YjQTpyIkqmf2ZvVfteFbXmJ6xPo4q1RvQX07PhtJbBiMq7VSuP4Jxlw4A4CTH55B/zc28/PdRwBg/2Y1Ux+6DoDHb53FCcmKQ5dqC/XsZau5Iuszq5YU+ry/BVeZnvnYHYmQmkriDv0dLz/mLWyL6rcS/1cFF38RP5Xbjk4wS22ulGaztIh0WzL/7v+JPrn8E946T8+6b/v4PJIbx+psg5hlNxCRg9ETg3PRyyh1TkWVUj8FKtvsno9S6k4R+SPaXv9S4K9WMKIXWhLDIVj8xakwtB8Kht1ed7zzyCy6vKJnRfZePQEtLA0dg4mH3F13XPTdvXHsSQh0UOEjIoPRq1Hno7dg3gBEKXViKPUEZWSqlFIishXYio7j0AX4j4jMV0rdGlLPY8wX793Cof9P3+aSeyM3O3lkQb51ACUX3xY4cyti4SvawG707+HVux/jwm467lDlaUP57lntOfvk4/ZzT2cH2Zsr6sqtd3ax/gJOJ/Y+vSlcO6NJ/eNOyqQ2XX/fG8bb+Oq3j7Nmwos6cQI8dq8O7z7/90fDnoom5Q2xxzU8Nv+Hc9JL9d9JWtmgVulltxGvXMPA9THpQmSJkPBpzqBeRMRKPwUdyfQSpdS3gcqKSFe0UMgFSoBzlVJ7rLTbgcsBF3CtUqqoUXtzgAFKqQP9dPkX4H/Ab5RSxVaZkCNbNutDRUSuFZGlwCPAV8BBSqkr0FHvgnafbTAYDO0FUcF/AtZTb1BfAAwHzheR4Y2yFQB51mcq8FwQZacDHyul8oCPrXOs9MnACCAfeNaqx9Of3wLlzdz+WeiJyKci8ncRGUcIDkU9BDPz6Q78trFVq1LKLSIRcbNgaB30f/UvAAyMcz8MhjZBZLTdgjGonwS8bO21LxCRziLSCz2r8Vd2EnCCVf4l4DPgNuv660qpamCtiBRbffhGRNLRXlenAm/667BS6l3gXRFJA85Ah9POFpHngHeVUvOCufFg9nz+FCBthb+01kT2O1pN+K6rDuTeHtFW1GsbLHn+Rkb/HtznaiXGPXvqPX//+v9sdPlfAuW9O/HtrKaz6d+s9V/vx5/c3uC8aO0L9EipqTu/sav2KXfjO6v53dqTWfiTFnVZ34C9sw4PvvGcA1p2U4YW0Xd2AiVHltEvIXbe3wHsohdefrnoGaZPDN/OLNaEoHDQXUSWeJ3PtkxFwLdB/ZhG5f0Z3Qcqm62U2gKglNoiIj286lrgoy6A+4DH0Ut7zaKUqgBeBV61lvnOQc+wIiN8DIZR18/EaWlEJVRCz1lLWHvXaAAGvrAJgMI1kffGVGcAaxFq1FWDIaoEL3x2KqVG+0lr1qA+QJ5gygbVnqW1NkgpdYOI5AasQOQkpdQn1nF/pdRapdRu4G8isrOZ9uvoEMKncMcsAEb/fgZ/uu/7yFTaDjRdljx/IwW5embTeWB63Q7g3gFJpOxyRaSN/1d0Gak99RLyhYMX8doL4wEoG1WNI62WYX8qqcvr+T95GN/efYO1Ej6ZP51BDyfx80VPx60PK0qD9sTfOghiPydImjWoD5DHEaDsNhHpZc16egHbm6nrSOAwESlBy4UeIvKZUuoEH31+DDjUOn7b6xjgTutas5gw2hEgv/vUeHfBYDDEmsiE0a4zqBcRB1oZYE6jPHOAi0QzFii1ltQClZ2Ddo2G9fc9r+uTRSTJMuLPAxYppZ5TSvVWSuUCxwCr/AgeaDh7ajyT6liRTIOly/L9OInMiP6GsfMBKDxmQETqixeFJc2HXM9+djFbrjm87tjWrQsD/1pM4ZZn9DgnACXT6uMaFQy/gz7soHD5gxQMuJmyg3tSeqz/78/m0G5+lCsy/zODf/KeWceC8/TrYHRSdUzbfmL3QbgvstdfKIlp8y1GIhBMzp9BvYhMs9Jnoe0dTwGK0fsxlwYqa1X9EPCmiFwOrEfvx2DV/SZaKcEJXKWUCvUHpvwc+zr3S4cSPrFi4mHaYK5oaRsxlosAA994AIBT8pazcLtWGEh/XPt7S/xS/x5s/fVsvyCvVZuGGQwxxZdBvSV0PMcKuKpxOX9lreu7gHF+yjwAPBCgPyWAPxsfgAGWLZB4HWOd9w9QrgEdSviU9U/BFSHPQJdmrgJgbpdDm8nZ9mkQJO7xltdz/Qf/5YkJlna+3caX7wSO9lpU9WrLGzOEROH6Jxh1nR40/O+25mfDkWBhtdaum3fLcTAcPv+wjQ1K2sG+bwvxjmncWNMoaM2jDiV8DAaDISJETuGgzWHFdPOJ5Vw0KDqE8KlT2T1/bNTbKuh1FfTurts96n4A5n/9x6i3G2uOOUsPcHYeaEfZgSHx7Y/BEHM6qPCxPCKci7YPmquU+slyOHAHkAIcEkw9HUL4eMi//QtqicAuoTc+lvGOfuk7/jf1iMi208pZcOyzANgQ0H5BOaWX9iPXJdGzkb2P1VN7xqV/hubpOftbAD69rivHJAdtrtFi7ik+XR/csIuvxscn0F1YdFDhA7yAVtdeBDwlIuvQqtrTlVL/DbaSDiV8osnQu/U6edC7bW2UiWkXsXPKIUTCC3/BkOkUrnyo+YwhMt6uQ0HYMzMAmLvn+Yi3YejYCJHRdmujjAZGWi7WkoGdaAPVraFU0iGET8KAXACmdnmV6giNVk5erG17+m71+JrpDcCpn64gL2krX5ePjExDrZQXntAerTMsS7Hauu+1/gvef3Yp7q86U36l/i6qVREvn/tXAO6bcG6sumoIkqLKVwAdevu627sBMH/8E6RFKWDnnoooBRKKBR14zweoUUq5AZRSVSKyKlTBAx1E+BgMBkPE6bjCZ6iIeJxkCjDQOhfArZQaFUwlcRE+IvIo8BugBvgVuFQptddKCxhrokVUa8eWG5xJZNurWlzNt9V6v+Jvk88gZ6kO0ucGbCkpPHupVssfkLCPF/aM5fBXf+a+g94Nr9+tkOpjR/D23Y/UndcE+AG+f+jsBo43ahT0tlcC8OtF0XOn4jpB73e6gE8+nh61dtor3tFPjznLxr+eDEO/PgCp72fWn5wRlSaiS8cVPsN8XBO0q547gq0kXu515gMHKqVGAquA26H5WBMGg8HQWohEPJ+2iFJqneeDDix6FTpkw32EEH06LjOfRvEeFgBnW8d+Y03EuIsGL/JH3cWGU7oCkBXhukfeOJNlM4IPgpjf5feAUSIwtALaoWAJhpiG0Y4yl6E7D4FjTTRARKaigx7Rr1+/gA0UbngSgAEzD6Dw7Ja5/t/kyuBvZ1uqobgbuPcfbz+PV3YcBcD0nkW88+pxANwXeUWuuHPKjM+oVuHvQN933ms88OyUCPSoKWapLXJ8+fbN3PD9EqZ1+yKi9a6q7U7Wh79GtM6Yojq0tltswmi3FBH5SER+8vGZ5JXnTrRzO48flaDjUyilZiulRiulRmdlRXo8bjji4hlMX3Y205ed3XzmMBl7wQzGXjAjpDIFw25vPpMhbBrHVDJ4ERmv1m2RmIXRbhFKqZMDpYvIxcBpwDjLcR4EF9uixQx5aiPvTRzFxPSfQi7bzbaffcMyfaY5Tz6U67O1CnGFSiDn8+ZCoLct7nv/FRxe3sCrVPjbcEMdW3n62ue47Y/TQionaamQEduImx0dp9sekf+5h73uFG79x6Xwe1hxX8gD5lZDe9zPCYZIhdGOi8KBiOSj44mfrpTyDtnqM9ZEPPpoiD3Hn/pI85kMhtZCx535ADqMtlLqVaXUaeiJwvfoMNpBEa89n78CScB8EQFYoJSaFqFYE34pXDuD4z6+heMH/RJy2RplpzbF98wyecUWqpT+Kv+6dRzZT6zjlTFte0N8zHVLObPTUgDcSqiKwqOSLLVMuuMjAObeckJQZQo3PhXxfhj8k5B7ABd0e7vu+Q6VGmvGdOOKc+FN7fOw88sL6cMC5rveCFS0ddPOBUuoeMJoW5+giJe226AAaQFjTRgMBkO8ETruslukMGG0Da2KgqxpFGSFtg9kMMSDSNn5iEi+iKwUkWIRabJsZYXPfspKXyYihzZXVkS6ish8EVlt/e3ilXa7lX+liEz0uj5XRH4QkZ9FZFa0bSxbg6p1TNnzQW8qrnWEVOa9PYfyy9ShdKOMosV3N0l3btrCrVdeAUDqyh36YnHYXY0r7/80kjOPWxL1dsakanXbPQ+m8ePJ3aLeniF0KlRovxdvpi29oO541T+0l3P+EW6PWgkRmPlYL/hngPFohavFIjJHKbXcK1sBev87DxgDPAeMaabsdOBjpdRDllCaDtzWyJC/N/CRiAy2tjfOVUrtE70X8h906O3Xw79L35iZj8FgMLSEyFYXWZ4AABeeSURBVCgcHAEUK6XWKKVq0C/7SY3yTAJeVpoFQGcR6dVM2UnAS9bxS9Q7MKoz5FdKrUUPk48AUErts/IkAI6geh8GHW7m0/uLUnZcqVWmU23VQZXZXNmZzCe38tZRz/lMb9Mbp/7Yb6fCnRSz5o7NWMmi0Rc0n9EQU/YP6xHWc1C9JzmCvWlFhOY6p7uIeC8jzFZKzbaOc4ANXmkb0bMbb3zlyWmmbLZSaguAUmqLiPTwqsuvIb+IFKGFUSF69hM1OpzwMRgMhogQvPDZqZQa7SctGMN6f3mCNsoPtj2l1EQrRs+rwEloP5xRocMJn6Il9zD6cj3zueX2fwfMu6ZGe07YfaflvufTqHatVTHk+XL2nxy7mQ+Azdlx/ZW0VkrOgv0tnPlMX/hbhl6zrP7C7yPUqVZChNzrBGNY7y+PI0DZbSLSy5r19AK2B9ueFaNnDnqJLmrCx+z5GAwGQwuIkLbbYiBPRPqLiAOtDDCnUZ45wEWW1ttYoNRaUgtUdg5wsXV8MfCe1/Umhvwikm4JKUQkATgF7cMtanS4mY/BYDCETYSMTJVSThG5GigC7MCLlrH9NCt9FjpMwSlo5YD9wKWBylpVPwS8KSKXA+vRmmv4M+S3XOXMEZEkq65PgFnh36F/OqTw6f6O1mJ8ZsoJXNLXd7SGVVU9+eRx7amafrDopRtj1b1Wgfvbn3m8eDwA0wY09Wi8tlovSX6+fRCbl/RGJehfYt+DN3NR3wVN8jfH8socHEvbsJfjdoq4pEUKBw9+n0/eJT+ggHm1UdPWjS8R0gVTSn1Iozg4ltDxHCt0zJygylrXdwHj/JRpYsivlNoGHB5q38OhQwofg8FgCAfj4SB8OqTw8QQiG/DEDMpytCrogr0D2FuTwi/LtHLBoH9X0on9zP/qj3HrZzyxJSUzrOs2AMrcyby76WAA1m/rim1jCoMe+BGAVHZRXPpoXbnxtnN47I6zAHAdWsaVw4OLA/PddaNgJHz8WdBReA2xINFNmTt4delN1dqQ/oDn7LiOGcXHn7bf0BfiNtInHDqk8DG0T8ad9BcAPv6k/b7wYkV+96n64LGB8e1Ia8U4Fg2bDi18Bv/pR+Y+pUf07i3bgErWVuqop1wdv361BpTLxeYb+gOwpscwUuf+oI8rrUiwt/ouN9/9FhNTLgTA1r0rf598KmdepGc/qfamRr1/W3I8AHlfRN+VjyF0Jhy4nP3uwO51Pto2FADH5TaorNTHbKRwyzNR7188Mctu4dGhhY/BYDC0GCN8wsIIH0O7oyC3PjpmYcnMOPbE0J4xM5/wiKvwEZGbgUeBLKXUTuva7cDlgAu4VilVFK325+5rL+51I8+8mtdaXLao8pUG57mv6OW2s0d+R61b2zWn2Gt58+dD6TVXP4Llk8fyzWs3tbhNAMeWfZQN60bLfTAbPGy6eBgAhzn+R7kzsMJBpkP/fyv36L8ehZ52jxE+YRE34SMifdGuwNd7XQvk7ttgMBhaBypi7nU6LPGc+cxEb1u/53Wtzt03sFZEPO6+fVuCGtoEJRfqGFf5I/+I7NejY3eXdPJWrASgqPwlv2VDYffhWeResYrtd/ePSH0dmX3DnPqvM8VvnnkfjKbXVzWkFO8EQLqlUlj8qN/87Qlj5xM+cRE+InI6sEkp9YOOW1RHQHffBoMvjpzyuD5I8OWw12CIEspIn3CImvARkY+Anj6S7gTuACb4Kubjms//sIhMBaYC9OvXr4W9NMSSucvup2CYtsGp6JeGI3NYROvv/GMpW/4yEJLhizm3RLTujoY49U+x3NVwB219eRe2ztNOkQe8vAaAwk1Px7ZzrQQz8wmPqAkfpdTJvq6LyEFAf8Az6+kDfCsiRxCce3FP/bOB2QCjR482j0EHY4JjCgClk0cb1+wRJr/rH+DBIfHuRuvGGJmGTcyX3ZRSPwKeqHqISAkwWim104oh8ZqIzEArHOQBi2LdR0P0cK7So+VNV3an39zIiI1tR+hReuY7xaSsgKKqVyNSb0cmM0dHVC6vbajptmZTFoyoYu0Ft8Nf4tGz1oNROAiPVmXn48/dd5y7ZTAYDE0wwic84i58lFK5jc6buPs2GLwZf9T9PjcHDYaYoTAKB2ESd+Fj6FjMd71Rf3JNy+tZ9cJIABzrhLzp9X7hzJJb+Ej3ruzbnQZAeadSAEp2dQUg591EnemCuHStVWEUDsLD7NUaDAZDS1BBfppBRPJFZKWIFIvIdB/pIiJPWenLROTQ5sqKSFcRmS8iq62/XbzSbrfyrxSRida1VBH5QER+EZGfReShFn4rQWNmPoY2R8LOMnr/n9ZZyZy/HFJTmbv773HuVftBJSXQq9ceACpqtaq1a0UmABlfF8etX62JSBmZiogdeAbt7WUjsFhE5iillntlK0ArX+UBY4DngDHNlJ0OfKyUesgSStOB2/x5kbHaeUwp9amIOICPRaRAKVUY/l36xgifVkb+iDsBUL+WAGYZyWBolSgVqWByRwDFSqk1ACLyOtrTi7fwmQS8bIXTXiAinUWkF5AboOwk/n97dx8sVX3fcfz9vcu9l4sIF0QUUQJGtFWiBlEypipNUC60HdI/TBiniWNNGFNtop1OJOM0JR0zg7aamkkmDElozYMxTU0r6QhoMrGdTn1AE0WpCkRsRFEkBB+A+7jf/nF+994F9uFcdvf89u79vGZ29jzufvd3z93vnnN+53tgUVj/XuBR4FZKVJFx98eAXyQfzXvN7Jckl7rUjZKPjDobt98RO4SmtuvqaUwaeAuA/oHkyPxAe/JF6z290eJqOLU55zMTeLVgfDfJ3k2lZWZWWPcUd98D4O57zGzw8paKVWTMrBP4E+CekX6YkVDyERE5DiM47DbNzArvlrguXCQP6aq6lFomdUWYFK+VzDQbB/wQ+NrgHlW9KPmIiIyUA+kPu+1z9wUl5qWp6lJqmbYy675pZjPCXs8MYG/K91sH7HD3fyz/kaqn5NNgNm3TJU4SV64bDvUO13Tr2zaZuXc8NzS+6e31McJqPLU57LYFmGtmc4DXSDoDXHPUMhuAm8I5nYXA2yGpvFVm3Q3AtcCa8PxgwfSiVWTM7HZgMvDpmnyyCpR8RKQqXVM/MzQ8lnod1qK3m7v3m9lNwGYgB6wPlV5uCPPXAg8By4CdwCHgunLrhpdeA/yLmV1Pcs+0q8M6RavImNnpJEWfXySptQnwdXev250BlXxE5Ain//0TbL8nOUrUNv0w/R3O9tXzePmWv4ocWWOpUW833P0hkgRTOG1twbADN6ZdN0z/LfDREuscU0XG3XdT/HxQ3Sj5iEgql378LiY/uRuAja98NXI0kamqddWUfETkCA/33c+SBasB2PWxyfjsbmZsaCu5/Fg61DYouchU2acaSj4iIsdDVa2rouQjInIctOdTHSUfETnG5qdWA9A1bSXW0QHAxlfresH76KJzPlWLVtXazP4yVFXdZmZ3Fkw/puKqiEhjSWq7pXlIcVH2fMzsD0kK3J3v7j2DdYdKVVzV3UxF4ti0b13lhcYqHXarSqw9n88Ca0JlVdx9sPTDUMVVd99FclHVJZFiFBEpzpPbaKd5SHGxks/ZwGVm9oSZ/aeZXRyml6reKiLSWNzTPaSouh12M7OfAacWmXVbeN8pwIeAi0nKQJzJCKq0mtlKYCXArFmzahGyiEh6yitVqVvycffFpeaZ2WeBn4SyEU+aWR6YRroKr4Ovv46kAisLFizQZiAimbK8jqlVI9Zht38HPgIQbuHaBuwjqbi6wszaQ6XWoYqrIiINw0kuMk3zkKJiXeezHlhvZs8DvcC1YS+oaMXVSDGKSBlLT76BgTmnkXstuespYU9g455vRIwqG4brItMqRUk+7t4L/FmJecdUXBURaThKPlVRhYM6WXLCp5KB35szNM36kp24TVtvr9v7XtW6gpYJE4bG+z54Fv0nJH/m8W8cTP06LW/sx99Lv/wgm9CBHzp8xLSBgtd5ZOBHI35NaRxLZ93M/iuSDj5TI8cSnZJPVZR8ROS4HFg8lxMfeJoB4OHe+2KHk63Bcz5y3KKV1xERGc0sn0/1qPg6Zl2hnNhOM1tVZL6Z2dfC/K1mNr/SumY21cweMbMd4XlKwbyiJczM7Ctm9qqZvVdVw6SkPZ+UrmpdAVY8V+dmzYTWpCm7Z3XS39HChKJLJhZ13XHMtLb9PeTeHT5cNbBzF5DcW2XQko5PAtAy5wwOz5oMQN/EHD2dw3E1+qGQBZ++G4DO7YdpeXxbyeXG3C/pBnZly9W0XDQP0K/VYbW5gNTMcsA3gCtJLjXZYmYb3P1/CxZbStLzdy6wEPgmsLDCuquAn7v7mpCUVgG3Vihh9lPg68COqj9YCko+dbL/ExcC0D214LrZPJyw1xm/r69u77vv+oX0TSwx8/c7R/BKI1m2vLZ3kmcVWRz93lz2Pqbft3VofNNY/ZHg1OqczyXATnd/GcDM7icpM1aYfJYD3w09gh83s04zmwHMLrPucmBRWP9e4FHgVgpKmAG7zGywhNlj7v54eJ1afK6KmiL57HjuVRZfXr8Ocvn2XLSGuvQTdw0NnxgphhgWX/4VWvrCIYunkj2kwr1AqZ+uC79E/6R2AA6fMr7sXvzSk29g41trswms0dTmnE+xkmILUywzs8K6p7j7HgB33zNYvDms83iR18pcUySfLLy94hJ6Okf+iyDXc+R49xSje0qxWxK3UZhexs+fRkv/kUvs/fPkUK9XOPZx9HvGNtA+ODTcfvsumAAXXHzEcoNFGDv25Zmwpzub4KSi/KL59LXnkuFW+N3yDwAwddP2mGFFN4LrfKaZ2VMF4+tChRZIV1Ks1DKpy5GleK3MKfnIqDF4vqj9QJ4TfnMIgJadyQ+/Tfu/FS2usW5JxyfZfPh7scMAkj0xgMMXzcFb63yGKn3y2efuC0rMS1NSrNQybWXWfdPMZoS9nhnA4J0DUpcwq7emSD4DE9v43TnlDg7URi7DH+N9E479x2npze79Y+qZ1ELPpOG/p52nu2pk6eCcSRyalhsa7+4cPzSc64b+8cmP5/1Lz6bzx7/KPL6G4A4DNTnutgWYG8qJvUbSGeCao5bZANwUzuksBN4OSeWtMutuAK4F1oTnBwum32dmd5N0OIhWwqwpko+ISOZq0OHA3fvN7CZgM5AD1rv7NjO7IcxfCzwELCO5v9kh4Lpy64aXXkNyt4Drgd8AV4d1SpYwC3eUvgaYYGa7gW+7++qqP2QJSj4ikSydfQsAvXOmV1jySK1PD59r2fTOP9F13m0A9E0v1c1xhDpba/M6ddI1bSUANn0aHnpmvTMvuchgcpaB1KjCgbs/RJJgCqetLRh24Ma064bpvwU+WmKdoiXM3P0LwBdGEns1miL5ODCuW914x4Jxh/Pkep3ce8kx0P4DByJH1HwGWi3V/1Pni+/CvLPYvOVvM4gqnQNXzgUgn6tzd2EHdOlAVZoi+YjENHTx76SJMKFjeMb4dvpPSvZGrH8AzOiZOtT1r2wX5rSu+KM76ai82Ki2dNbNyUBbK4fOOblBPq+Dq75ONZoi+Vg++UUszWvSY/8HQP8bSaedjU1QoPTgBacByZ7GSPRfcW4yEH55Hzw7qZySz5VaY2Qs74w7XPlXve16rTZvOEKHPpzs3RS2mw14eE7GW+rde9ipVYeDMaspko9IPXRNug4AO7FG51LGsKtaV6RaznI5Wk6aCh1JD7v8lImNW9JHVa2r0hTJx/LOuEP6FdLUWsJXUAMe6ui79DwABtrLfU2Gk/gF31fWl4yM66vuS6ylJ2mTzL+k+/orL1MDh+adNjzitWu3qin5VCVK8jGzC4G1wHiS7n5/4e5PhnlfBK4HBoDPufvmGDFKY1uyYHXd3yObClejXxZ/i8ZTm8KiY1msPZ87gS+7+0YzWxbGF1WouFpScs5Hd9tuZr1nnQpAbuZJWH82ez8DH3g/APm28idTBgumjrVtMH/O7FS3DADIX3o+3nJ86bwh29UZum24HJ9Yh1MdmBSGJzNc3mGo4qq77yK5qEqXt4tI43FP95CiYu353AxsNrN/IEmAl4bpqSuumtlKYCVAe3vtyv+LiFRWs/I6Y1bdko+Z/Qw4tcis20iuvL3F3R8ws48D3wEWM4KKq6Eq7DqASRNneu5w/e6RI43Fx2Wzw27hyyWnbvzFtYC3pPtbtPRk0zkhMw7egJ1fRpO6JR93X1xqnpl9F/h8GP0x8O0w3DAVV0VEylKFg6rEOuz2OnAFyd31PsLwbVuPq+KqDeRpOdhgN7ERkeam8zlViZV8PgPcY2bjgG7CuZtyFVdFRBqGu3q7VSlK8nH3/wYuKjGvaMXVsvKOac9HRLKkPZ+qNEWFAxGRbDk+oIMy1VDyEREZKd1SoWrNkXzyeTh4MHYUIjKWqKt1VZoj+YiIZMgB155PVZoj+eQHyL+nPR8RyYjrZnLVao7kIyKSMXU4qI55E3QXNLO3gIPAvtixHGUaiimtRoxLMaUz2mJ6n7ufXM2Lm9mm8B5p7HP3rmrerxk1RfIBMLOn3H1B7DgKKab0GjEuxZSOYpLj0bB3qBURkeal5CMiIplrpuSzLnYARSim9BoxLsWUjmKSEWuacz4iIjJ6NNOej4iIjBKjPvmYWZeZvWRmO81sVeRYXjGz58zsGTN7KkybamaPmNmO8DylzjGsN7O9ZvZ8wbSSMZjZF0PbvWRmSzKMabWZvRba6hkzW5ZxTGeY2S/M7AUz22Zmnw/To7VVmZiitZWZjTezJ83s2RDTl8P02NtUqbiiblcyAu4+ah9ADvg1cCbQBjwLnBsxnleAaUdNuxNYFYZXAXfUOYbLgfnA85ViAM4NbdYOzAltmcsoptXAXxdZNquYZgDzw/CJwPbw3tHaqkxM0dqK5Nb2E8NwK/AE8KEG2KZKxRV1u9Ij/WO07/lcAux095fdvRe4H1geOaajLQfuDcP3Ah+r55u5+38B+1PGsBy439173H0XsJOkTbOIqZSsYtrj7r8Mw+8CLwAzidhWZWIqJYuY3N3fC6Ot4eHE36ZKxVVKJnFJeqM9+cwEXi0Y3035f9Z6c+BhM3vazFaGaae4+x5IvlyA6RHiKhVD7Pa7ycy2hsNyg4dtMo/JzGYDHyT59dwQbXVUTBCxrcwsZ2bPAHuBR9y9IdqpRFzQINuVlDfak48VmRaz+96H3X0+sBS40cwujxhLGjHb75vA+4ELgT3AXTFiMrOJwAPAze7+TrlFi0yrS1xFYoraVu4+4O4XAqcDl5jZvDKLZ9ZOJeJqiO1KKhvtyWc3cEbB+OnA65Fiwd1fD897gX8j2a1/08xmAITnvRFCKxVDtPZz9zfDl0ce+BbDh0Ayi8nMWkm+5H/g7j8Jk6O2VbGYGqGtQhwHgEeBLhpomyqMq1HaSiob7clnCzDXzOaYWRuwAtgQIxAzO8HMThwcBq4Cng/xXBsWuxZ4MEJ4pWLYAKwws3YzmwPMBZ7MIqDBL67gT0naKrOYzMyA7wAvuPvdBbOitVWpmGK2lZmdbGadYbgDWAy8SORtqlRcsbcrGYHYPR6qfQDLSHoF/Rq4LWIcZ5L0pnkW2DYYC3AS8HNgR3ieWuc4fkhyuKGP5Nfe9eViAG4LbfcSsDTDmL4HPAdsJflimJFxTH9ActhlK/BMeCyL2VZlYorWVsD5wK/Cez8PfKnSdp3R369UXFG3Kz3SP1ThQEREMjfaD7uJiMgopOQjIiKZU/IREZHMKfmIiEjmlHxERCRzSj7SlMxstpldUzC+yMz+I2ZMIjJMyUea1WzgmkoLiUgcSj4y6pnZxaGQ5PhQaWIb8K/AZeGeLrfEjlFEjjQudgAi1XL3LWa2Abgd6AC+DzxGcl+XP4bksFu8CEXkaEo+0iz+jqTWXzfwOeCyuOGISDlKPtIspgITSW4qNj5yLCJSgc75SLNYB/wN8APgDuBdkltRi0gD0p6PjHpm9img393vM7Mc8D/AT4F+M3sW+GeSCsgi0iBU1VpERDKnw24iIpI5JR8REcmcko+IiGROyUdERDKn5CMiIplT8hERkcwp+YiISOaUfEREJHP/D8yVkjtuzc6lAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ds_grid_spec.AREA_LND.where(ds_grid_spec.AREA_LND > 0.0).plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x7f1c4a2adad0>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAELCAYAAADawD2zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnX28HGV1+L8nNxEC4UUSIJCQhJcQwUsQGomFVFAJRLBiiCCWtL5Q0ypViKCC0Ma0pEUForaVmvpuQAVChB8icKFFG5FgeAsXQwxo0LyBRAMEY0huzu+PZ2b32bmzu7O7szuzu+d7P/u5M8/MPHN2dvfMmfOc5xxRVQzDMIzOZ0jWAhiGYRitwRS+YRhGl2AK3zAMo0swhW8YhtElmMI3DMPoEkzhG4ZhdAmm8A3DMDJGRC4SkX4ReVJELm7WeUzhG4ZhZIiI9AIfAk4AjgXeISITm3EuU/iGYRjZchTwoKr+UVV3Aj8GZjbjREOb0WkWjBo1SidMmJC1GIZhtAEPP/zwC6q6fyN9nP6WPXXz7weSnW/l9ieBP3lNi1R1UbDcDywQkZHANuAMYEUjspWjYxT+hAkTWLGiKdfIMIwOQ0SebbSPzb8f4KG7xyXat+egNX9S1Slx21R1lYh8FugDtgKPAzsblS8Oc+kYhmHUgQK7Ev5V7Uv1a6p6vKq+Gfg9sKYZMneMhW8YhtFKFGWHJnPpVENEDlDV50VkHHA28OepdBzBFL5hGEadJLHeE7Ik8OHvAC5U1T+k1bGPKXzDMIw6UJSBlNLLq+pfpNJRFUzhG4Zh1Mku2queiCn8DmLXpiMLy0NG/zJDSYy84H8nwH0vom31UOn7lUb/jdKK778CA6bwDcMwugOz8GtEROYCf4u7YT4BfADYA/g+MAFYC5ybxiBGvZZHUmvhmEsWlqw//onrATjjtHO5856bUj1XNZK8V3sK6FxmTL4SgDvvKW1Py/rOgxWfNQrsaLMSsZnG4YvIGOBjwBRV7QV6gPOAy4D7VHUicF+wblRhyOhfmhI3ynLGaedWXK9G1t+vM047t/DKA4oykPCVFzK38HEyDBeRHTjLfgNwOXBKsP1bwP3Apxo5iVkkg5k+5JxBbX27bs5AEqMeQiu+VdhvKILCQH50eSIyVfiqul5ErgF+g8shcY+q3iMiB6rqxmCfjSJyQNzxIjIHmAMwblyyKc4hoZVQzdVSi0XzxLVzIy1u/a6VMGMyg87XLGupmgvp9IOPbcp5jfwRZw37bUl/B7W4JZtFq92i1XAzbduLrF06rwXOAg4FDgb2FJHZSY9X1UWqOkVVp+y/f0N5kAzDMGpEGEj4ygtZu3ROBX6tqr8DEJFbgROB50TkoMC6Pwh4Pkshm0XcI3Jocd218qpWiwPEuwnSlsU/R1bvM29MH3IOd294PDf+aaM6Cuwyl05N/AZ4k4jsgXPpvA2XFvQV4H3A1cH/2xo9kR9/XO3xtBmPhEXFdlXTfaF3rbyq4EICGOhfTU/vpMJ6T69ri+LvY3QOWyfuW1j2byh33nNT1RtM+LvJ2p1TjiwNJAVebbP8k1n78JeLyC3AI7h0oI8Ci4ARwE0icgHupjB4dLGN8MM1x/RtHhQqlwatHsBLSjW5otvN4jfaiV2aH3dNErK28FHVecC8SPN2nLWfKqHlftdKgOwUS6UnjGJ7uvKFFn1oxZs1nz/cYHrp55QGw5cuB2DbzKkl1j7ED4RGn3DjnngtYiecaWsK3wiITsSqhTxZ7LXIYhZ6bRy1dD4AYzOWw8dSdCRDEQbMpWMk4YzTzo1Vjs1Q9Fla8248IT83r3YiOvZS67FRQks/5PSlpeG5d294PHH/5W4E/tyOWvpLQnw48epM546YS8cwDKMLUIRXtSdrMWrCFH5A6H4ZvfABNs09cdD2wZOq0qcTLOFG30Ol4+txF82YfGXB2q3VEjRX1uDZ2DYTu4ibeGUuHSOGMX2bAVg/fSSQ7g1kztIfsmjmmcxZ+kMAFs08M7W+jdbQ0ztpkBsm2laLe8ffN869E0faM7DD/tJ27eQp4MAGbduAuBwyxFj1PuUGYOMU9/Qh5/DhNU8Df1O2v1ufOQ6A6yceUXiiGFNRAqNZT0Bp9Bv7nQrYNnMqAJtnvwLA+Hk7ihvnD2v43K3Ef5/dbu2rCgNqFn5bE+fOqZf5H/l2ceUj7t+8L/9NaXsKhBZ9sy37rRP3ZcSaLU09RzNotkWY1IKOw1f+cXI2KnsjstVD366bK978GiF8UjjjNHdN8uBG25WShR+XJl5V/5RK5x5dpfALfvoU+wwtdXDKPO3+jfSYMfnKgvvLaJxut/ZdHH7jFr6XJv5oVd0mIjfh0sR/s+HOI3SVwi/HprknMnrhAwCBK6ZIqMSjzP/It8tu2zT3RK6fWHxa2DntJRZM/kH8cXMpnJuINddIWF4thJNxqlnv7WjdJyG0FKOuHf/6V/Klx31G/v5hOOT4NbV/lvVY6HHjAWly94bHW5ZxNS78002czB5F2KGpqdC4NPGpYwrfMAyjTgZSiMMvlya+4Y5j6BqFP33IOVUHZo3upBPCYbsBP/lhuSeuekNw66HGmbajRGSFt75IVRfBoDTxW4CbRWS2qi5OVWBAtM1qMpZjypQpumLFirLbowNJSQZnw8HVcq6baoShmNXw89w08ijeDPfPsxWiSEqiTXJMGK4KrQtZTfI5+p/3uiW9g7aPndWfulzlzt8ovoKN/tbqDcuMunPSVPgi8rCqTqlLsIDDj9lT/3Xp0Yn2PW/iirLnE5FzgBmqekGw/jfAm1T1I43IF0fHWvjNihQwDKN9CJ/emhHRk9agLeXTxKdOxyr8NEhi2UcnVPltSXLOV7OwooPI1088oqpMjRJn1Sex5sPjWmX5Jx1sXjTzzMIPfhHNcd+Uy0Za6fP1t42d1T/Iyl+3pLcpVn5pbYTavo9xRK3pcD0s6pKUaonaituvLJHbz+vfyklZiqTlwy+XJj51TOEbhtE1+GHUjaJKalE6ZdLEp07mCl9E9gW+CvTinpI+iEsK/n1gArAWOFdV/1DvOeL89WEoZL0TraL++dELH0gUrleJgf7VJRa9b82HvueewAhspl86tNB9Sz+JL7/VPv16wkTTyFAahmuG1rBfmnCgf3VDpQqHLtsbgH3WDvDihJ7C93OftQODsl3WQyXfdq3XoZqfvG/XzezadGTVFAu1pmAe/BleVTZUs9L3tnEktYlXrSIP84K/CNylqq8DjgVWAZcB96nqROC+YN0wDCM3KDCgQxK98kKmkojI3sCbga8BqOqrqroFF6L0rWC3bwHvykZCwzDKccwlCxsq8tMJDDAk0SsvZO3SOQz4HfANETkWeBi4CDhQVTcCqOpGETkg7mARmQPMARg3blzZk0TdN1H3y5i+zSWl35K4Cp6dP6xqHpQkDD7u6Yj75unoIYBz44QuhQ+vaU66gCzDLkOX1dmHPwokczekHYlRaaZz6LqB0jDBsM0vHVhuVqrvFgrPs3PaSwBsngYjF+9Z2LcRd06aMel9u26uSckPGf1L+na55dC947t2ml1Ra9XMUrf4FStD2/GxhvtWpO0KoGR96xkKHA9cr6rHAa9Qg/tGVRep6hRVnbL//vs3S0YjQ6bNuoZps67JWoy2Js2JZXHKPs9W/lFL53PrM8dx6zPHpR6qrbhB2ySvvJC1JOuAdaoami+34BT+cyJyUGDdHwQ8X2vHfn4ccBaUvx5aVmNwllWtA4Bx1m+lySBJngC2Ttx3kKUeXfdz3od9LpqZn/zgaXLdJ88vWQ+fwpYtuTTV85QbyC33mYWW+ekHH0vfrpsHDRhGC3yXy00/0L+65Klz/LxJkdBW7ztZQ377bTOn1jWgvX76SEZXiAJtNJOsb+1nRTgong5i+fBrQVU3ichvRWSSqq7GTTj4RfB6H3B18P+2DMU0uohOTbPgv6/wxvnxz90AFN1mjXDMJQtbUhUuTyiwK0cDsknI2sIH+Chwg4i8BvgV8AGcq+kmEbkANwut5mexJ66dy1GBPzScvBI3GSZqfUd985WIq1BUqT0Of+wgXA+ts+g28KzeicW2atZc3vLYJ/G1z5gMIzwLd6B/NUyc2kyxBlHpc+zpnVT2ffh+6Thr3/fbF8J7g/bx89w5Kn1mlVIihMVW6uGJa+cy3XsKjlr0fkbZetONZMm2mVML1zut6Wxm4deIqj4GxOWYeFurZTEMw0iKqpiFnxVPPfMc02Zdk7p/1+hs/FwrnerOqUY4+/S2zcc31E/eCqLMmHxl00tI5inGPgkdo/CrUS3ELiR8pE5yXK3tPnHummrbhi9dXnhkH7FmS2G/uJwyvksgT+6cpITvyQ9HHL50eUsVSfg5NhLuOWT0LwuDuuENxXf9Rd1E4XdxeD+DCuJEZQvzN+2c9lLBDRl+R+K+Q+G1/PjS+FBfKA1dLrdt3pfjXT15ZPy8HSzCBT74LsI0cAVQelLts9m01+0pAeFkEMuWmV+61ZI2Ogs3aCuJXnmhYyz8nu07eXFC6d3WH6z1B2OjA2fRQdxWU84CHzRw1zupaPF61l/c8WlZ9UkzUparFdtIDdnoZKNmW/fNLooduo1Cy3wMRSs9HEz0Lexq4byFAd++YltP76SSpz+obdLW8e9xv5kNC+MnLIbrebbqQxoJyEh8jjazmTtG4RvtRTPzlBvZE07G6uRQzXacadsxFa+GH3GwnrzHB8tur3R3j8tlnwb+JJroFPqQqBWyfvrIijLE+WZb5aNvRnhnJYUfuuXyMADYDHzXVtTSD9vq/T7+8m9fC8DhFz9YcT8/jHPDyUXlVe24KNHvcEiWCj/8/oSyrZ8+spC64qmzP9NwxavRr99PZ984PdG+177hpobPlwZm4RuGYdSBajpFzFtJezmgDMNoKzp5gF4Rdu7qSfTKCx1j4b9mg/LsDcWY25GL9yxxP1QarNk57aVCZsJGHqP9Ygvj5+0oDBo/u6SXsbNce1wonk80c2eULEMso+eOvt+06VRXjk/o6hvdH/9dqJcjvxrUC6qSg8cf0D18aeU+owOgIatmzitR7Gm6RRvBL7UIQd6sYID7qZTO0VEzbUXkpSrHC7BRVY+ssp9hGB3KMZcsZEzWQmRAGJbZTlSz8J8J0haXRUQaz7yUAq8eLAUrc/30kbw4AV6cUBw88ien+FZ8aKmE1uvm2cNYz+ABtCT4Vq7r19UjHLps75LBsegEKZ8Ra7aUteIH+ldXnIzTSiqVjis3QO7PZrXonNJrkGfXR1zuHv+3BPl+L6Gln75c6aRWEJFJuJKuIYcB/6SqX2i48wjVFP6sBH0k2ccwOGrpfMYHy6b4jU4gjZq2QabgNwCISA+wHqjiYKuPigpfVX9VrYMk+7SaMX2bB1mZz04b5oVJFi33ocv2LrHk3TGDLftac62UWrilfvmoVR/6UcMp8b5f1U+lEFchKStGLt6TzbNfKaz71zvOum9uMWkjScW2ctlcw21+OGW0doT/3auWFz+vN/G7Vl7FUUvnu5WzG+9PFXakPyD7Npxn5dm0O4aEg7YicjbwWeAAnN9eAFXVNKsJ5Jq8PaYahlGZ2N9sioZHjROvRonICm99kaouitnvPOC7DQtXhqRROp8D/lJVVzVLkEZ5zQaFI0rbfIsFiuPPvjU6pm/zoIlXUWu00QiUSlE30eRgUGq9R621JJb9QP/qkjGD4UuXp/5EsHn2K4OuS2m1piLRpyOz9tOnJH1Clc867rswKDotJrqnUg2AdiH8bqYVpVODS+eFahOvgpog7wQub1SuciQdcXiumcpeRHpE5FERuSNY309E+kRkTfD/tc06t2EYRj00IXna24FHVPW5ZslcLSwz9HStEJHvAz8AtofbVfXWlOS4CFgFhC6iy4D7VPVqEbksWP9USucyDMNIhZQLoLyXJrpzoLpL5y+95T8Cp3nrCjSs8EVkLHAmsAD4eNB8FnBKsPwt4H6qKPxXDy69i66aOY8Z87xcJRTzaEDUreC1Tys/maiWiUZ+3plotsJtM6fGZjCs1e1SqQSfP6DK7N6CvM/OH1aYEFbPOcu5bcK2u1ZeBTMHH9furoC806jLLi4EOZpFthMmwhV+wykM2pJi6mMR2QOYDvxdKh2WoVqUzgcCYU5S1Z/620TkpJRk+ALwSWAvr+1AVd0YyLBRRA5I6VxGBsT571fNnJeVOIaRCgrsTMnCV9U/AiOr7tggSQdt/x2I1j+La6sJEXkH8LyqPiwip9Rx/BxgDsDQ/feBPcrvO6Zvc4n1PnSZ8x6FaRVCizhpSOGz84cxcvGeJQOtoUUUWvBbg4HTzfN7S6zqcvnJq+U/D/epRk/vpLIZOsfP21FxUC+aETP6VNOMFApG/mlH636gfzXrlvQ27TvbcTNtReTPgROB/UXk496mvYE0AlBPAt4pImcAuwN7i8hi4DkROSiw7g8Cno87OAhrWgQuPbK/zcIo84E/2cowOo2OUvjAa4ARwX6+y+Ul4N2NnlxVLycIQQos/EtVdbaIfB54H3B18P+2an3FhWX6uNDIol/b9+eXm0AUEq771vW2mVPdcbN7S/YZtB03SckPk/SppRoR1OarDfetJY99uf2qPflUs6KiybaMZMTlyW8V7R6GuW5Jb8l62uHA7VgApZoP/8cisgw4RlXnt0gmcIr+JhG5APgNYAVqDcPIHWmkVmglVX34qjogIvs1WxBVvR8XjYOqbsZNMTbajGmzrimNEDKMTkU7z6UT8qiI3A7cjOcXSTEOPxUqlSp8cUIPQ5ftnazkWhBWGI4DFB8Fewt586NKzc+OOXaWc9OMDdIflXPnAE0dVArx3TSNFBUPB8HKETej2F2v7HL4twtxuZpKZ4o3Tq0lKtvZnQNF1+G0xdcAMHJx8beQxkxbBXbuaq8aUkkV/n64jGJv9dpSicNvNTMmX1n4IYXKf8bkK83HbBgdRuEGWiG1SSN0nA8/JIzHzzO6bXtphr/I9riKQtWIs3BCa2H8vB1lQxjXLektCdn0ibaNXzO4qLlPKHOlIuy1UIt1P9C/mpGLpwbnf4XwCSYO996LTz2hrFsnOsty1cxL65K3m4ha+dGC5nHtEB/OGz5tJbXok4QEtyvNrBKnnajwg9mw/44Lo1RgGXCRqq5romxNp1zo5jGXLCxZDyN6xs7qhwruGcMwuouOG7QN+AZwI8VomdlB2/RmCNVsohaSbz0P9K+umO+74Ltf7CwoP2zOT1+wefYrjMTdHIYvXV7ix/dTK2yduG+hT39yli8b+BbzvoU+07DIfFdWtA6wG18IKoExOKTNd30V8owHtLv/Nyv863bMJQtLvqv+sv8k6H8PKln2cW2dmD6hEkkmLiZF23DQNumIw/6q+g1V3Rm8vgns30S5DMMwco4wsGtIoldeSGrhvyAisylmcnsvcWWhDMMwMqaVs+w70ocPfBD4D2Ahzof/QNCWK6qVcCuHn3cGSsu7+e6dTXNPZOysB0raw33Xe/uNnfUAPb2T2DrRrUezY0aXx5apXrluSW8h58/66e6RvvBYXqZARa2EYaYhfmjlyMUU3sPIxTEHe1kxLbKpOZSbaVvpsx6xZkvhc9xwsnDNmcUPb9HMM3NTJrPVpOnOgQ7MpROiqr/BVWIxDMPILS3NoaXOj99OJI3S2R/4EDDBP0ZVc2flQ9FyCa2j0QsfKDvIFYajRdt8Qivbt/xHL3yATXNPLHkCCLeHbfusLQaHRksOVpPdP281klhq5QZ/fWvQyBdPXDuX6UOCOImkWVQnFoveFweALynsd/bK0uMK/XcwaVv2Pp0apXMb8H/AvQwOcTcMw8iULLLjKp3rw99DVXNdYnDHgXuWWPTgJlvFEbX2y1n/YV/lwjSj26L7hZb8prknlsTyR1Mt+BZ/dIKV77vdOnHfkn2rWfblsgP6Fr1Z9/kmDJNMqtDsia25Fn0pwsCu9lL4SeOF7ghy1huGYRgBqpLolReSKvyLcEp/m4i8JCIvi8hLVY8yDCO3xD01TJt1TQaSNEZW4xCq6Sl8EdlXRG4RkadEZFVQfCp1kkbp7FVpu4i8XlWfTEek+hj23Cslg6pJiT7+1fo4GLpcKmU2HL3wARfQGhAdtPVnN/ozVndOewn63HL4mF5u8Dc627Jc/hBXmrG4vnn2K4NCM8tRy75Gc/HdNiPWbCkpdlML9R6XN1rnxiklxbDMLwJ3qeq7ReQ1VCzYWj9JffjV+A4N1rc1DKO5TB9yTsFgWLakfDK70MqvtI/hSCMsU0T2Bt4MvN/1qa8Crzbe82DSUvh13eZE5BDg28BoYBewSFW/GBRc+T4uDHQtcK6q/qFaf5UmXjVCkieHMX2bS6z86DHRiVrlWDVzHjPmpRNx4Of28Ynm8vefBuIsvfCHf9TS+VbcJCOiuYmmzbqm8Ln5mS7TyrFkVEcRdiVPmzBKRFZ464uCmtwAhwG/A74hIscCD+OSU6b+Y0tL4dd7n9sJXKKqj4jIXsDDItKHu9Pdp6pXi8hlwGVArqOEjM4m9BO3U3Kx0FKvNz1wnD/frP5SalB8L6jqlDLbhuI8JB9V1eUi8kWczvvHhgWMOVFmqOpGYGOw/LKIrALGAGcBpwS7fQtX+rCqwm+1H88P6YybgOUTnZQFlK2+VU+x5VCOWsP4oLJV72PpE/LDsiWXFj5j36Kv5JOPhgPHPTWEWHhnAjS1OPx1wDpVDQflbsEp/NRJK41bw/4mEZkAHAcsBw4MbgbhTeGAMsfMEZEVIrJiB9sbFcEwDKM2NOGrUheqm4Dfikh4534b8ItmiJs0tcJJwGOq+kqQNfN44Iuq+mwg8JsaEUJERgBLgItV9SWRZHfNwAe2CGBv2a/Nslq0H9Gc92bxdx/+U4C5d1KdaftR4IYgQudXQFOqDCZ16VwPHBsMKHwS+BpusPXkRgUQkWE4ZX+DVxT9ORE5SFU3ishBwPONnidtkszI9feppSC1r0jLDeCOWLNlUI4L378cVxTbx8/Ts3l2cUqFhV3Gk0fffeiS8T/nWksb+kQVeLvE5G+ae2JdIdlpkFbyNFV9DCjn40+NpAp/p6qqiJyFs+y/JiLva/Tk4kz5rwGrVPU6b9PtwPuAq4P/tzV6LsMwioRRPr6f3iz22lAFzVFxkyQkVfgvi8jlwF8DfyEiPRBT8652Tgr6fEJEHgvaPo1T9DeJyAXAbyiWVsyccOArtKD8UNByefSjlBusrUR4Pv8H6g/Cxf1Yi/uWRneNn7eD9WWKU26e/Yq5atqMu1ZeVdYajxYnr2Uwtl1uAEmt+2jerDToyPTIwHuAvwI+qKqbRGQc8PlGT66qyygfw/+2Rvs3mstRS+fbzcHobjpR4QdKfgkQ1D/iBaBMnab2Y92S3tgC4lG2zZw6yDdaKCreX3lyVT1WfRxRS7+SFeZv8wdct07ct5Cr/8UJPTWdP/Tx2wSsfFHue+BP0Ir+jy63K0kt9/TDtvOVGC0JSaN0PgTMAfYDDsfFyv8XZoV3PeGNxCz9zsNPSpYkJUNXkoGFH4x9ng8cpqr/HHhcRqvqQ9WOTerSuRA4ARcjj6quEZHY2Ph2wfd/D11W2cqNTliBwRNbts2cWsxd3zuJJ66dm0pRhjAS46il80vy5DdimfnW+cjFe/IiLmInzNlfjlUz53EU8yvuY+SLOAUd3qQ7JSJr68R9C5W+ookJo9F0qVr56U28qpUv41LRvBX4Z+BlXKTjG6sdmFThb1fVV8P4eBEZStt5rwzDMFImG4U/VVWPF5FHAVT1D0H8flWSKvwfi8ingeEiMh34CPD/6pO1O8ii5FoziU66MjqPGZOvLHlyHR6zTxgNFFrSXe/qycbs3RFESioUao7vSnJgUoV/GXAB8ATwd8CdwFdrlzM/RN0aPnEunJCoKycuXLIZpJlFE0rfc6HY+tqgbWZqpzFySLGQfbEtzkVY6XeQJ5YtubRwI4rLFBotfZoq2Sj8L+GCZg4UkQXAu4FEyiFplM4uEVkM/ERVs6k0YMQSPklEE2GlRVLL3kI0u5tbnzmusHz24Y9mKEkLUTJx6ajqDSLyMMWgmXep6qokxyaaJiYi7wQeA+4K1t8gIrfXI2yriN7pN809sfACZ+WMn7eDsbP6S0ITe3onMWLNlsSDolsn7luw7tdPH1lTCoVaGehfXXjVS2jdbZ79SuxjeNc+mncRd628KtZAaOfsmMuWXDrouxu+n9ELH2D0wgeaUifAlTms/moCewA9OB0e532LJem84Hm4KJ0tUMj7MKE2+Yy8YX55oxm0Sw6eVNglyV4pIiL/hEsbvx8wClc4JT2XDi6XzotJs1hmgQzfDf5UXI9awYWQSYBI/vAxfZsHJSKDwX7NaB3R6PqINW55/fSRpedLCT+BV62DwqG7JToOYKkUuhffIp4x+cqqaRf88apobeXrOL95gtaAnzRwxJot7rfI4HoBPJHO+SQbH/57geNU9U8AInI18AhQ1a+bVOH3i8hfAT0iMhH4GJBNejrDMIw8kCDXfZNYC+xO0cTdDXgmyYFJFf5HgSuA7cCNwN0kuJsY+cB33YyvcX/DqIeoW6czx4Ykqzj87cCTQTlYBaYDy0TkSwCq+rFyB1ZV+EG853xV/QRO6bcF1fJrhNt6eieVzNSLI4w5HkHyQZ9oycG0qTcqp57yiT6++8e/MZhbqL0Jv0+V/O/l3D3RHFNZD/76dQLG9G0uBFI0w82akYW/lNJcZvcnPbCqwlfVARH5szqEMtqM6UPOYd2S3rqObXZ4qNF8umqwNS2yUfibgTtVNdFkK5+kLp1HgzDMm/GSq3sVqjJn4tFj6FsRb01XG+CsNDgbxX8yKBe6uc/aAbbNnJqbx9hVM+eVWOMLJv+AeV/+m+IOdU60KgwEd9isYqN2srbq46hkfLj5Sg2ipB6Bk5DzgC8GGYy/kTQGH5Ir/P1wd5W3em0KNE3hi8gM4Iu4WNOvqurVzTqXkZyS8QAvmZvR3ph1Xx9pRemIyFpcErQBXFRk2XKHqjpbRPbGRet8Q0QU+AbwXVV9udJ5ks60bUpB3XIE4wb/iRuMWAf8XERuV9WGKrn7/uvx88r746PpEnp6Y8K6Ygit/+H9+ayBGnLFynex4CPf9lqDy0BWAAAgAElEQVRcrv5qE1PGz9vR8BiA0T5UmnwYrZ2bRwu/JaTr0nmLqr6Q6LSqLwUW/nDgYtxz+idE5Euq+u/ljkuaD/9LMc0vAitUtRn1Zk8AnlbVXwXn/x5wFtCQwjcMw2h3gswHH8DVJvkOcIKqPi8iewCrgMYUPi7m83U4Hz7ALOBJ4AIReYuqXlyv8GUYA/zWW18HDAqjEZE5uMIsjBs3LmURuoew0EUaU89vfeY4rlj5LmDw2IFF8nQv0aLpeRnfapQaXDqjRGSFt75IVRd56wrcE7hnvhLZFuXdwEJV/YnfqKp/FJEPVhIiqcI/Anirqu4EEJHrgXtwLpeU5qyVEDcSMujSBhdlEcCUKVOqXvrQ5zzQv7ow2zYs8uyHcPqKb/PsVxgR6KmtE/ct+5jb7DDMRgmVbaiAQ6W8YPIPCvs8O39YiV9+oH91SdTOs/OHlZQ4DN074TFzlv6wie/AyIq4733Xu3JCksfhv1DJLw+cpKobgsJSfSLyVFShe2yMbhORz6rqp1T1vkpCJFX4Y4A9cW4cguWDg5DN7Qn7qIV1wCHe+lhgQxPOYzSRabOugdnFdT/feqdYeEZ9zJh8JQP9q0tSNLTdd0JJmIU+QVeqG4L/z4vIUpxbu5zCnw58KtL29pi2QSRV+J8DHhOR+3HW95uBfxWRPYF7E/ZRCz8HJorIocB6XBjSX6XVedSij07QCi3XkYtd3vj10102zX3WDhS2DV22t8vB4x2bV+veZ9XMecyYfGXhfVyx8l0sWBNa+U/DUrh+4hEAfHjN01yxsjQu368j4D8NPDt/GItmngnAgqU/4LpPnj9oHyO/RAuaVCp0XotVX+7JoFNII0on0KNDVPXlYPk0XOnC6H4fxhWfOlxEVnqb9gJ+muRcSaN0viYid+LuOgJ8OrwjAZ9I0kctqOpOEfkHXAqHHuDrqvpk2ucxsiP06cJgBdJ2lp6RCuFYUjsYTgXSidI5EFgaJKccCtyoqnfF7Hcj8CPg33BFqUJeVtXfJzlR0igdwSXbL1RJF5ETklRJrxdVvRNXWath7lp5VeHLBLBuSS9jZyU7dvPsVxi6zBX5fnFCT2E5nKbd0zuprWeXjp+3g+smFjMdLltyKWcHj6nTZl3Dgs/dABR9/qXVwYrWmh+yed0nzy88CXx88g8LbTA4y6iRH8pVuAp9+NFQzHL7VlofsWZLbMqTaPHxtiEFhR9EIx6bYNce4CXgwugGEdkvidJP6tKpu0q6YRhGJyLa8vTID1O8xURHixU4rFoHSRV+3VXSjfZi2qxrCtb5SPassnfzsNw8RlvQwtQKqnpoo30kVfh1V0nPK2Gpw7CwcRh+OHLxnoxcXNxv5OI9eXFCcd3PuNfMcobNxFei0xZfU7Gg+3WfPN+5eQ4PGirk3Qmn5/v9XbHyXSWhn3HuHAvxywe+W2VrxL1T7bMpuGmqZZ3tnRQ736NdCqZHaaWFLyKvU9WnROT4uO2q+ki1PpIq/LBK+gG1Vkk32h8/nNJXzr5it0lVRlfSWpfOx3ETTa8tI8lbY9pLSBql41dJF2qokp4X+nbdXJh0FIZUAoXShuFgLF6xwxcn9LBz2kuRjty/drXuoyxbcmnhuriB11cKVv7wpctrmn0bFoMfumxvnrh2bqH91md+wMeDwd9FM88say2OWLOl8PRhrpz8EA7Yhjf4Ed69vfBZxlj21Z7c3ES+dGVtKS324avqnOD/W+rto2IRcxHZL3wBzwPfxYUGPRe0GQYAx1yyMJV+LNWy0VZowleKiMgwEfmYiNwSvP5BRBJlNaxm4YejwgKMA/4QLO8L/AZoeBChlRQnAQ2ufBP68tvVl9gIBXdM6J9PkB9/2qxrSvz9/jhHlLMPf7SwPG+6uzE8ce1cU+45om/XzSWhy8OXLi/8FqIWejQFR5Ry1a+iT4tjZ/W3/+8tmwIo1wPDcNGTAH8dtP1ttQMrKvxwVFhE/gu4PYiNR0TeDpzagMCGYRhtT4vDMkPeqKp+3P7/iMjjSQ5MOmj7RlX9+3BFVX8kIv9Si4Tdhj/QCd03ezSagTPqk58+5JxUsnMaRqZko/AHRORwVX0GQEQOwx98rEBShf+CiFwJLMa9xdnE+UU6AP9Rdp+1A2yeVrq9WAAkMpjbZURvYNOHnDMo1LUc4faB/tXtNY2+w4lz64REax0XM6XWNnO6NANrb2HQNnQDtdVgfesnXoVcCvyviPwqWJ+Ay49flaQK/73APFxopuKyuL23NhmzpZX+4m6xXn3lYBhdSTYKfyTQi1P0ZwEnUsxkXJGkYZm/By6qU7i2I7RsXM6d/kJ71MqpRqdlBqxGaLlvmntiiRXvM6Zvc7JnTyNXJP0dDF+6vDBpa/jS5SXhmnE5dEK2Tty3Pd2e2Sj8f1TVm4O6ttNxcfnXE1MkKkpFhS8in1HVzzS6TzfQLREn9Vr14SzcEZH2GZOvTP0xPq4gd1sqEyPXCJm5dEKb6Uzgv1T1NhH5TJIDq1n4fysilZzVgstVn+hkWRNaGL4FWg0/d/7QZXsPnohlAO46+b7cSv78bnB3dSLh51YaTrmFzbOLk/WqhVlumzm1ZLJVWz8FK0g2CWbWi8hXcJGSnxWR3agypyqkmsL/b1xy/Wr7GIZhdB/ZWPjnAjOAa1R1i4gcRMK6JNXi8OdX2t4IIvJ54C+BV4FngA+o6pZg2+XABbhHl4+p6t2Nni+sXQvJ/chjZ/WX7Odb96FFE+cq6EYLdqB/NcMDN+9wipbetplTy+Y6b1V6CnPnJCcuaiouCCH8TAeAsbOSTVjc5vv26ZBJjhkofFX9I3Crt74R2Jjk2ESPAU2iD+hV1cnAL4HLAUTkaJyb6PW4u9iXg0ydhmEYuSLMiV/tlReShmWmjqre460+iMvACS7M6Huquh34tYg8jSut+LMWi5gq0UHddok3jubIGd2kc/jJ1gyjbciRMk9CZgo/wgeB7wfLY3A3gJB1QdsgRGQOLl0o48aNq3qScPA1qctl09wTS/Lfj53VzzNfeBMAG06GX33skqp9hIq9W6J4QsLH9k1zT4SZU1ua895cOOnSt+vmst9fP6ghDj/javjZ+66ctq6BkPKgbeDJWAGsV9V3pNdzkaQ1bY/ExXkeqKq9IjIZeKeqVjRTReRe4o3CK1T1tmCfK4CdwA3hYTH7x95HVXURsAhgypQpbXavzT9pZcCs9VxJrP0w9NKUu5Ep6Wqdi4BVwN7VdqyXpBb+f+NGgb8CoKorReRGoKLCV9WKCdZE5H3AO4C3qWp46dYBh3i7jQU2JJSzLP60cX8At5J14lv3AM984U0cfrH38PGx4qLfTycO2kbrByRl9MIH2DZzallLbp+1A4U8+ka+Kfdb6emdVMiVP37ejsJnHYbqhk93ftoSn3a+aaflnxeRsbi4+gW4QidNIemg7R6q+lCkbWcjJxaRGcCncE8Kf/Q23Q6cJyK7icihwEQgem7DMFqEpdCoQPJ8+KNEZIX3mhPp6QvAJ2ly6dhakqcdTrGm7btJGAZUgf8AdgP6RATgQVX9e1V9UkRuAn6Bu6lcqKqZzMb3J2o56v8s2mWQthzrp48c9MSTlHKW3T5rB5yvP9hWi6XfzlZhu+I/uUZLXvr44zVhtSygYr3btqS24iYvqOqUuA0i8g7geVV9WEROSUe4eJIq/AtxvvLXich64Ne4jJl1o6pHVNi2APdoYxiGkUuE+AHHOjgJeKeInAHsDuwtIotVtSEdG0fS5Gm/Ak4VkT2BIar6ctqCGEZIWGPXCqMbeSeNKB1VvZziPKRTgEuboeyhevK02MGDwAWDql7XBJmaRjiL0PdJVgsr810Rh1/8YCEssxJhf52U6z105yQZ7K6G77p5seAuc1iuovwR/kZ8l07otgm/ByMXF3Pr+AP0cYP1HeWOa7PYwGoWfphHZxLwRtyAKriUCD9pllCGYRhtQcoKX1XvB+5Pt9ciiXLpiMg9wPGhKydIxdnW5muSSVjhtoLVudBZ+XGWeydZ842yKWK1j174QOL8KUOXeSHICYqpG80n/K2E/7d5E+nC9ZA4i97ft92DF0rIWdqEJCQdtB2HS3IW8iqu2ophGEb30qEK/zvAQyISljicCXy7aVI1mWjtznJW/rPzhzF+3hY24+p3rlvS21UDiX4u+zAutmC9Ty9a8dXCNf36A9HQy33WFiNuX5zQU3fop9E8ohlPhy9dXvg+xD2xRdNotHX6hCpklA+/bpJG6SwQkR8BfxE0fUBVH22eWIZhGPmnI106IjIOeAFXxLzQpqq/aZZghmEYuaa2iVe5IKlL54cU39pw4FBgNS5nfccyft6OYKCqtuLlnUhcWchK7he/oPmYvs2sjwzkhvguHnPntD/Dly6HwEXasYO1Pp2o8FX1GH9dRI4H/q4pEhmGYbQBGRYxr5u68uGr6iMi8sa0hWkl/iSsaqGZ3TRQ69O36+aStMW+lR+9Zr4F5w+Ij174wCDrPhyojQ7ghiUPrRhK/unxrHh/ULbTsmFWpRMVfmTG7RDgeOB3TZHIMOogzI8f0tFKpkUUbtydUHu2GSjIrvbS+Ekt/L285Z04n/6S9MVpPTZhqjJRv32S6+WHvYZhnH4/0YlZRvvhT1yMWvl+CGen/7461aXzC1Ut+eRE5BzafLatYRhGQ3Sowr+cwco9rs0wcoE/jtDpVqaRHR1l4YvI24EzgDEi8iVv0940WPHKaA/qVZbR46aVGaiNUu+A7bIllxb8+CPWbKm5FKMxmHJ5j/wMswP9qwftt23m1O4ZQ2kzhV+txOEGXBX1PwEPe6/bgdPTEEBELhURFZFRXtvlIvK0iKwWkVTOYxiGkSpB8rQkr7xQLVvm48DjInKDqqZu0YvIIcB04Dde29HAebhJXQcD94rIkVmVOTTSIbT4/DDPkLTCMOPK7RnpUS58OTpo2y3WvdB+uXQqWvhBbVmAR0VkZfSVwvkX4gr3+vfAs4Dvqep2Vf018DRwQgrnMnJK3E0giz4Mo2ZUk71yQrVB24uC/+9I+8Qi8k5gvao+HlbQChgDPOitrwva4vqYA8wBGDduXNoiGk3giWvnFkoYQiT/fYr09E7q3On8GeD77aOW/kD/apat7M6B8Ty5a5JQzaWzMVj8iKp+yt8mIp8FPjX4qJJ97gVGx2y6Avg0cFrcYXGilJFvEa64OlOmTGmzS28YRlvTwcnTpjNYub89pq0EVT01rl1EjsElYAut+7HAIyJyAs6iP8TbfSxu8NgwYjF3jpEVafjwRWR3XMnY3XA6+RZVbUo+l2phmR8GPgIcFvHZ7wX8tN6TquoTwAHeedYCU1T1BRG5HbhRRK7DDdpOBB6q91xGfhk7qz+1Wbdj+jYXcvGA5eNJm2gRc3OZOVIatN0OvFVVt4rIMGCZiPxIVR+sdmCtVLPwbwR+BPwbcJnX/rKq/j5tYQBU9clgsPgXuFj/Cy1Cx4hyzCULi+mUPUVvpMOtzxwHnJ+1GPlGSWVAVlUV2BqsDgteTXEWVfPhvwi8CLwXQEQOAHYHRojIiLQKoKjqhMj6AmBBGn0b+SPMPjqdcxi98IGGZ8KGin9M32azOpvA8KXLK2aU7WbSGrQVkR7cHKcjgP9U1eVVDqmLahOvQmH+UkTWAL8GfgysxVn+hmEY3YsmfMEoEVnhveaUdKM6oKpvwI1ZniAiTam6lHTQ9irgTcC9qnqciLyFwOo3jHppxLKPq45l1n26+FkvjcHUWADlBVWdUm0nVd0iIvcDM4D+uoUrQyILH9ihqpuBISIyRFX/F3hD2sIYhmG0DUknXVXx84vI/iKyb7A8HDgVeKoZIie18LeIyAhc6NANIvI8ljytaVimRyMrdm06MljaM1M52oWUonQOAr4V+PGHADep6h2p9BwhqcI/C5dAbS5u6H4f4J+bIZBhJMHcN83l+olHcPeGxwE447RzM5Ymv6QxaKuqK4HjGu+pOkmLmL/irX6rSbIYhtFmzJh8JQBzlv6w0Hb24Y9mJU5rUaCTShyKyMvEx4MKLny0OYlQDMPInNMPPhaAnph4EXvCCmgvfV81Dn+vStsNwzC6mY5KnmZkgw3UGlkxZPQvg6Vi4EA0rYJPpW1dQY5SHychaVimYRhGAT+SrGtRF6WT5JUXzMI3DGMQfbtuLlHqYS78uP0c3efTdxOv2svCN4VvGIZRLzmy3pNgCt8wDKNOzMI3DKMjsOCBKnRwxSvDMAyjBEXabOJVplE6IvJREVktIk+KyOe89stF5Olg2+lZymgYhlGWFJKntZLMLPwgxfJZwGRV3R4UV0FEjgbOA16PK3F4r4gcaVWvDMPIFZqvkMskZGnhfxi4WlW3A6jq80H7WcD3VHW7qv4aeBo4ISMZDcMwytNmFn6WCv9I4C9EZLmI/FhE3hi0jwF+6+23LmgbhIjMCSvI/O53v2uyuIZhGBGSV7zKBU116YjIvcDomE1XBOd+La6S1huBm0TkMNx8hiixl0xVFwGLAKZMmZKjy2oYRjdgYZkeqnpquW0i8mHg1qBi+0MisgsYhbPoD/F2HQtsaKachmEYNaPAgCn8pPwAeCtwv4gcCbwGeAG4HbhRRK7DDdpOBB7KTErDMAZRrIwVTzEJW+ciqFn4NfB14Osi0g+8CrwvsPafFJGbgF/gyiheaBE6hmHkElP4yVDVV4HZZbYtABa0VqL0qWQFpWkBxZ1nyOhfVrXCaiUqc639d4PV1w2k/b1qa1JQ+CJyCPBt3HjnLmCRqn6x4Y5jsJm2hmEY9aCklTxtJ3CJqj4iInsBD4tIn6r+IpXePUzhG4Zh1EkaPnxV3QhsDJZfFpFVuFB0U/hZEn2UbcRtkvS4et0ozXjsbrTPeo43N1A+MDdOHAq70p1qKyITgOOA5al2HGAK3zCMVOi6m7NSiw9/lIis8NYXBfOICojICGAJcLGqvpSOkKV0jsLf0d9yK6QV5+t2y8rC/1pLJ37fmvqekhv4L6jqlHIbRWQYTtnfoKq3piBZLJ2j8A3DMFpMGj58ERHga8AqVb2u4Q4rYArfaGt8662brf1arNhy16ldrftM5U4nDv8k4K+BJ0TksaDt06p6Zxqd+5jCNwzDqAdVGGh80FZVlxGfQyx1TOEbhmHUi820NQwjjry4TJolR9hvs11rebmOgCl8wzCMrkCBNqtpawrfMOogV1Zmzkh7ID2/11pB26vGoSl8o+3o5mgcI2eYS8cwmkt+LT4jSkd/VkoqUTqtxBS+YRhGvbSZhZ9ZEXMReYOIPCgijwWFyE/wtl0uIk+LyGoROT0rGQ3DMMqjTuEneeWELC38zwHzVfVHInJGsH6KiBwNnAe8Hlfi8F4ROdKqXhmGkSuU1LNlNpvMLHzc5do7WN6HYqHys4Dvqep2Vf018DRwQszxhmEY2WIWfmIuBu4WkWtwN54Tg/YxwIPefuuCtkGIyBxgDsC4MTYcYRhGi8mRMk9CU7WkiNyLq9MY5QrgbcBcVV0iIufissWdSnxOidirGuSTXgQw5djd2+vKG4bR3qiiA+3laW6qwlfVU8ttE5FvAxcFqzcDXw2W1wGHeLuOpejuMQzDyA9tNtM2Sx/+BuDkYPmtwJpg+XbgPBHZTUQOBSYCD2Ugn2EYRmXMh5+YDwFfFJGhwJ8IfPGq+qSI3IQr4LsTuNAidAzDyB2afk3bZpOZwg9yQP9ZmW0LgAWtlcgwDKNGcmS9J8FCWwzDMOrCBm0NwzC6A0uPbBiG0UW0WXrkLKN0DMMw2hYFdJcmelVDRL4uIs+LSH8zZTaFbxiGUQ8aFEBJ8qrON4EZzRXYXDqGYRh1k8R6T9SP6k9EZEIqnVVAtM3CisohIr8DXgFeyFqWCKMwmZKSR7lMpmS0m0zjVXX/RjoXkbuCcyRhd9x8o5BFQWoYv78JwB2q2tuIXJXoGAtfVfcXkRWqOiVrWXxMpuTkUS6TKRndKJOqNt0FkzbmwzcMw+gSTOEbhmF0CZ2m8BdV36XlmEzJyaNcJlMyTKYGEJHvAj8DJonIOhG5oCnn6ZRBW8MwDKMynWbhG4ZhGGUwhW8YhtEldITCF5EZIrJaRJ4WkcsylmWtiDwhIo+JyIqgbT8R6RORNcH/1zZZhkHTtCvJICKXB9dutYic3kKZPiMi64Nr9ZiInNFimQ4Rkf8VkVUi8qSIXBS0Z3atKsiU2bUSkd1F5CEReTyQaX7QnvV3qpxcmX6vco2qtvUL6AGeAQ4DXgM8DhydoTxrgVGRts8BlwXLlwGfbbIMbwaOB/qryQAcHVyz3YBDg2vZ0yKZPgNcGrNvq2Q6CDg+WN4L+GVw7syuVQWZMrtWuDrTI4LlYcBy4E05+E6VkyvT71WeX51g4Z8APK2qv1LVV4HvAWdlLFOUs4BvBcvfAt7VzJOp6k+A3yeU4Szge6q6XVV/DTyNu6atkKkcrZJpo6o+Eiy/DKwCxpDhtaogUzlaIZOq6tZgdVjwUrL/TpWTqxwtkSvPdILCHwP81ltfR+UfSLNR4B4ReVhE5gRtB6rqRnA/aOCADOQqJ0PW1+8fRGRl4PIJXQItlymY1n4czkrMxbWKyAQZXisR6RGRx4DngT5VzcV1KiMX5OR7lTc6QeFLTFuWsaYnqerxwNuBC0XkzRnKkoQsr9/1wOHAG4CNwLVZyCQiI4AlwMWq+lKlXWPamiJXjEyZXitVHVDVNwBjgRNEpFK+l5ZdpzJy5eJ7lUc6QeGvAw7x1scCGzKSBVXdEPx/HliKe2R8TkQOAgj+P5+BaOVkyOz6qepzwQ92F/DfFB+vWyaTiAzDKdYbVPXWoDnTaxUnUx6uVSDHFuB+XCrf3HynfLnycq3ySCco/J8DE0XkUBF5DXAecHsWgojIniKyV7gMnAb0B/K8L9jtfcBtGYhXTobbgfNEZDcRORSYCDzUCoFCZREwE3etWiaTiAjwNWCVql7nbcrsWpWTKctrJSL7i8i+wfJw4FTgKTL+TpWTK+vvVa7JetQ4jRdwBi6a4RngigzlOAwXBfA48GQoCzASuA9YE/zfr8lyfBf3KLsDZ9VcUEkG4Irg2q0G3t5Cmb4DPAGsxP0YD2qxTNNwj/QrgceC1xlZXqsKMmV2rYDJwKPBufuBf6r2vW7R51dOrky/V3l+WWoFwzCMLqETXDqGYRhGAkzhG4ZhdAmm8A3DMLoEU/iGYRhdgil8wzCMLsEUvmEYRpdgCr+LEJGt1fequc93SpCSWkTeJSJH19HH/SIypcb9V4vIO2O2TRAvBXOnIyLvF5GDvfUbROT3IvLuLOUy8okpfKMhVPV2Vb06WH0XLgVtKzhfVZs6o1pEeprZf0q8HygofFU9n4xmmhv5xxR+FyKOz4tIv7hiLe8J2k8JrOdbROSpwFqUYNsZQdsyEfmSiNwRtL9fRP5DRE4E3gl8Pig6cbhvuYvIKBFZGywPF5HvBdkMvw8M92Q7TUR+JiKPiMjNQRKxau/nz8QVwfgZcKHX3hO8z58H5/q7oH2IiHxZXNGMO0TkztAiFlfA5p9EZBlwTvA+7hKX/fT/ROR1wX77i8iSoO+fi8hJQfvJUiy88WiYaqOM3J/wZJvvtf8gON+TEmRcDd7LN73PbG4g8xTghuB8w8udyzCAzkitYK9kL2Br8H8W0IcrHnMg8Btc4Y1TgBdxSaWGAD/DTfXfHZdW9tDg+O8CdwTL7wf+I1j+JvBu73z3A1OC5VHA2mD548DXg+XJwE6c4hoF/ATYM9j2KYLp8pH3Ueg3WF8JnBwsf56gyAowB7gyWN4NWIErfPFu4M7gPY4G/hDKjStg80mv7/uAicHyVOB/guUbgWnB8jhc7huA/4fLmAowAhha5rM4DViEy+A4BLgDeHOwbb/g/3BcyoCRwJ/h0v+Gx+8bdy3iPgd72St8DcXoRqYB31XVAVzGwx8DbwReAh5S1XUA4vKMTwC2Ar9SVzQCnMKfM6jX5LwZ+BKAqq4UkZVB+5twLqGfBg8Wr8HddMoiIvvglN+Pg6bv4FJTg1Oqkz1/9j64hFnTgJvVZVPcJCL/G+n2+0HfI4ATgZsDecDdOMAl6jraa987sOZ/ClwnIjcAt4bXMobTgtejwfqIQLafAB8TkZlB+yFB+2rgMBH5d+CHwD2VrothxGEKvzuJywsest1bHsB9RyrtX4mdFN2Gu0e2xSVxEpwV+94aziFl+gq3fVRV7y5pFDmzSp+vBP+HAFvU5VuPMgT4c1XdFmm/WkR+iEt49qCInKqqT5WR7d9U9SsR2U7B3Uz+XFX/KCL3A7ur6h9E5FjgdJzb6lzgg1Xeh2GUYD787uQnwHsCv/D+OIu7UprYp3DW5YRg/T1l9nsZV4c1ZC3OFQHOjeKf/3wAcQUrJgftDwInicgRwbY9ROTISm9EXR70F0VkWtB0vrf5buDD4vLLIyJHiktbvQyYFfjyD8S5suL6fgn4tYicExwvgdIFZ2H/Q7iviLwh+H+4qj6hqp/FuZBeV0b0u4EPhmMUIjJGRA7APYX8IVD2r8M99SAio4AhqroE+EdcfWAYfM0Noyym8LuTpTi/9+PA/+B81pvK7RxYsR8B7goGM5/D+fqjfA/4RDBYeThwDU7hPoDzz4dcD4wIXDmfJLjZqOrvcGMC3w22PUh5henzAeA/g0Fb3+L+KvAL4BFxoZpfwT2xLMGlaA7blpd5P+BuIBeISJjyOqyX/DFgSjDg+gvg74P2i4OB1ccDWX4U16mq3oMbB/iZiDwB3IJT3HcBQ4P3/y/BNQBXiu/+wM32TeDyoP2bwH/ZoK2RBEuPbCRCREao6lZxTuv/BNao6sKMZLkfuFRVVzTQR/h+RuJuOCdVuum1EyLyTdyg+i1Zy2LkC7PwjaR8KLAun8S5Hb5SZf9m8nvgmxIz8aoG7gjez/8B/9JByv4G4GTgT1nLYuQPs/ANo4mIyDG4yCGf7ao6NQt5jO7GFL5hGEaXYC4dw7fMz94AAAAeSURBVDCMLsEUvmEYRpdgCt8wDKNLMIVvGIbRJfx/VPqez1Ey2nUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ds_soil.soiltype.where(ds_soil.soiltype > 0.0).plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "def fill_miss2d(ds, lat_i, lon_j):\n",
    "# fist cycle\n",
    "    cycle = 1\n",
    "\n",
    "    while True:\n",
    "        bl_i, bl_j = lat_i - cycle, lon_j - cycle\n",
    "        tr_i, tr_j = lat_i + cycle, lon_j + cycle\n",
    "        # Bottom left i\n",
    "        if bl_i < 0:\n",
    "            bl_i = 0\n",
    "        # Bottom left j\n",
    "        if bl_j < 0:\n",
    "            bl_j = 0\n",
    "        # Top right i\n",
    "        if tr_i > ds.lat.shape[0]-1:\n",
    "            tr_i = ds.lat.shape[0]-1\n",
    "        # Top right j\n",
    "        if tr_j > ds.lon.shape[0]-1:\n",
    "            tr_j = ds.lon.shape[0]-1\n",
    "        \n",
    "        lats_i = np.arange(bl_i,tr_i)\n",
    "        lons_j = np.arange(bl_j,tr_j)\n",
    "        for lats_i_idx in lats_i:\n",
    "            for lons_j_idx in lons_j:\n",
    "                if ds[lats_i_idx, lons_j_idx]>0.0:\n",
    "                    ds[lat_i, lon_j] = ds[lats_i_idx, lons_j_idx]\n",
    "                    return\n",
    "        cycle = cycle + 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [],
   "source": [
    "def count_mismatch(ds_grid_spec, ds_data):\n",
    "    count_lgs_odata = 0.0\n",
    "    count_ogs_ldata = 0.0\n",
    "    for i in range(ds_grid_spec.shape[0]):\n",
    "        for j in range(ds_grid_spec.shape[1]):\n",
    "            if ds_grid_spec[i,j].values > 0.0 and (not ds_data[i,j].values > 0.0):\n",
    "                count_lgs_odata = count_lgs_odata + 1.0\n",
    "            if ds_grid_spec[i,j].values == 0.0 and (not ds_data[i,j].values == 0.0):\n",
    "                count_ogs_ldata = count_ogs_ldata + 1 \n",
    "    return count_lgs_odata, count_ogs_ldata"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [],
   "source": [
    "def fill_mismatch(ds_grid_spec, ds_data):\n",
    "    ds_data_o = ds_data.copy()\n",
    "    count = 0.0\n",
    "    for i in range(ds_grid_spec.shape[0]):\n",
    "        for j in range(ds_grid_spec.shape[1]):\n",
    "            if ds_grid_spec[i,j].values > 0.0 and (not ds_data[i,j] > 0.0):\n",
    "                count = count + 1.0\n",
    "                fill_miss2d(ds_data_o, i, j)\n",
    "            if ds_grid_spec[i,j].values == 0.0 and (not ds_data[i,j].values == 0.0):\n",
    "                count = count + 1.0\n",
    "                ds_data_o[i,j] = 0.0\n",
    "    return ds_data_o"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [],
   "source": [
    "ds_soiltype_o = fill_mismatch(ds_grid_spec.AREA_LND, ds_soil.soiltype)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x7f1c4a6b4ad0>"
      ]
     },
     "execution_count": 55,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAELCAYAAADawD2zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnX28VWWV+L+LCyqKSOILCgKKV9Ku+DIkjTJpJUbaaEiajcxYOTGjTippviQzxIzMWKlkM5PTnV6s0ExD0p8Vis1YQyaGb3gVCV/QAPEFQ8UI4bJ+fzx7n/ucffc5Z59z9jn7vKzv/ezP3a/PXmeffdZeez3rWUtUFcMwDKP1GZC1AIZhGEZ9MIVvGIbRJpjCNwzDaBNM4RuGYbQJpvANwzDaBFP4hmEYbYIpfMMwjIwRkYtEpEdEnhSRi2t1HlP4hmEYGSIiXcBngWOAI4CPikhnLc5lCt8wDCNbDgUeVNU/qup24JfAtFqcaGAtGs2CvfbaS8eOHVt8p209ffODumoqj9E8rH5qXW6+87CRGUpi1IuHH374NVXdu5o2PvyB3XTj673Jzrdi65PAn7xV3araHcz3APNEZDiwBTgZWF6NbIVoGYU/duxYli8vfo12bDgkNz9gRE2up9GETJ0wOze/ePnVGUpi1AsReaHaNja+3stD94xOtG/Hfqv/pKoT47ap6koR+TKwBNgMPA5sr1a+OMylYxiGUQEK7Ej4V7It1W+r6tGq+n7gdWB1LWRuGQs/CQNG/C5rEYwGZPEKs+qN8lGUbZrMpVMKEdlHVV8RkdHA6cCfp9JwhLZS+IZhGGmSxHpPyMLAh78NuEBV/5BWwz6m8A3DMCpAUXpTSi+vqn+RSkMlMIVvGIZRITtornoipvANo03YseEQTj7pTH52721VtZNGG2mSVd+cAr2m8A3DMNoDs/DLRERmAX+Le2A+AXwa2BX4ETAWWAOcWatODD82P45KrIfDL5kPwONfuLFfO6XOV+k5oxQ7j0UrtTb+uIJGssRrRdy9Xo97XIFtTVYiNlOFLyIjgQuBw1R1i4jcBpwFHAb8QlWvEZErgCuAyzMU1TCanqgrplzXzMknnQlk9xAJz5+lDD6KmkunAgYCg0VkG86yXw9cCZwQbP8ecD9VKvyoFdB381bTanMyZcAZBbct2XF7HSUxqsG35I0MUOhtLn2frcJX1XUici3wIi6HxL2qeq+I7KuqLwX7vCQi+8QdLyIzgZkAo0cnG+JcD564blYwN6vfNnOnGPXEt4rj1iW12huho7bU+ev923IjbZuLTFMriMi7gNOAA4H9gd1EZEbS41W1W1UnqurEvfeuKg+SYRhGmQi9CadGIWuXzonA86r6KoCI3AEcC7wsIvsF1v1+wCtZCllrfHdT1JKqt9USdROkkXYgbNNSGBQnDJtMgl3L7FFgh7l0yuJF4H0isivOpfMhXFrQt4FzgGuC/3dWeyJfcU6dMLvhXg/TppD8UyfMpqMLentW9dvW0TW+1mIZNWTxiqtj/fqbO4fl5qMdn0kfMFm7cwoRyp/FA1CBd5os/2TWPvxlIvJj4BFcOtBHgW5gCHCbiJyLeygU7mVsEsJQzZFLNgLp/4CKdeD5yr1cpT51wuyqfky+XEk6Gc1yTUbctbRO3PqzQxvHXZOErC18VHUOMCeyeivO2q8JTqk0pmIJHwS1fsMwa77x+PD+RwDu4Zzm9zN40TIAtkyblGftQ7zhUeje86O77ln/eGryNStupK0pfCNCaN2XQyNZa+XIYhZ6eRy6aC4Ao6b3lNizfuQXCmpu12YtUYRec+kYSTj5pDMLKse0lX2W1nwhv3LcfkY+vT2rKv7u4vpoQks/5MOLjshbLmW1JxmjUcu3APcGFGVVpmNHzKVjGIbRBijCO9qRtRhlYQo/IHS7jJj/ABtmHVtwv75BVcnbLIdGcuVUQqXy1+otYOqE2Tlrt1xL0FxZ/Udl20jsPtzAK3PpGDGE0TnrpgwHyntwlGLmop/SPe0UZi76KQDd005JrW2jPnR0je/nhomuK8e94+8b596JI95lUjlhe2m7dsLP1ggPWeu0bQJic8kUsep9ola7r7j9bXPP/z7wNwXbuePZowC4sfPg3BtF+FAw4qnV208a7RbLT7Rl2iQANs54G4Axc7b1bZw7qOpz15Pwc5qlD6pCr5qF39QUc+eUy9zzv8/p4x51yv18t27ON/4meBikR2jR19qy39w5jCGrN9X0HLWg1p3WSS3oOHzlHydntbJXI1slLNlxe9GHXzWEbwonn9Q4IcU7UrLw49LEq+qfUmnco60Ufs5Pn2KboaXu6G/RH37JfOaen+IJjYqZOmF2zv1lVIf59sM4/Oot/CJp4m+quvEIbaXwC7Fh1rGMmP8AAOetfiZv25xvxLtl5p7//YLbIN9Vc9WKjzFvwk/ij5tF7txErLlqwvLKIRyMU8p6b0brvhr861/Mlx73Hfn7h+GQY1aX/11WYqHH9QekyT3rH0/d35+EcEzA4hV1P3UsirBNU1OhcWniU8cUvmEYRoX0phCHXyhNfNUNx9A2Cn/KgDMSd8wahtF4hCOATz7pzIJvXJWG4FZCmSNt9xKR5d5yt6p2Q7808ZuA20VkhqouSFVgQLTJajIWYuLEibp8+fKC26M+xySds2HnajHXTTGSRt2EN2m1r+K1cP+8UCSKJC/apIEJw1WhfiGrSb5H//teu7Cr3/Zap1tI0/XjK9job63SsMxoWoc0Fb6IPKyqEysSLGDc4bvpvy46LNG+Z3UuL3g+ETkDmKqq5wbLfwO8T1VT7/1raQu/VtEChmEYaXXaUjhNfOq0tMKvliSWfXRAlb8uSc75UhZWtBP5xs6DS8pULXFWfRJrPjyuXpZ/0s7m7mmn5AbpdNO/wEsacfj+W5r/v9j3628bNb2nn5W/dmFXTax8/x4s936MI2pNh8tTBpyRyLpPmqCtb7/ZeXL7ef3rmTdKkbR8+IXSxKdOyyp8s+4Nw6glqqQWpVMgTXzqZO7DF5FhwLeALtxb0mdwScF/BIwF1gBnquofirUT9eH7Cj/OXx+GQlY60Crqn680hDJqUfkW/Y2dB+eWo77nevili/nvfZrBl59kGH5SSz/8rsPv7p71j+cszN6eVf2Wk7J2YRcDlw4FYI81vbwxti8x1x5revtlu6yEqDUe/czlyJukY3THhkNKpljIIgVzGj78sV276+w7jky072fHL636fGnQCOOCbwAWq+q7gSOAlcAVwC9UtRP4RbBsGIbRMCjQqwMSTY1Cpi4dERkKvB/4FICqvgO8IyKnAScEu30PuB+4vP4SGu1As2cozYpw5HqaiQCbDSuAUh4HAa8C3xWRI4CHgYuAfVX1JQBVfUlE9ok7WERmAjMBRo8eXfAkUffNiPkP5LlfRi7ZmFf6LcmI0hfmDiqZByUJ/Y97Juei6ehyy3F0Tzsl9/p93urapAvI0lUTuqxOH/cokE093GJuutB1A/lhguE6v3RgoVGpvlsoPM/2yW8CsHEyDF+wW27fSt05pTJLRteX6vtasuP2stJ+DxjxO5bscPOhe8d37TRzRS1FrABKBec/Gvhc0FN9A2W4b4KBC93gfPi1EdHIksnTrwVcVXujcsIHZrGH4uGXzC+ZZypO2R9+yfyGtvLjMtOmgZJep229yFratcBaVQ3Nlx/jFP7LIrJfYN3vB7xSbsN+fhxwlo6/HFpWI3EWULl5YuKs32KDQZK8AWzuHNbPUi/UWds97ZRcm93TGid7YFLiFI9vxXdPO6Vfwe1weenCS+siS6HvLFw/YMTtwfH99/HrwhbKTd/bsyrvrXPMnPGR0Fbvniwjv/2WaZNqkveoWmXpW/utgVg+/HJQ1Q0i8nsRGa+qq3ADDp4KpnOAa4L/d2YoptFGtLo/f+qE2XkP0rQeno1q5R+6aC7zJtSmbQV2NFCHbBKytvABPgfcLCI7Ac8Bn8ZFD90mIufiRqGVHVT/xHWzILgBQ79k3GCYqPUd9c0XI65CUbH1ccRZsaF1Ft0GcP1lZ7uZzr51pay5RsljX8rH7m+fOmF2nsy9Paugc1LNZIuj2PdY6rPk+6b7HiJRv30uvDdYP2aOO0ex76xYSoSw2EolPHHdLKZ4b8FRi97PKFtpupEs2TJtUu56pzWczSz8MlHVx4C4+NQP1VsWwzCMpKiKWfhZ8fSzLzN5+rW5V9TDL5mfi3gYlaVgRg7fXdII9UhDWt2NUwy/gM+dG4+uqi0/wqcRCqJMnTC75iUkGynGPgkto/BLUSrELiR8pU5yXLnrfeLcNaW2DV60LPfKPmT1ptx+cTllfJdAI7hzyiX8TH444uBFy+qqSDq6xqfyYArbCB8svusvmtcmvBcH99CvIE5UtjB/0/bJb+bckOE9EncPhdfy84viQ30hP3S50LY534h39TQiY+ZsoxsX+DCEdH8HrgBKR+kdG4iWU/jlxAgb2ZEkTNAwKiHaZ1crXKet+fAzoWPr9rzcIyPmPwCB7u/oGp/XGRvtOIt24tabQhZ4v467rvF9Fq93M8cdn5ZVnzQjZaFasdXUkI0ONqqVdV/vh05omY+kz0oPOxN9C7tUOG+uw3dJ37qOrvF5b39Q3qCtoz/hujPXz48fsBguN7JVH1JNQEbic9hIW8MwjNanGUfaZp4tMy0GH7y/Hr/rZwpuL/Z0j8tlnwb+IJroEPqQqBWybsrwojLE+Wbr5aOvRXhnO7t0/M7iqKUfrqv0fvzd374LgHEXP1h0Pz+Mc/3xfcqr1HFRovd1I32v4XVeN2V4LpDj6dO/VHX2yhHv2VNn3DIl0b7XHXlbQ2TLNAvfMAyjAlTTKWJeT0zhG4aROu3QKa8I23dYlE4m7LReeeHmvpjb4Qt2y3M/FOus2T75zVxmwmpeo/2CIWPmbMuVp3thYRejprv1caF4PtHMnVGyDLGMnjv6eY3yCV19I3ri74VKOeRbQb2gEjl4/A7dcYuKtxntAA0ZuHRoqq7QtMkLiw06uJ9Oqe1mG2lbtItZRN4sMb0lIs2b39QwjJrSyoPawrDMJFOjUMrCf1ZVjyq2g4g8mqI8FfPO/pKzMtdNGc4bY+GNsX2dsf7gFN+KDy2V0HrdOGMQ6+jfgZYE38p17bqi1AOXDs3rHIsOkPIZsnpTQSu+t2dV0cE49aRY+cNCHeR+wfBWftVPSjR3UKMSl7vH/y35/yHdwIfGJp3UCiIyHlfSNeQg4J9U9WtVNx6hlMKfnqCNJPsYBocumsuYYN4Uv9EK7EjBpRNkCj4SQEQ6gHVACQdbZRRV+Kr6XKkGkuxTb0Yu2djPynxh8iAvTLLPAon6H90x/S0U3zpNQr6Fm++Xj1r1oR81HBLv+1X9VApxFZKyYviC3dg44+3csn+946z7pAXRjcpIUrGtUDbXcJtvpUdrR/j3Xlxe/EZMjRxl8YqrOXTRXLdwevXtqcK29DttP4TzrLyQdsOQsNNWRE4HvgzsA0gwqaoOrYVQhmEY1RJroKVoeJQ58GovEVnuLXcHFfuinAX8sGrhCpA0SucrwF+q6spaCVItO61XODh/nW+xwJu5Od8aHblkYz//Y9QaXTltTlWyFYu6iSYHg3zrPWqtJbHse3tW5fUZDF60LPU3go0z3u5nyedXa+oj+nZk1n765KVPKPFdx90L/aLTItE9vT2rGiIDZrWE92ZaUTpluHReKzXwKqgJcipwZbVyFSJpj8PLtVT2ItIhIo+KyN3B8p4iskREVgf/31WrcxuGYVRCDaJ0PgI8oqov10rmohZ+4MoBWC4iPwJ+AmwNt6vqHSnJcRGwEghdRFcAv1DVa0TkimD58pTOVTaNHEFhGEZ2pFwA5ZPU0J0DpV06f+nN/xE4yVtWoGqFLyKjgFOAecDng9WnAScE898D7qeEwn9n//yn6Jg521i8wnUkhQMufFdCvlvBc/dMLjyYqJyBRn7emWi2wi3TJsVmMCzX7VKsBJ/focqMrpy8L8wdlBsQVsk5C7ltwnWLV1wN04q3kdvPSI1qXXZxoZRZZI6tNbnfcAqdtqQYYy8iuwJTgL9LpcEClIrS+XQgzHGq+mt/m4gcl5IMXwMuA3b31u2rqi8FMrwkIvukdC4jA+wNyWhFFNiekoWvqn8EhpfcsUqSdtr+OxCtfxa3rixE5KPAK6r6sIicUMHxM4GZAAP33oMXvjqoaAdruO3QRXMZuNR5j8K0CqFFnDSk8IW5gxi+YLe8jtbQIgot+M1Bx+nGuV15VnWh/OSl8p+H+5Sio2t8wQydY+ZsK9qpF82IGX2rqTaFgnXYNifN2GHb27OKtQu7apb2o+UKoIjInwPHAnuLyOe9TUOBNAJQjwNOFZGTgV2AoSKyAHhZRPYLrPv9gFfiDg7CmrrBpUeG4tZkLgbXqBv+YCvDaDVaSuEDOwFDgv18l8ubwMerPbmqXkkQghRY+Jeq6gwR+SpwDnBN8P/OUm3ttF6LPsldaGSfX3v75Df7LH76HgRxaQHCZd+63jJtknsrmNGVt0+/7bhBSn6YpE851YigPF9tuG85eewL7VfqzaeUFeXvW22YazsRlye/XjSjVe+zdmFX3nLab5fNWACllA//lyKyFDhcVetpHl8D3CYi5wIvAmfU8dyGYRiJSCO1Qj0p6cNX1V4R2bPWgqjq/bhoHFR1I26IsdFkTJ5+bX6EkGG0Ktp6Lp2QR0XkLuB2PL9IinH4qVAsY98bYzsYuHRobM4P38Vw6KK5/UIP+14Fu3J586NKzc+OOWq6c9OMCtIfFXLnADXtVArx3TTVFBUPXVrhq360TyRuRLG7Xtnl8G8W4nI15Y8Ur55alKhsZMLf9eQF1wIwfEHfbyGNkbYKbN/RmkXM98RlFPugty6VOPx6M3XC7NwPKVT+UyfMtuiRMrDOb6MZyD1Ai6Q2qYaW8+GHhPH4jYxu2Zqf4S+yPa6iUCFyFn/MAKLQWhgzZ1vBEMa1C7vyQjZ9ouvGrO5f1NwnlLlYEfZyKMe67+1ZxfAFk4Lzv034BhOS92bEXPxO8VDWzZ3Oslw57dLKhW4TovnxowXNQ6Jvr3HhvOHbVlKL3m+j1QbF1fKtRltR4QejYf8dF0apwFLgIlVdW0PZak6hEM7DL5mftxxWuh81vQeKuGcMw2gvWq7TNuC7wC30RcvMCNZNqYVQtSZqIfnWc2/Pqth83yE53/0CZ0H5YXO+z3/jjLcZjns4DF60LM+P76dW2Nw5LNemPzjLlw18i3lYrs00MmD6rqxoHWDXvxBUAqO4y2vltDl5rp5WsxLrzeIVV3P4JfPz7lV/3n8T9O+DYpZ93Dq/nVbJiFmMNNNFaBN22ibtcdhbVb+rqtuD6SZg7xrKZRiG0eAIvTsGJJoahaQW/msiMoO+TG6fJK4slGEYRsbUM3dTS/rwgc8A/wHMx/nwHwjWNRSlSrgVws87A/nl3Xz3zoZZxzJq+gN568N913n7jZr+AB1d49nc6Zaj2TGj86MKVK9cu7Arl/Nn3RT3Sp97LY8UqCj1GQsRhpmG+KGVwxeQ+wzDF8QcHOnUthG06VNopG2x73rI6k2573H98cK1p+R/eTd2HlyyjVYk7eyfLZdLJ0RVX8RVYjEMw2hY6pqZVZ0fv5lIGqWzN/BZYKx/jKo2nJUPfZZLaB2NmP9AwU6uMBwtus4ntLJ9y3/E/AfYMOvYvDeAcHu4bo81fcGh0ZKDpWT3z1uKJJZaoc5f3xo0Ggv/fitVtjB3z3bmR5Ht/0vl9AsfzVt3+o6++SkDWj9rSS3z+rdqlM6dwP8B99E/xN0wDCNTsqi5oLSuD39XVc2sxGAStu27W55FD26wVRxRa7+Q9R+2VShMM7otup9vyYeZ+0ZN7+mXasHfLzrAyvfdbu4clrdvKcu+0Ohh36I3675x8UMkkyo0/41t6cL2HOxWv0pdQu+O5lL4SeOF7g5y1huGYRgBqpJoahSSKvyLcEp/i4i8KSJvicibJY8yDMOoMVn1Q6imp/BFZJiI/FhEnhaRlUHxqdRJGqWze7HtIvIeVX0yHZEqY9DLb+d3ciUk+vpX7utg6HIpldkwrsRh3KhGf7Tq9slvwhI3H76mF+r8jY62LJQ/xJVm7FveOOPtfqGZhShnXyNdotk0fVfckNWb8ordlEOlxzUaWRVcTzEs8wZgsap+XER2AnZNq2GfpD78UvyAKuvbtiOhZdLqw9mNxmDKgDPyDIY4o2Dy9Gtz8+3aB1AOaYRlishQ4P3Ap1yb+g7wTvUt9ycthV/RY05EDgC+D4wAdgDdqnpDUHDlR7gw0DXAmar6h1LtFRt4VQ1J3hxGLtmYZ+VHj4kO1CrEymlzmDonnYiD/vn8HdFc/v4Pv5SlZ8VNssPPTzR5+rW5783PdJlWjiWjNIqwI3nahL1EZLm33B3U5AY4CHgV+K6IHAE8jEtOmfqPLS2FX+lzbjtwiao+IiK7Aw+LyBLck+4XqnqNiFwBXAE0dJSQ0do049tYaK1Xaqn71n6IWf35lKH4XlPViQW2DcR5SD6nqstE5AaczvvHqgWMOVFmqOpLwEvB/FsishIYCZwGnBDs9j1c6cOSCr/efjw/pDNuAJZPdFAWEFt9CyJ50RNa+9FqVOXEJcdZ9fbDbmyWLrw09x37Fn0xn3w0HHhz57C879lX8DYgLwGaWhz+WmCtqoadcj/GKfzUSSuNW9X+JhEZCxwFLAP2DR4G4UNhnwLHzBSR5SKyfBtbqxXBMAyjPDThVKwJ1Q3A70UkfHJ/CHiqFuImTa1wHPCYqr4dZM08GrhBVV8IBH5fNUKIyBBgIXCxqr4pkuypGfjAugGGyp5NltWi+YiWNrRkac1BrUv9tTMpxth/Drg5iNB5DqhJlcGkLp0bgSOCDoXLgG/jOluPr1YAERmEU/Y3e0XRXxaR/VT1JRHZD3il2vOkTZIRuf4+aRekHrJ6U78cF75/Oa4oto+fp2fjjL4hFRZ2GU+z+O6TljaMc9ktXXhpzq3TTO6cDbOOrSgkOw3SSp6mqo8BhXz8qZFU4W9XVRWR03CW/bdF5JxqTy7OlP82sFJVr/c23QWcA1wT/L+z2nMZhtGHH+Vj1ckqQxW0gYqbJCGpwn9LRK4E/hr4CxHpgBI175JxXNDmEyLyWLDuizhFf5uInAu8SF9pxcwJO77CH0tH13gWr7iaKQPOKJhHP0qhztpihOfzLS+/Ey7OYuvbNz+6a8ycbawrUJxy44y3zVXTJMQp6mhkTbTAeanO2GbsrE9q3UfzZqVBS6ZHBj4B/BXwGVXdICKjga9We3JVXUrhGP4PVdt+vZg8/VoGZy1EBhy6aK49HIz2phUVfqDkFwJB/SNeAwrUaWo+1i7sii0gHmXLtEn9fKO+tVRscFUlVn0cUUu/mEXmb/M7XDd3Dsvl6n9jbEdZ5w99/DYAqzGJ3g++6yb6v1VIarmnH7bdWInRkpA0SuezwExgT2AcLlb+v2giKzxr/A7UevhM+w2amVGb84QPErP0Ww8/KVk0hr8ZXT81IQMLP+j7PBs4SFX/OfC4jFDVh0odm9SlcwFwDC5GHlVdLSKxsfHNgn8DD1xa3MqN3uzQP4phy7RJfbnru8bzxHWzUinKED4cDl00Ny9PfrlWWqiQD100N886H75gN97ARexsn1w6AapZ9s1FVDFPnTA7l2qjVSKyNncOy1X6ilaTi0bTpWrlpzfwqly+gUtF80Hgn4G3cJGO7y11YFKFv1VV3wnj40VkIE3nvTIMw0iZbBT+JFU9WkQeBVDVPwTx+yVJqvB/KSJfBAaLyBTgfOD/VSZre5BFyTXDqIS4nDulAhFCd0/49tu2Lp5szN5tQaSkQq7m+I7ihziSKvwrgHOBJ4C/A34GfKt8ORuHqFvDJ86FExJ15cSFS8ZRrd8+LotmsTZL/QD9z5wrtr4mWDetsAzR0bZGc+EPyNvcWXzfYr+DRsIfMBaXKTRa+jRVslH4X8cFzewrIvOAjwOJLMykUTo7RGQB8CtVzabSgNGPqRNm26AZIzXismMm5Y5nj8rNnz7u0TTEaXyUTFw6qnqziDxMX9DMx1R1ZZJjEw0TE5FTgceAxcHykSJyVyXC1ovok37DrGNzE7jBR2PmbGPU9J680MSOrvEMWb0pcafo5s5hOet+3ZThqadQ8OntWZWbKiXs+N044+2Cw+uN1mbxiqv7GQrNHqq5dOGl/e7d8Hc5Yv4DjJj/QE3qBLgyh6WnGrAr0IHT4YmHASV16czBRencDy7vQ5Dd0mhiDl00l+EUj9QwF45RLvUOQc6UHfW38EXkn3DZBxbiBq5+V0RuV9WSF7ucXDpvJM1imQUyeGf4U99y1ArOhUwCRPKHj1yysV8iMuhv9UTriEaXh6x28+umDM+dL80b3k/glVancLmpFPzwTqO58S3iqRNml0y74PdXRWsrX8/ZffvRGG8LQ1Zvcr9F+tcL4Il0ziHZ+PA/CRylqn8CEJFrgEeA1BR+j4j8FdAhIp3AhUA26ekMwzAagQS57mvEGmAX+kzcnYFnkxyYVOF/DrgK2ArcAtxDgqeJ0Rj41viYCo8zjEqotsxiYyNZxeFvBZ4MysEqMAVYKiJfB1DVCwsdWFLhB/Gec1X1Czil37CEmSvBvaIW69wMt3V0jc8bqRdHOHpvCMk7fWpdbrFSV1G0oHmlWCqF1iK8n4pF6hRy90RzTGWZSz+vPOiE2YxcsjEXSJHn1k2LbCz8ReTnMrs/6YElFb6q9orIn1UglNGEVGLVh/0JLd9B1+JUE5bZtmSj8DcCP1PVRIOtfJK6dB4NwjBvx0uu7lWoypzOw0ayeHmfwok+6YtRrHM2iv9mUCiUbY81vWyZNqlhXmOjA6bmTfgJc77xN307FBhoZRhJacQKWcUMEDdeqUqUTKJ0gLOAG4IMxt9NGoMPyRX+nrinyge9dQrUTOGLyFTgBlys6bdU9ZpanctwTJ0wG0q4fPL6A7xkbuHxZuU3J2bdV0ZaUToisgaXBK0XFxVZsNyhqs4QkaG4aJ3viogC3wV+qKpvFTtP0pG2NSmoW4ig3+A/cZ0Ra4HfishdqlpVJXfffz1mTmF/fDRdQkdXTFhXDKH1P7in8Wqg+uHZXHUAAAAf+klEQVSUV634GPPO/763NVmu/jFztqXWB2A0PsUGY0Vr5zaihV8X0nXpfEBVX0t0WtU3Awt/MHAx7j39CyLydVX990LHJc2H//WY1W8Ay1W1FvVmjwGeUdXngvPfCpwGVKXwDcMwmp0g88GncbVJfgAco6qviMiuwEqgOoWPi/l8N86HDzAdeBI4V0Q+oKoXVyp8AUYCv/eW1wL9wmhEZCauMAujR49OWYT2Icx8mMbQ8zuePcq9QUz4CQBXrfhYbptF9mSD767Jql8prLyVpFJbM1GGS2cvEVnuLXerare3rMC9gXvmm5FtUT4OzFfVX/krVfWPIvKZYkIkVfgHAx9U1e0AInIjcC/O5ZLSmLU84npC+l3a4KJ0A0ycOLHkpQ99zr09q3KjbcMiz34Ypa/4Ns54myGBntrcOazga254fKO5cgoRKuJQMYNzefl++d6eVaxd2JW33S9xGLp3wmNmLvppzeU26k/cfd/2rpyQ5HH4rxXzywPHqer6oLDUEhF5OqrQPV6KbhORL6vq5ar6i2JCJFX4I4HdcG4cgvn9g5DNrQnbKIe1wAHe8ihgfQ3OY9SRqRNm91Me1snbnoTjZPwUDU1n9SsJs9AnaEp1ffD/FRFZhHNrF1L4U4DLI+s+ErOuH0kV/leAx0Tkfpz1/X7gX0VkN+C+hG2Uw2+BThE5EFiHC0P6q7Qaj1r00UFSfSXgXN74dVNcNs091vTmtg1cOtTl4PGObQbrfuW0OXll7q5a8THmrXZW/vWXTWLmopu5sfNgAM5b/QxXrejKOz6sIzB8wW55ivuFuYPonnYKAPMW/YTrL3O5VcY0eRbGViCJIg0HF4YK2P9uywlbjlLozaBVSCNKJ9CjA1T1rWD+JFzpwuh+5+GKT40TkRXept2BXyc5V9IonW+LyM9wTx0Bvhg+kYAvJGmjHFR1u4j8Ay6FQwfwHVV9Mu3zGNkTWv2FaDqrz6iYsC+pGQynHOlE6ewLLAqSUw4EblHVxTH73QL8HPg3XFGqkLdU9fUkJ0oapSO4ZPu5KukickySKumVoqo/w1XWqprFK67O3UwAaxd2MWp6smM3znibgUtdke83xnbk5sNh2n46h2Zh8Yqrc/H0Y+Zs4/rOvkyHp497lOunuU6+6y+bxLyv3Az0+fz9SlmbO4flK+RgANfk6cNybwKfn/DToC13Dr/jrtWsvWanUIWr0A2X5HuLPrzjKsTFpR2JFh9vGlJQ+EE04hEJdu0A3gQuiG4QkT2TKP2kLp2Kq6QbhmG0IqJ1T4/8MH2PmGhvsQIHlWogqcKvuEq60VzYiEvDKIM6plZQ1QOrbSOpwq+4SnqjEpY6DAsbh+GHwxfsxvAFffsNX7Abb4ztW/Yz7tWynGGtCWPiJy+4tl9Bd79j9vrLzmbpwks5fVywQ4K8O/l+9zl59U7j/PVxIX7m7qk/vltlc8S9U6qjNuemKZV1tmt87HiPZimYHqWeFr6IvFtVnxaRo+O2q+ojpdpIqvDDKun7lFsl3WgNQsvfV87+gwJqM7Aqmviu2fpLjBanvi6dz+MGml5XQJIPxqzPI2mUjl8lXSijSnqjsGTH7bmOyjCkEsiVNgw7Y/GKHb4xtoPtk9+MNOT+NbN177N04aV5HbgvzH071zGbZkfa54PO3+5ppxS0Foes3tRPoadVytEozpIdt+cFNviEHbbhA36I91zPfZcxln2pwVluIF8VQmdNnX34qjoz+P+BStsoqvBFZE9v8RXgh/62pKFArU6oMC11gGG0GRnkwxeRQcB5uPFQ4AqgfFNVtxU8KKCUhR/2CgswGvhDMD8MeBGouhOhnvSlDehf+Sb05ZfrS+z3BtCE5B5UoX8+xk8fWvu595/OSQxfsFuiOPnTxz2am58zZT4AT1w3y6z3BmbwomW530LUQo+m4IhSqPpV1Hc/anpP0/ruc2RTAOVGYBAuehLgr4N1f1vqwKIKP+wVFpH/Au4KYuMRkY8AJ1YhsGEYRtNT57DMkPeqqh+3/z8i8niSA5N22r5XVf8+XFDVn4vIv5QjYbvRjmX/Dr/EWe/h2xL0WXW5KI4gOmrKgDNSyc5pGJmSjcLvFZFxqvosgIgchN/5WISkCv81EZkNLMB9xBnE+UVaAP9Vdo81vWycnL+9rwBI87tyyiHN4e7hA6G3Z1XJdtvpgZk10Y5bv9Pez5oKeJlSN5WVWyc/A2tXrtM2dAM11fdd/4FXIZcC/ysizwXLY3H58UuSVOF/EpiDC81UXBa3T5YnY7bU01/sW6+tbOn3i+oIrHfDaBuyUfjDgS6coj8NOJa+TMZFSRqW+TpwUYXCNR2hZeNy7vTk1ketHCOf0HLfMOvYPCveZ+SSjcnePY2GIunvYPCiZblBW4MXLcsL14zLoRPSLy9Ts5CNwv9HVb09qGs7BReXfyMxRaKilArL/JKqfqnafdqFdog6KRSrXSm1KnweTRHRlMrEaGiEzFw6oc10CvBfqnqniHwpyYGlLPy/FZFizmrB5apPdLKsCS0M3wIthZ87f+DSof3CMP3skdC+6QA6usbn+XKjqSv8POvWWdu4hH0qcQ/28HvLD6fcxMYZfYP1SoVZbpk2KW+wVVOn0FCQbBLMrBORb+IiJb8sIjsDA5IcWErh/zcuuX6pfQzDMNqPbCz8M4GpwLWquklE9iNhXZJScfhzUxAuFhH5KvCXwDvAs8CnVXVTsO1K4Fzcq8uFqnpPtecLa9dCcj/yqOk9efv51r1v2S9deGmeO8e3YFuxszaO3p5VDA7cvIPps/S2TJtUMEVDvdJTmDunfPzoqbgQWn8g3qjpyQYsbvF9+zRvwrQ8MlD4qvpH4A5v+SXgpSTHJnoNqBFLgC5VnQD8DrgSQEQOw7mJ3oN7in0jyNRpGIbRUIQ58UtNjULSsMzUUdV7vcUHcRk4wYUZ3aqqW4HnReQZXGnF39RZRMMjN6iqRm0/cd2sRPv6nbFmuRuZ00DKPAmZKfwInwF+FMyPxD0AQtYG6/ohIjNx6UIZPXp00RP4g0qSdhpumHVsXv77UdN7ePZr7wNg/fHw3IWXJGqnHQlf2zfMOhamTSqZOTFt7GGQHkt23F4wAs0Paohj8KJlud9b+N37rpx63Q81IeVO28CTsRxYp6ofTa/lPpLWtD0EF+e5r6p2icgE4FRVLeqgFpH7iDcKr1LVO4N9rgK2AzeHh8XsH/scVdVuoBtg4sSJTfasbXxCq76e50pq6YOz9k2x1x5Lg1GEdLXORcBKYGipHSslqYX/37he4G8CqOoKEbkFKKrwVbVogjUROQf4KPAhVQ0v3VrgAG+3UcD6hHImwu/ALWad+NY9wLNfex/jLvZePi7MbzOkFX8c0foBSRkx/wG2TJtU0JLbY00vb4y1LppmoNBvpaNrfC5X/pg523LfdRiqG77d+WlLfJr5oZ2Wf15ERuHi6ufhCp3UhKSdtruq6kORddurObGITAUux70p/NHbdBdwlojsLCIHAp1A9NxGi1LPtwrDqBpNOMFeIrLcm2ZGWvoacBk1Lh1bTvK0cfTVtP04CcOAivAfwM7AEhEBeFBV/15VnxSR24CncA+VC1Q1ldH4xQaVxOEP1HKU9120UkjmuinD+73xJKWQZbfHml7n6w+2JbH0m9kabHb8N9fQio8bOOX314TVsoDYqlhN/X32KfMkvKaqE+M2iMhHgVdU9WEROSEd4eJJqvAvwPnK3y0i64DncRkzK0ZVDy6ybR7u1cYwDKMhEeI7HCvgOOBUETkZ2AUYKiILVLUqHRtH0uRpzwEnishuwABVfSttQQwjxEpGGs1CGlE6qnolfeOQTgAurYWyh9LJ02I7DwIXDKp6fQ1kqiulwsp8V0TYYVsqh3ux9pqV0J2TpLO7FL7r5o1ISuVWKBnZaoS/Ed+lE7ptwvtg+IK+3Dp+B320s76pXThxNFlsYCkLP8yjMx54L65DFVxKhF/VSijDMIymIGWFr6r344qS14REuXRE5F7g6NCVE6TiTK8EUgaElkmxEMpwW2h1rp3cVdDNkGZFqGZnQ8RqHzH/gcT5UwYu9UKQY4qpG/Un/K2E/7d4A+nC5ZC48NumzohZjAZLm5CEpJ22o3FJzkLewVVbMQzDaF9aVOH/AHhIRMISh9OA79dMqhoTrd1ZyMp/Ye4gxszZxEYGxW5vdfyaAWFcbM56n9JnxZcK1/TrD0RDL/dY0xdx+8bYjopDP43aEc14OnjRstz9EPfGFk2j0dTpE0qQUT78ikkapTNPRH4O/EWw6tOq+mjtxDIMw2h8WtKlIyKjgddwRcxz61T1xVoJZhiG0dCUN/CqIUjq0vkpfR9tMHAgsAqXs75lGTNnG709qxg13Tpl48pCFnO/+AXNRy7ZyLpIR26I7+Ixd07zM3jRMghcpH5nbSuNOs+jFRW+qh7uL4vI0cDf1UQiwzCMJiDDIuYVU1E+fFV9RETem7Yw9cTPq1MqNLNlrZME+CGWvpUfvWb+NfI7xEfMf6CfdR921EY7cMOSh+WkSDayocOz4v1O2VbLhlmSVlT4kRG3A4CjgVdrIpFhVIhVw0qX3IO7FWrP1gIF2dFcGj+phb+7N78d59NfmL442RBNE9Du/nqfqN8eSl8fP+w1fEPw24kOzDKaD3/gYtTK90M4W/231KounadUNe+bE5EzaPLRtoZhGFXRogr/Svor97h1htEQ+P0IrW5lGtnRUha+iHwEOBkYKSJf9zYNpcqKV42CKYPCVHNtosdOLtBRGyWNDtshqzeVXYrR6E+hvEd+htnenlX99tsybVL79KE0mcIvVeJwPa6K+p+Ah73pLuDDaQggIpeKiIrIXt66K0XkGRFZJSKpnMcwDCNVguRpSaZGoVS2zMeBx0XkZlVN3aIXkQOAKcCL3rrDgLNwg7r2B+4TkUPSKnNoZENo8UVr1oaZSNModtKyWRkbhELhy9FO23ax7oXmy6VT1MIPassCPCoiK6JTCuefjyvc6z8DTwNuVdWtqvo88AxwTArnMloYK35uZIJqsqlBKNVpe1Hw/6Npn1hETgXWqerjYQWtgJHAg97y2mBdXBszgZkAo0ePTltEowY8cd2sXAlD8PLfp5z7vt0HzKWN77ePWvq9PatYuqI9+8IayV2ThFIunZeC2fNV9XJ/m4h8Gbi8/1F5+9wHjIjZdBXwReCkuMPiRCkgXzeuuDoTJ05ssktvGEZT08LJ06bQX7l/JGZdHqp6Ytx6ETkcl4AttO5HAY+IyDE4i/4Ab/dRuM5jw8jD3Di1445njwLOzlqMhicNH76I7IIrGbszTif/WFWr79SKoVRY5nnA+cBBEZ/97sCvKz2pqj4B7OOdZw0wUVVfE5G7gFtE5Hpcp20n8FCl5zIal1HTe1IbdTtyycZcLh6wfDxpEy1ibi4zR0qdtluBD6rqZhEZBCwVkZ+r6oOlDiyXUhb+LcDPgX8DrvDWv6Wqr6ctDICqPhl0Fj+Fi/W/wCJ0jChm3RuZo6TSIauqCmwOFgcFU02cRaV8+G8AbwCfBBCRfYBdgCEiMiStAiiqOjayPA+Yl0bbRuORC8Gs0jqK5s8fuWSjWZ01YPCiZUUzyrYzaXXaikgHbozTwcB/quqyEodURKmBV6Ewfykiq4HngV8Ca3CWv2EYRvuiCSfYS0SWe9PMvGZUe1X1SFyf5TEi0lULcZN22l4NvA+4T1WPEpEPEFj9hpEFcdWxzLpPFz/rpdGfMgugvKaqE0vtpKqbROR+YCrQU7FwBUhk4QPbVHUjMEBEBqjq/wJHpi2MYRhG05B00FUJP7+I7C0iw4L5wcCJwNO1EDmphb9JRIbgQoduFpFXaJHkaY2In+kRLMGbYTQqKUXp7Ad8L/DjDwBuU9W7U2k5QlKFfxougdosXHDuHsA/10Igw0iCuW9qy42dB3PP+scBOPmkMzOWpnFJo9NWVVcAR1XfUmmSFjF/21v8Xo1kMQzDaB4UaKUShyLyFvHxoIILHx1aE6naHHPhGI3Ah/c/AoCOmHgRe8MKaC59XzIOf/di2w3DMCBMxeA4fdyjGUpSX1oqeZphGIafVsEs+wgNlPo4CabwDcNIjNUK9tDmK4BiCt8wjH4s2XF7nnIPc+EXop3cOCFu4JVZ+IZhGO2BWfiGYRjtgVn4hmG0BKGP3vz2BWjhileGYRhGHoq00sCrWiMinwP+AZeX56eqelmw/krgXKAXuFBV78lOSsNob8yqL4K5dJIRpFg+DZigqluD4iqIyGHAWcB7cCUO7xORQ6zqlWEYDUUThmUmTY9cC84DrlHVrQCq+kqw/jTgVlXdqqrPA88Ax2Qko2EYRmFSSI9cT7JU+IcAfyEiy0TklyLy3mD9SOD33n5rg3X9EJGZYQWZV199tcbiGoZhREhe8aohqKlLR0TuA0bEbLoqOPe7cJW03gvcJiIH4cYzRIm9ZKraDXQDTJw4sYEuq2EY7YCFZXqo6omFtonIecAdQcX2h0RkB7AXzqI/wNt1FLC+lnIahmGUjQK9pvCT8hPgg8D9InIIsBPwGnAXcIuIXI/rtO0EHspMSsMwCrJjwyGx6weM+F2dJak/gpqFXwbfAb4jIj3AO8A5gbX/pIjcBjyFC9e8wCJ0DMNoSEzhJ0NV3wFmFNg2D5hXX4nSp5D1A7W1gIqdtxFoB+uvlWn0+6uupKDwReQA4Pu4/s4dQLeq3lB1wzHYSFvDMFKh7R7kSlrJ07YDl6jqIyKyO/CwiCxR1adSad3DFL5hGEaFpOHDV9WXgJeC+bdEZCUuFN0Ufpak+SpbTlu+5dQKr9NJP0PbWYwNTCvcd+mjsCPdobYiMhY4CliWasMBpvANw6iKtn0wK+X48PcSkeXecncwjiiHiAwBFgIXq+qb6QiZT+so/G09LWuFtOrnKkXc525b5VJnyrnndmw4pH2/l+QG/muqOrHQRhEZhFP2N6vqHSlIFkvrKHzDMIw6k4YPX0QE+DawUlWvr7rBIpjCN5qKrEJdjeI0kpVf1zfidOLwjwP+GnhCRB4L1n1RVX+WRuM+pvANwzAqQRV6q++0VdWlxOcQSx1T+IZhGJViI20Nw0hKte6H0I3Sbh37DfN5TeEbhtFuNIr/vq4oYDVtDaP9yMribBhLt8Y05udU0OaqcWgK32ha2tKqNBoLc+kYRuvQmJZlc9Jy11JJJUqnnpjCNwzDqJQms/AzK2IuIkeKyIMi8lhQiPwYb9uVIvKMiKwSkQ9nJaNhGEZh1Cn8JFODkKWF/xVgrqr+XERODpZPEJHDgLOA9+BKHN4nIodY1SsjSsu5CJqYtvwulNSzZdaazCx83OUaGszvQV+h8tOAW1V1q6o+DzwDHBNzvGEYRraYhZ+Yi4F7RORa3IPn2GD9SOBBb7+1wbp+iMhMYCbA6JHWHWEYRp1pIGWehJpqSRG5D1enMcpVwIeAWaq6UETOxGWLO5H4nBKxVzXIJ90NMPGIXZrryhuG0dyoor3N5WmuqcJX1RMLbROR7wMXBYu3A98K5tcCB3i7jqLP3WMYhtE4NNlI2yx9+OuB44P5DwKrg/m7gLNEZGcRORDoBB7KQD7DMIzimA8/MZ8FbhCRgcCfCHzxqvqkiNyGK+C7HbjAInQMw2g4NP2atrUmM4Uf5ID+swLb5gHz6iuRYRhGmTSQ9Z4EC20xDMOoCOu0NQzDaA8sPbJhGEYb0WTpkbOM0jEMw2haFNAdmmgqhYh8R0ReEZGeWspsCt8wDKMSNCiAkmQqzU3A1NoKbC4dwzCMiklivSdqR/VXIjI2lcaKINpkYUWFEJFXgbeB17KWJcJemExJaUS5TKZkNJtMY1R172oaF5HFwTmSsAtuvFFId5Aaxm9vLHC3qnZVI1cxWsbCV9W9RWS5qk7MWhYfkyk5jSiXyZSMdpRJVWvugkkb8+EbhmG0CabwDcMw2oRWU/jdpXepOyZTchpRLpMpGSZTFYjID4HfAONFZK2InFuT87RKp61hGIZRnFaz8A3DMIwCmMI3DMNoE1pC4YvIVBFZJSLPiMgVGcuyRkSeEJHHRGR5sG5PEVkiIquD/++qsQz9hmkXk0FErgyu3SoR+XAdZfqSiKwLrtVjInJynWU6QET+V0RWisiTInJRsD6za1VEpsyulYjsIiIPicjjgUxzg/VZ31OF5Mr0vmpoVLWpJ6ADeBY4CNgJeBw4LEN51gB7RdZ9BbgimL8C+HKNZXg/cDTQU0oG4LDgmu0MHBhcy446yfQl4NKYfesl037A0cH87sDvgnNndq2KyJTZtcLVmR4SzA8ClgHva4B7qpBcmd5XjTy1goV/DPCMqj6nqu8AtwKnZSxTlNOA7wXz3wM+VsuTqeqvgNcTynAacKuqblXV54FncNe0HjIVol4yvaSqjwTzbwErgZFkeK2KyFSIesikqro5WBwUTEr291QhuQpRF7kamVZQ+COB33vLayn+A6k1CtwrIg+LyMxg3b6q+hK4HzSwTwZyFZIh6+v3DyKyInD5hC6BussUDGs/CmclNsS1isgEGV4rEekQkceAV4AlqtoQ16mAXNAg91Wj0QoKX2LWZRlrepyqHg18BLhARN6foSxJyPL63QiMA44EXgKuy0ImERkCLAQuVtU3i+0as64mcsXIlOm1UtVeVT0SGAUcIyLF8r3U7ToVkKsh7qtGpBUU/lrgAG95FLA+I1lQ1fXB/1eARbhXxpdFZD+A4P8rGYhWSIbMrp+qvhz8YHcA/03f63XdZBKRQTjFerOq3hGszvRaxcnUCNcqkGMTcD8ulW/D3FO+XI1yrRqRVlD4vwU6ReRAEdkJOAu4KwtBRGQ3Edk9nAdOAnoCec4JdjsHuDMD8QrJcBdwlojsLCIHAp3AQ/UQKFQWAdNw16puMomIAN8GVqrq9d6mzK5VIZmyvFYisreIDAvmBwMnAk+T8T1VSK6s76uGJute4zQm4GRcNMOzwFUZynEQLgrgceDJUBZgOPALYHXwf88ay/FD3KvsNpxVc24xGYCrgmu3CvhIHWX6AfAEsAL3Y9yvzjJNxr3SrwAeC6aTs7xWRWTK7FoBE4BHg3P3AP9U6r6u0/dXSK5M76tGniy1gmEYRpvQCi4dwzAMIwGm8A3DMNoEU/iGYRhtgil8wzCMNsEUvmEYRptgCt8wDKNNMIXfRojI5tJ7ld3mqRKkpBaRj4nIYRW0cb+ITCxz/1UicmrMtrHipWBudUTkUyKyv7d8s4i8LiIfz1IuozExhW9UhareparXBIsfw6WgrQdnq2pNR1SLSEct20+JTwE5ha+qZ5PRSHOj8TGF34aI46si0iOuWMsngvUnBNbzj0Xk6cBalGDbycG6pSLydRG5O1j/KRH5DxE5FjgV+GpQdGKcb7mLyF4isiaYHywitwbZDH8EDPZkO0lEfiMij4jI7UESsVKf58/EFcH4DXCBt74j+Jy/Dc71d8H6ASLyDXFFM+4WkZ+FFrG4Ajb/JCJLgTOCz7FYXPbT/xORdwf77S0iC4O2fysixwXrj5e+whuPhqk2Csj9BU+2ud76nwTne1KCjKvBZ7nJ+85mBTJPBG4Ozje40LkMA2iN1Ao2JZuAzcH/6cASXPGYfYEXcYU3TgDewCWVGgD8BjfUfxdcWtkDg+N/CNwdzH8K+I9g/ibg49757gcmBvN7AWuC+c8D3wnmJwDbcYprL+BXwG7BtssJhstHPkeu3WB5BXB8MP9VgiIrwExgdjC/M7AcV/ji48DPgs84AvhDKDeugM1lXtu/ADqD+UnA/wTztwCTg/nRuNw3AP8PlzEVYAgwsMB3cRLQjcvgOAC4G3h/sG3P4P9gXMqA4cCf4dL/hscPi7sWcd+DTTaF00CMdmQy8ENV7cVlPPwl8F7gTeAhVV0LIC7P+FhgM/CcuqIR4BT+zH6tJuf9wNcBVHWFiKwI1r8P5xL6dfBisRPuoVMQEdkDp/x+Gaz6AS41NTilOsHzZ++BS5g1GbhdXTbFDSLyv5FmfxS0PQQ4Frg9kAfcgwNcoq7DvPVDA2v+18D1InIzcEd4LWM4KZgeDZaHBLL9CrhQRKYF6w8I1q8CDhKRfwd+Ctxb7LoYRhym8NuTuLzgIVu9+V7cPVJs/2Jsp89tuEtkW1wSJ8FZsZ8s4xxSoK1w2+dU9Z68lSKnlGjz7eD/AGCTunzrUQYAf66qWyLrrxGRn+ISnj0oIieq6tMFZPs3Vf1mRLYTcA+TP1fVP4rI/cAuqvoHETkC+DDObXUm8JkSn8Mw8jAffnvyK+ATgV94b5zFXSxN7NM463JssPyJAvu9havDGrIG54oA50bxz382gLiCFROC9Q8Cx4nIwcG2XUXkkGIfRF0e9DdEZHKw6mxv8z3AeeLyyyMih4hLW70UmB748vfFubLi2n4TeF5EzgiOl0DpgrOw/yHcV0SODP6PU9UnVPXLOBfSuwuIfg/wmbCPQkRGisg+uLeQPwTK/t24tx5EZC9ggKouBP4RVx8Y+l9zwyiIKfz2ZBHO7/048D84n/WGQjsHVuz5wOKgM/NlnK8/yq3AF4LOynHAtTiF+wDOPx9yIzAkcOVcRvCwUdVXcX0CPwy2PUhhhenzaeA/g05b3+L+FvAU8Ii4UM1v4t5YFuJSNIfrlhX4POAeIOeKSJjyOqyXfCEwMehwfQr4+2D9xUHH6uOBLD+Pa1RV78X1A/xGRJ4AfoxT3IuBgcHn/5fgGoArxXd/4Ga7CbgyWH8T8F/WaWskwdIjG4kQkSGqulmc0/o/gdWqOj8jWe4HLlXV5VW0EX6e4bgHznHFHnrNhIjchOtU/3HWshiNhVn4RlI+G1iXT+LcDt8ssX8teR24SWIGXpXB3cHn+T/gX1pI2d8MHA/8KWtZjMbDLHzDqCEicjgucshnq6pOykIeo70xhW8YhtEmmEvHMAyjTTCFbxiG0SaYwjcMw2gTTOEbhmG0Cf8fRGtNsr6doHUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ds_soiltype_o.where(ds_soiltype_o > 0.0).plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.0 0.0\n"
     ]
    }
   ],
   "source": [
    "count_lgs_odata,count_ogs_ldata = count_mismatch(ds_grid_spec.AREA_LND, ds_soiltype_o)\n",
    "print(count_lgs_odata,count_ogs_ldata)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "115.0 251.0\n"
     ]
    }
   ],
   "source": [
    "count_lgs_odata,count_ogs_ldata = count_mismatch(ds_grid_spec.AREA_LND, ds_soil.soiltype)\n",
    "print(count_lgs_odata,count_ogs_ldata)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    "ds_soil['soiltype'] = ds_soiltype_o"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x7f1c49d0d710>"
      ]
     },
     "execution_count": 59,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAELCAYAAADawD2zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnX28VWWV+L+LCyqKSOILCgKKV9Ku+DIkjTJpJUbaaEiajcxYOTGjTippviQzxIzMWKlkM5PTnV6s0ExD0p8Vis1YQyaGb3gVCV/QAPEFQ8UI4bJ+fzx7n/ucffc5Z59z9jn7vKzv/ezP3a/PXmeffdZeez3rWUtUFcMwDKP1GZC1AIZhGEZ9MIVvGIbRJpjCNwzDaBNM4RuGYbQJpvANwzDaBFP4hmEYbYIpfMMwjIwRkYtEpEdEnhSRi2t1HlP4hmEYGSIiXcBngWOAI4CPikhnLc5lCt8wDCNbDgUeVNU/qup24JfAtFqcaGAtGs2CvfbaS8eOHVt8p209ffODumoqj9E8rH5qXW6+87CRGUpi1IuHH374NVXdu5o2PvyB3XTj673Jzrdi65PAn7xV3araHcz3APNEZDiwBTgZWF6NbIVoGYU/duxYli8vfo12bDgkNz9gRE2up9GETJ0wOze/ePnVGUpi1AsReaHaNja+3stD94xOtG/Hfqv/pKoT47ap6koR+TKwBNgMPA5sr1a+OMylYxiGUQEK7Ej4V7It1W+r6tGq+n7gdWB1LWRuGQs/CQNG/C5rEYwGZPEKs+qN8lGUbZrMpVMKEdlHVV8RkdHA6cCfp9JwhLZS+IZhGGmSxHpPyMLAh78NuEBV/5BWwz6m8A3DMCpAUXpTSi+vqn+RSkMlMIVvGIZRITtornoipvANo03YseEQTj7pTH52721VtZNGG2mSVd+cAr2m8A3DMNoDs/DLRERmAX+Le2A+AXwa2BX4ETAWWAOcWatODD82P45KrIfDL5kPwONfuLFfO6XOV+k5oxQ7j0UrtTb+uIJGssRrRdy9Xo97XIFtTVYiNlOFLyIjgQuBw1R1i4jcBpwFHAb8QlWvEZErgCuAyzMU1TCanqgrplzXzMknnQlk9xAJz5+lDD6KmkunAgYCg0VkG86yXw9cCZwQbP8ecD9VKvyoFdB381bTanMyZcAZBbct2XF7HSUxqsG35I0MUOhtLn2frcJX1XUici3wIi6HxL2qeq+I7KuqLwX7vCQi+8QdLyIzgZkAo0cnG+JcD564blYwN6vfNnOnGPXEt4rj1iW12huho7bU+ev923IjbZuLTFMriMi7gNOAA4H9gd1EZEbS41W1W1UnqurEvfeuKg+SYRhGmQi9CadGIWuXzonA86r6KoCI3AEcC7wsIvsF1v1+wCtZCllrfHdT1JKqt9USdROkkXYgbNNSGBQnDJtMgl3L7FFgh7l0yuJF4H0isivOpfMhXFrQt4FzgGuC/3dWeyJfcU6dMLvhXg/TppD8UyfMpqMLentW9dvW0TW+1mIZNWTxiqtj/fqbO4fl5qMdn0kfMFm7cwoRyp/FA1CBd5os/2TWPvxlIvJj4BFcOtBHgW5gCHCbiJyLeygU7mVsEsJQzZFLNgLp/4CKdeD5yr1cpT51wuyqfky+XEk6Gc1yTUbctbRO3PqzQxvHXZOErC18VHUOMCeyeivO2q8JTqk0pmIJHwS1fsMwa77x+PD+RwDu4Zzm9zN40TIAtkyblGftQ7zhUeje86O77ln/eGryNStupK0pfCNCaN2XQyNZa+XIYhZ6eRy6aC4Ao6b3lNizfuQXCmpu12YtUYRec+kYSTj5pDMLKse0lX2W1nwhv3LcfkY+vT2rKv7u4vpoQks/5MOLjshbLmW1JxmjUcu3APcGFGVVpmNHzKVjGIbRBijCO9qRtRhlYQo/IHS7jJj/ABtmHVtwv75BVcnbLIdGcuVUQqXy1+otYOqE2Tlrt1xL0FxZ/Udl20jsPtzAK3PpGDGE0TnrpgwHyntwlGLmop/SPe0UZi76KQDd005JrW2jPnR0je/nhomuK8e94+8b596JI95lUjlhe2m7dsLP1ggPWeu0bQJic8kUsep9ola7r7j9bXPP/z7wNwXbuePZowC4sfPg3BtF+FAw4qnV208a7RbLT7Rl2iQANs54G4Axc7b1bZw7qOpz15Pwc5qlD6pCr5qF39QUc+eUy9zzv8/p4x51yv18t27ON/4meBikR2jR19qy39w5jCGrN9X0HLWg1p3WSS3oOHzlHydntbJXI1slLNlxe9GHXzWEbwonn9Q4IcU7UrLw49LEq+qfUmnco60Ufs5Pn2KboaXu6G/RH37JfOaen+IJjYqZOmF2zv1lVIf59sM4/Oot/CJp4m+quvEIbaXwC7Fh1rGMmP8AAOetfiZv25xvxLtl5p7//YLbIN9Vc9WKjzFvwk/ij5tF7txErLlqwvLKIRyMU8p6b0brvhr861/Mlx73Hfn7h+GQY1aX/11WYqHH9QekyT3rH0/d35+EcEzA4hV1P3UsirBNU1OhcWniU8cUvmEYRoX0phCHXyhNfNUNx9A2Cn/KgDMSd8wahtF4hCOATz7pzIJvXJWG4FZCmSNt9xKR5d5yt6p2Q7808ZuA20VkhqouSFVgQLTJajIWYuLEibp8+fKC26M+xySds2HnajHXTTGSRt2EN2m1r+K1cP+8UCSKJC/apIEJw1WhfiGrSb5H//teu7Cr3/Zap1tI0/XjK9job63SsMxoWoc0Fb6IPKyqEysSLGDc4bvpvy46LNG+Z3UuL3g+ETkDmKqq5wbLfwO8T1VT7/1raQu/VtEChmEYaXXaUjhNfOq0tMKvliSWfXRAlb8uSc75UhZWtBP5xs6DS8pULXFWfRJrPjyuXpZ/0s7m7mmn5AbpdNO/wEsacfj+W5r/v9j3628bNb2nn5W/dmFXTax8/x4s936MI2pNh8tTBpyRyLpPmqCtb7/ZeXL7ef3rmTdKkbR8+IXSxKdOyyp8s+4Nw6glqqQWpVMgTXzqZO7DF5FhwLeALtxb0mdwScF/BIwF1gBnquofirUT9eH7Cj/OXx+GQlY60Crqn680hDJqUfkW/Y2dB+eWo77nevili/nvfZrBl59kGH5SSz/8rsPv7p71j+cszN6eVf2Wk7J2YRcDlw4FYI81vbwxti8x1x5revtlu6yEqDUe/czlyJukY3THhkNKpljIIgVzGj78sV276+w7jky072fHL636fGnQCOOCbwAWq+q7gSOAlcAVwC9UtRP4RbBsGIbRMCjQqwMSTY1Cpi4dERkKvB/4FICqvgO8IyKnAScEu30PuB+4vP4SGu1As2cozYpw5HqaiQCbDSuAUh4HAa8C3xWRI4CHgYuAfVX1JQBVfUlE9ok7WERmAjMBRo8eXfAkUffNiPkP5LlfRi7ZmFf6LcmI0hfmDiqZByUJ/Y97Juei6ehyy3F0Tzsl9/p93urapAvI0lUTuqxOH/cokE093GJuutB1A/lhguE6v3RgoVGpvlsoPM/2yW8CsHEyDF+wW27fSt05pTJLRteX6vtasuP2stJ+DxjxO5bscPOhe8d37TRzRS1FrABKBec/Gvhc0FN9A2W4b4KBC93gfPi1EdHIksnTrwVcVXujcsIHZrGH4uGXzC+ZZypO2R9+yfyGtvLjMtOmgZJep229yFratcBaVQ3Nlx/jFP7LIrJfYN3vB7xSbsN+fhxwlo6/HFpWI3EWULl5YuKs32KDQZK8AWzuHNbPUi/UWds97ZRcm93TGid7YFLiFI9vxXdPO6Vfwe1weenCS+siS6HvLFw/YMTtwfH99/HrwhbKTd/bsyrvrXPMnPGR0Fbvniwjv/2WaZNqkveoWmXpW/utgVg+/HJQ1Q0i8nsRGa+qq3ADDp4KpnOAa4L/d2YoptFGtLo/f+qE2XkP0rQeno1q5R+6aC7zJtSmbQV2NFCHbBKytvABPgfcLCI7Ac8Bn8ZFD90mIufiRqGVHVT/xHWzILgBQ79k3GCYqPUd9c0XI65CUbH1ccRZsaF1Ft0GcP1lZ7uZzr51pay5RsljX8rH7m+fOmF2nsy9Paugc1LNZIuj2PdY6rPk+6b7HiJRv30uvDdYP2aOO0ex76xYSoSw2EolPHHdLKZ4b8FRi97PKFtpupEs2TJtUu56pzWczSz8MlHVx4C4+NQP1VsWwzCMpKiKWfhZ8fSzLzN5+rW5V9TDL5mfi3gYlaVgRg7fXdII9UhDWt2NUwy/gM+dG4+uqi0/wqcRCqJMnTC75iUkGynGPgkto/BLUSrELiR8pU5yXLnrfeLcNaW2DV60LPfKPmT1ptx+cTllfJdAI7hzyiX8TH444uBFy+qqSDq6xqfyYArbCB8svusvmtcmvBcH99CvIE5UtjB/0/bJb+bckOE9EncPhdfy84viQ30hP3S50LY534h39TQiY+ZsoxsX+DCEdH8HrgBKR+kdG4iWU/jlxAgb2ZEkTNAwKiHaZ1crXKet+fAzoWPr9rzcIyPmPwCB7u/oGp/XGRvtOIt24tabQhZ4v467rvF9Fq93M8cdn5ZVnzQjZaFasdXUkI0ONqqVdV/vh05omY+kz0oPOxN9C7tUOG+uw3dJ37qOrvF5b39Q3qCtoz/hujPXz48fsBguN7JVH1JNQEbic9hIW8MwjNanGUfaZp4tMy0GH7y/Hr/rZwpuL/Z0j8tlnwb+IJroEPqQqBWybsrwojLE+Wbr5aOvRXhnO7t0/M7iqKUfrqv0fvzd374LgHEXP1h0Pz+Mc/3xfcqr1HFRovd1I32v4XVeN2V4LpDj6dO/VHX2yhHv2VNn3DIl0b7XHXlbQ2TLNAvfMAyjAlTTKWJeT0zhG4aROu3QKa8I23dYlE4m7LReeeHmvpjb4Qt2y3M/FOus2T75zVxmwmpeo/2CIWPmbMuVp3thYRejprv1caF4PtHMnVGyDLGMnjv6eY3yCV19I3ri74VKOeRbQb2gEjl4/A7dcYuKtxntAA0ZuHRoqq7QtMkLiw06uJ9Oqe1mG2lbtItZRN4sMb0lIs2b39QwjJrSyoPawrDMJFOjUMrCf1ZVjyq2g4g8mqI8FfPO/pKzMtdNGc4bY+GNsX2dsf7gFN+KDy2V0HrdOGMQ6+jfgZYE38p17bqi1AOXDs3rHIsOkPIZsnpTQSu+t2dV0cE49aRY+cNCHeR+wfBWftVPSjR3UKMSl7vH/y35/yHdwIfGJp3UCiIyHlfSNeQg4J9U9WtVNx6hlMKfnqCNJPsYBocumsuYYN4Uv9EK7EjBpRNkCj4SQEQ6gHVACQdbZRRV+Kr6XKkGkuxTb0Yu2djPynxh8iAvTLLPAon6H90x/S0U3zpNQr6Fm++Xj1r1oR81HBLv+1X9VApxFZKyYviC3dg44+3csn+946z7pAXRjcpIUrGtUDbXcJtvpUdrR/j3Xlxe/EZMjRxl8YqrOXTRXLdwevXtqcK29DttP4TzrLyQdsOQsNNWRE4HvgzsA0gwqaoOrYVQhmEY1RJroKVoeJQ58GovEVnuLXcHFfuinAX8sGrhCpA0SucrwF+q6spaCVItO61XODh/nW+xwJu5Od8aHblkYz//Y9QaXTltTlWyFYu6iSYHg3zrPWqtJbHse3tW5fUZDF60LPU3go0z3u5nyedXa+oj+nZk1n765KVPKPFdx90L/aLTItE9vT2rGiIDZrWE92ZaUTpluHReKzXwKqgJcipwZbVyFSJpj8PLtVT2ItIhIo+KyN3B8p4iskREVgf/31WrcxuGYVRCDaJ0PgI8oqov10rmohZ+4MoBWC4iPwJ+AmwNt6vqHSnJcRGwEghdRFcAv1DVa0TkimD58pTOVTaNHEFhGEZ2pFwA5ZPU0J0DpV06f+nN/xE4yVtWoGqFLyKjgFOAecDng9WnAScE898D7qeEwn9n//yn6Jg521i8wnUkhQMufFdCvlvBc/dMLjyYqJyBRn7emWi2wi3TJsVmMCzX7VKsBJ/focqMrpy8L8wdlBsQVsk5C7ltwnWLV1wN04q3kdvPSI1qXXZxoZRZZI6tNbnfcAqdtqQYYy8iuwJTgL9LpcEClIrS+XQgzHGq+mt/m4gcl5IMXwMuA3b31u2rqi8FMrwkIvukdC4jA+wNyWhFFNiekoWvqn8EhpfcsUqSdtr+OxCtfxa3rixE5KPAK6r6sIicUMHxM4GZAAP33oMXvjqoaAdruO3QRXMZuNR5j8K0CqFFnDSk8IW5gxi+YLe8jtbQIgot+M1Bx+nGuV15VnWh/OSl8p+H+5Sio2t8wQydY+ZsK9qpF82IGX2rqTaFgnXYNifN2GHb27OKtQu7apb2o+UKoIjInwPHAnuLyOe9TUOBNAJQjwNOFZGTgV2AoSKyAHhZRPYLrPv9gFfiDg7CmrrBpUeG4tZkLgbXqBv+YCvDaDVaSuEDOwFDgv18l8ubwMerPbmqXkkQghRY+Jeq6gwR+SpwDnBN8P/OUm3ttF6LPsldaGSfX3v75Df7LH76HgRxaQHCZd+63jJtknsrmNGVt0+/7bhBSn6YpE851YigPF9tuG85eewL7VfqzaeUFeXvW22YazsRlye/XjSjVe+zdmFX3nLab5fNWACllA//lyKyFDhcVetpHl8D3CYi5wIvAmfU8dyGYRiJSCO1Qj0p6cNX1V4R2bPWgqjq/bhoHFR1I26IsdFkTJ5+bX6EkGG0Ktp6Lp2QR0XkLuB2PL9IinH4qVAsY98bYzsYuHRobM4P38Vw6KK5/UIP+14Fu3J586NKzc+OOWq6c9OMCtIfFXLnADXtVArx3TTVFBUPXVrhq360TyRuRLG7Xtnl8G8W4nI15Y8Ur55alKhsZMLf9eQF1wIwfEHfbyGNkbYKbN/RmkXM98RlFPugty6VOPx6M3XC7NwPKVT+UyfMtuiRMrDOb6MZyD1Ai6Q2qYaW8+GHhPH4jYxu2Zqf4S+yPa6iUCFyFn/MAKLQWhgzZ1vBEMa1C7vyQjZ9ouvGrO5f1NwnlLlYEfZyKMe67+1ZxfAFk4Lzv034BhOS92bEXPxO8VDWzZ3Oslw57dLKhW4TovnxowXNQ6Jvr3HhvOHbVlKL3m+j1QbF1fKtRltR4QejYf8dF0apwFLgIlVdW0PZak6hEM7DL5mftxxWuh81vQeKuGcMw2gvWq7TNuC7wC30RcvMCNZNqYVQtSZqIfnWc2/Pqth83yE53/0CZ0H5YXO+z3/jjLcZjns4DF60LM+P76dW2Nw5LNemPzjLlw18i3lYrs00MmD6rqxoHWDXvxBUAqO4y2vltDl5rp5WsxLrzeIVV3P4JfPz7lV/3n8T9O+DYpZ93Dq/nVbJiFmMNNNFaBN22ibtcdhbVb+rqtuD6SZg7xrKZRiG0eAIvTsGJJoahaQW/msiMoO+TG6fJK4slGEYRsbUM3dTS/rwgc8A/wHMx/nwHwjWNRSlSrgVws87A/nl3Xz3zoZZxzJq+gN568N913n7jZr+AB1d49nc6Zaj2TGj86MKVK9cu7Arl/Nn3RT3Sp97LY8UqCj1GQsRhpmG+KGVwxeQ+wzDF8QcHOnUthG06VNopG2x73rI6k2573H98cK1p+R/eTd2HlyyjVYk7eyfLZdLJ0RVX8RVYjEMw2hY6pqZVZ0fv5lIGqWzN/BZYKx/jKo2nJUPfZZLaB2NmP9AwU6uMBwtus4ntLJ9y3/E/AfYMOvYvDeAcHu4bo81fcGh0ZKDpWT3z1uKJJZaoc5f3xo0Ggv/fitVtjB3z3bmR5Ht/0vl9AsfzVt3+o6++SkDWj9rSS3z+rdqlM6dwP8B99E/xN0wDCNTsqi5oLSuD39XVc2sxGAStu27W55FD26wVRxRa7+Q9R+2VShMM7otup9vyYeZ+0ZN7+mXasHfLzrAyvfdbu4clrdvKcu+0Ohh36I3675x8UMkkyo0/41t6cL2HOxWv0pdQu+O5lL4SeOF7g5y1huGYRgBqpJoahSSKvyLcEp/i4i8KSJvicibJY8yDMOoMVn1Q6imp/BFZJiI/FhEnhaRlUHxqdRJGqWze7HtIvIeVX0yHZEqY9DLb+d3ciUk+vpX7utg6HIpldkwrsRh3KhGf7Tq9slvwhI3H76mF+r8jY62LJQ/xJVm7FveOOPtfqGZhShnXyNdotk0fVfckNWb8ordlEOlxzUaWRVcTzEs8wZgsap+XER2AnZNq2GfpD78UvyAKuvbtiOhZdLqw9mNxmDKgDPyDIY4o2Dy9Gtz8+3aB1AOaYRlishQ4P3Ap1yb+g7wTvUt9ycthV/RY05EDgC+D4wAdgDdqnpDUHDlR7gw0DXAmar6h1LtFRt4VQ1J3hxGLtmYZ+VHj4kO1CrEymlzmDonnYiD/vn8HdFc/v4Pv5SlZ8VNssPPTzR5+rW5783PdJlWjiWjNIqwI3nahL1EZLm33B3U5AY4CHgV+K6IHAE8jEtOmfqPLS2FX+lzbjtwiao+IiK7Aw+LyBLck+4XqnqNiFwBXAE0dJSQ0do049tYaK1Xaqn71n6IWf35lKH4XlPViQW2DcR5SD6nqstE5AaczvvHqgWMOVFmqOpLwEvB/FsishIYCZwGnBDs9j1c6cOSCr/efjw/pDNuAJZPdFAWEFt9CyJ50RNa+9FqVOXEJcdZ9fbDbmyWLrw09x37Fn0xn3w0HHhz57C879lX8DYgLwGaWhz+WmCtqoadcj/GKfzUSSuNW9X+JhEZCxwFLAP2DR4G4UNhnwLHzBSR5SKyfBtbqxXBMAyjPDThVKwJ1Q3A70UkfHJ/CHiqFuImTa1wHPCYqr4dZM08GrhBVV8IBH5fNUKIyBBgIXCxqr4pkuypGfjAugGGyp5NltWi+YiWNrRkac1BrUv9tTMpxth/Drg5iNB5DqhJlcGkLp0bgSOCDoXLgG/jOluPr1YAERmEU/Y3e0XRXxaR/VT1JRHZD3il2vOkTZIRuf4+aRekHrJ6U78cF75/Oa4oto+fp2fjjL4hFRZ2GU+z+O6TljaMc9ktXXhpzq3TTO6cDbOOrSgkOw3SSp6mqo8BhXz8qZFU4W9XVRWR03CW/bdF5JxqTy7OlP82sFJVr/c23QWcA1wT/L+z2nMZhtGHH+Vj1ckqQxW0gYqbJCGpwn9LRK4E/hr4CxHpgBI175JxXNDmEyLyWLDuizhFf5uInAu8SF9pxcwJO77CH0tH13gWr7iaKQPOKJhHP0qhztpihOfzLS+/Ey7OYuvbNz+6a8ycbawrUJxy44y3zVXTJMQp6mhkTbTAeanO2GbsrE9q3UfzZqVBS6ZHBj4B/BXwGVXdICKjga9We3JVXUrhGP4PVdt+vZg8/VoGZy1EBhy6aK49HIz2phUVfqDkFwJB/SNeAwrUaWo+1i7sii0gHmXLtEn9fKO+tVRscFUlVn0cUUu/mEXmb/M7XDd3Dsvl6n9jbEdZ5w99/DYAqzGJ3g++6yb6v1VIarmnH7bdWInRkpA0SuezwExgT2AcLlb+v2giKzxr/A7UevhM+w2amVGb84QPErP0Ww8/KVk0hr8ZXT81IQMLP+j7PBs4SFX/OfC4jFDVh0odm9SlcwFwDC5GHlVdLSKxsfHNgn8DD1xa3MqN3uzQP4phy7RJfbnru8bzxHWzUinKED4cDl00Ny9PfrlWWqiQD100N886H75gN97ARexsn1w6AapZ9s1FVDFPnTA7l2qjVSKyNncOy1X6ilaTi0bTpWrlpzfwqly+gUtF80Hgn4G3cJGO7y11YFKFv1VV3wnj40VkIE3nvTIMw0iZbBT+JFU9WkQeBVDVPwTx+yVJqvB/KSJfBAaLyBTgfOD/VSZre5BFyTXDqIS4nDulAhFCd0/49tu2Lp5szN5tQaSkQq7m+I7ihziSKvwrgHOBJ4C/A34GfKt8ORuHqFvDJ86FExJ15cSFS8ZRrd8+LotmsTZL/QD9z5wrtr4mWDetsAzR0bZGc+EPyNvcWXzfYr+DRsIfMBaXKTRa+jRVslH4X8cFzewrIvOAjwOJLMykUTo7RGQB8CtVzabSgNGPqRNm26AZIzXismMm5Y5nj8rNnz7u0TTEaXyUTFw6qnqziDxMX9DMx1R1ZZJjEw0TE5FTgceAxcHykSJyVyXC1ovok37DrGNzE7jBR2PmbGPU9J680MSOrvEMWb0pcafo5s5hOet+3ZThqadQ8OntWZWbKiXs+N044+2Cw+uN1mbxiqv7GQrNHqq5dOGl/e7d8Hc5Yv4DjJj/QE3qBLgyh6WnGrAr0IHT4YmHASV16czBRencDy7vQ5Dd0mhiDl00l+EUj9QwF45RLvUOQc6UHfW38EXkn3DZBxbiBq5+V0RuV9WSF7ucXDpvJM1imQUyeGf4U99y1ArOhUwCRPKHj1yysV8iMuhv9UTriEaXh6x28+umDM+dL80b3k/glVancLmpFPzwTqO58S3iqRNml0y74PdXRWsrX8/ZffvRGG8LQ1Zvcr9F+tcL4Il0ziHZ+PA/CRylqn8CEJFrgEeA1BR+j4j8FdAhIp3AhUA26ekMwzAagQS57mvEGmAX+kzcnYFnkxyYVOF/DrgK2ArcAtxDgqeJ0Rj41viYCo8zjEqotsxiYyNZxeFvBZ4MysEqMAVYKiJfB1DVCwsdWFLhB/Gec1X1Czil37CEmSvBvaIW69wMt3V0jc8bqRdHOHpvCMk7fWpdbrFSV1G0oHmlWCqF1iK8n4pF6hRy90RzTGWZSz+vPOiE2YxcsjEXSJHn1k2LbCz8ReTnMrs/6YElFb6q9orIn1UglNGEVGLVh/0JLd9B1+JUE5bZtmSj8DcCP1PVRIOtfJK6dB4NwjBvx0uu7lWoypzOw0ayeHmfwok+6YtRrHM2iv9mUCiUbY81vWyZNqlhXmOjA6bmTfgJc77xN307FBhoZRhJacQKWcUMEDdeqUqUTKJ0gLOAG4IMxt9NGoMPyRX+nrinyge9dQrUTOGLyFTgBlys6bdU9ZpanctwTJ0wG0q4fPL6A7xkbuHxZuU3J2bdV0ZaUToisgaXBK0XFxVZsNyhqs4QkaG4aJ3viogC3wV+qKpvFTtP0pG2NSmoW4ig3+A/cZ0Ra4HfishdqlpVJXfffz1mTmF/fDRdQkdXTFhXDKH1P7in8Wqg+uHZXHUAAAAf+klEQVSUV634GPPO/763NVmu/jFztqXWB2A0PsUGY0Vr5zaihV8X0nXpfEBVX0t0WtU3Awt/MHAx7j39CyLydVX990LHJc2H//WY1W8Ay1W1FvVmjwGeUdXngvPfCpwGVKXwDcMwmp0g88GncbVJfgAco6qviMiuwEqgOoWPi/l8N86HDzAdeBI4V0Q+oKoXVyp8AUYCv/eW1wL9wmhEZCauMAujR49OWYT2Icx8mMbQ8zuePcq9QUz4CQBXrfhYbptF9mSD767Jql8prLyVpFJbM1GGS2cvEVnuLXerare3rMC9gXvmm5FtUT4OzFfVX/krVfWPIvKZYkIkVfgHAx9U1e0AInIjcC/O5ZLSmLU84npC+l3a4KJ0A0ycOLHkpQ99zr09q3KjbcMiz34Ypa/4Ns54myGBntrcOazga254fKO5cgoRKuJQMYNzefl++d6eVaxd2JW33S9xGLp3wmNmLvppzeU26k/cfd/2rpyQ5HH4rxXzywPHqer6oLDUEhF5OqrQPV6KbhORL6vq5ar6i2JCJFX4I4HdcG4cgvn9g5DNrQnbKIe1wAHe8ihgfQ3OY9SRqRNm91Me1snbnoTjZPwUDU1n9SsJs9AnaEp1ffD/FRFZhHNrF1L4U4DLI+s+ErOuH0kV/leAx0Tkfpz1/X7gX0VkN+C+hG2Uw2+BThE5EFiHC0P6q7Qaj1r00UFSfSXgXN74dVNcNs091vTmtg1cOtTl4PGObQbrfuW0OXll7q5a8THmrXZW/vWXTWLmopu5sfNgAM5b/QxXrejKOz6sIzB8wW55ivuFuYPonnYKAPMW/YTrL3O5VcY0eRbGViCJIg0HF4YK2P9uywlbjlLozaBVSCNKJ9CjA1T1rWD+JFzpwuh+5+GKT40TkRXept2BXyc5V9IonW+LyM9wTx0Bvhg+kYAvJGmjHFR1u4j8Ay6FQwfwHVV9Mu3zGNkTWv2FaDqrz6iYsC+pGQynHOlE6ewLLAqSUw4EblHVxTH73QL8HPg3XFGqkLdU9fUkJ0oapSO4ZPu5KukickySKumVoqo/w1XWqprFK67O3UwAaxd2MWp6smM3znibgUtdke83xnbk5sNh2n46h2Zh8Yqrc/H0Y+Zs4/rOvkyHp497lOunuU6+6y+bxLyv3Az0+fz9SlmbO4flK+RgANfk6cNybwKfn/DToC13Dr/jrtWsvWanUIWr0A2X5HuLPrzjKsTFpR2JFh9vGlJQ+EE04hEJdu0A3gQuiG4QkT2TKP2kLp2Kq6QbhmG0IqJ1T4/8MH2PmGhvsQIHlWogqcKvuEq60VzYiEvDKIM6plZQ1QOrbSOpwq+4SnqjEpY6DAsbh+GHwxfsxvAFffsNX7Abb4ztW/Yz7tWynGGtCWPiJy+4tl9Bd79j9vrLzmbpwks5fVywQ4K8O/l+9zl59U7j/PVxIX7m7qk/vltlc8S9U6qjNuemKZV1tmt87HiPZimYHqWeFr6IvFtVnxaRo+O2q+ojpdpIqvDDKun7lFsl3WgNQsvfV87+gwJqM7Aqmviu2fpLjBanvi6dz+MGml5XQJIPxqzPI2mUjl8lXSijSnqjsGTH7bmOyjCkEsiVNgw7Y/GKHb4xtoPtk9+MNOT+NbN177N04aV5HbgvzH071zGbZkfa54PO3+5ppxS0Foes3tRPoadVytEozpIdt+cFNviEHbbhA36I91zPfZcxln2pwVluIF8VQmdNnX34qjoz+P+BStsoqvBFZE9v8RXgh/62pKFArU6oMC11gGG0GRnkwxeRQcB5uPFQ4AqgfFNVtxU8KKCUhR/2CgswGvhDMD8MeBGouhOhnvSlDehf+Sb05ZfrS+z3BtCE5B5UoX8+xk8fWvu595/OSQxfsFuiOPnTxz2am58zZT4AT1w3y6z3BmbwomW530LUQo+m4IhSqPpV1Hc/anpP0/ruc2RTAOVGYBAuehLgr4N1f1vqwKIKP+wVFpH/Au4KYuMRkY8AJ1YhsGEYRtNT57DMkPeqqh+3/z8i8niSA5N22r5XVf8+XFDVn4vIv5QjYbvRjmX/Dr/EWe/h2xL0WXW5KI4gOmrKgDNSyc5pGJmSjcLvFZFxqvosgIgchN/5WISkCv81EZkNLMB9xBnE+UVaAP9Vdo81vWycnL+9rwBI87tyyiHN4e7hA6G3Z1XJdtvpgZk10Y5bv9Pez5oKeJlSN5WVWyc/A2tXrtM2dAM11fdd/4FXIZcC/ysizwXLY3H58UuSVOF/EpiDC81UXBa3T5YnY7bU01/sW6+tbOn3i+oIrHfDaBuyUfjDgS6coj8NOJa+TMZFSRqW+TpwUYXCNR2hZeNy7vTk1ketHCOf0HLfMOvYPCveZ+SSjcnePY2GIunvYPCiZblBW4MXLcsL14zLoRPSLy9Ts5CNwv9HVb09qGs7BReXfyMxRaKilArL/JKqfqnafdqFdog6KRSrXSm1KnweTRHRlMrEaGiEzFw6oc10CvBfqnqniHwpyYGlLPy/FZFizmrB5apPdLKsCS0M3wIthZ87f+DSof3CMP3skdC+6QA6usbn+XKjqSv8POvWWdu4hH0qcQ/28HvLD6fcxMYZfYP1SoVZbpk2KW+wVVOn0FCQbBLMrBORb+IiJb8sIjsDA5IcWErh/zcuuX6pfQzDMNqPbCz8M4GpwLWquklE9iNhXZJScfhzUxAuFhH5KvCXwDvAs8CnVXVTsO1K4Fzcq8uFqnpPtecLa9dCcj/yqOk9efv51r1v2S9deGmeO8e3YFuxszaO3p5VDA7cvIPps/S2TJtUMEVDvdJTmDunfPzoqbgQWn8g3qjpyQYsbvF9+zRvwrQ8MlD4qvpH4A5v+SXgpSTHJnoNqBFLgC5VnQD8DrgSQEQOw7mJ3oN7in0jyNRpGIbRUIQ58UtNjULSsMzUUdV7vcUHcRk4wYUZ3aqqW4HnReQZXGnF39RZRMMjN6iqRm0/cd2sRPv6nbFmuRuZ00DKPAmZKfwInwF+FMyPxD0AQtYG6/ohIjNx6UIZPXp00RP4g0qSdhpumHVsXv77UdN7ePZr7wNg/fHw3IWXJGqnHQlf2zfMOhamTSqZOTFt7GGQHkt23F4wAs0Paohj8KJlud9b+N37rpx63Q81IeVO28CTsRxYp6ofTa/lPpLWtD0EF+e5r6p2icgE4FRVLeqgFpH7iDcKr1LVO4N9rgK2AzeHh8XsH/scVdVuoBtg4sSJTfasbXxCq76e50pq6YOz9k2x1x5Lg1GEdLXORcBKYGipHSslqYX/37he4G8CqOoKEbkFKKrwVbVogjUROQf4KPAhVQ0v3VrgAG+3UcD6hHImwu/ALWad+NY9wLNfex/jLvZePi7MbzOkFX8c0foBSRkx/wG2TJtU0JLbY00vb4y1LppmoNBvpaNrfC5X/pg523LfdRiqG77d+WlLfJr5oZ2Wf15ERuHi6ufhCp3UhKSdtruq6kORddurObGITAUux70p/NHbdBdwlojsLCIHAp1A9NxGi1LPtwrDqBpNOMFeIrLcm2ZGWvoacBk1Lh1bTvK0cfTVtP04CcOAivAfwM7AEhEBeFBV/15VnxSR24CncA+VC1Q1ldH4xQaVxOEP1HKU9120UkjmuinD+73xJKWQZbfHml7n6w+2JbH0m9kabHb8N9fQio8bOOX314TVsoDYqlhN/X32KfMkvKaqE+M2iMhHgVdU9WEROSEd4eJJqvAvwPnK3y0i64DncRkzK0ZVDy6ybR7u1cYwDKMhEeI7HCvgOOBUETkZ2AUYKiILVLUqHRtH0uRpzwEnishuwABVfSttQQwjxEpGGs1CGlE6qnolfeOQTgAurYWyh9LJ02I7DwIXDKp6fQ1kqiulwsp8V0TYYVsqh3ux9pqV0J2TpLO7FL7r5o1ISuVWKBnZaoS/Ed+lE7ptwvtg+IK+3Dp+B320s76pXThxNFlsYCkLP8yjMx54L65DFVxKhF/VSijDMIymIGWFr6r344qS14REuXRE5F7g6NCVE6TiTK8EUgaElkmxEMpwW2h1rp3cVdDNkGZFqGZnQ8RqHzH/gcT5UwYu9UKQY4qpG/Un/K2E/7d4A+nC5ZC48NumzohZjAZLm5CEpJ22o3FJzkLewVVbMQzDaF9aVOH/AHhIRMISh9OA79dMqhoTrd1ZyMp/Ye4gxszZxEYGxW5vdfyaAWFcbM56n9JnxZcK1/TrD0RDL/dY0xdx+8bYjopDP43aEc14OnjRstz9EPfGFk2j0dTpE0qQUT78ikkapTNPRH4O/EWw6tOq+mjtxDIMw2h8WtKlIyKjgddwRcxz61T1xVoJZhiG0dCUN/CqIUjq0vkpfR9tMHAgsAqXs75lGTNnG709qxg13Tpl48pCFnO/+AXNRy7ZyLpIR26I7+Ixd07zM3jRMghcpH5nbSuNOs+jFRW+qh7uL4vI0cDf1UQiwzCMJiDDIuYVU1E+fFV9RETem7Yw9cTPq1MqNLNlrZME+CGWvpUfvWb+NfI7xEfMf6CfdR921EY7cMOSh+WkSDayocOz4v1O2VbLhlmSVlT4kRG3A4CjgVdrIpFhVIhVw0qX3IO7FWrP1gIF2dFcGj+phb+7N78d59NfmL442RBNE9Du/nqfqN8eSl8fP+w1fEPw24kOzDKaD3/gYtTK90M4W/231KounadUNe+bE5EzaPLRtoZhGFXRogr/Svor97h1htEQ+P0IrW5lGtnRUha+iHwEOBkYKSJf9zYNpcqKV42CKYPCVHNtosdOLtBRGyWNDtshqzeVXYrR6E+hvEd+htnenlX99tsybVL79KE0mcIvVeJwPa6K+p+Ah73pLuDDaQggIpeKiIrIXt66K0XkGRFZJSKpnMcwDCNVguRpSaZGoVS2zMeBx0XkZlVN3aIXkQOAKcCL3rrDgLNwg7r2B+4TkUPSKnNoZENo8UVr1oaZSNModtKyWRkbhELhy9FO23ax7oXmy6VT1MIPassCPCoiK6JTCuefjyvc6z8DTwNuVdWtqvo88AxwTArnMloYK35uZIJqsqlBKNVpe1Hw/6Npn1hETgXWqerjYQWtgJHAg97y2mBdXBszgZkAo0ePTltEowY8cd2sXAlD8PLfp5z7vt0HzKWN77ePWvq9PatYuqI9+8IayV2ThFIunZeC2fNV9XJ/m4h8Gbi8/1F5+9wHjIjZdBXwReCkuMPiRCkgXzeuuDoTJ05ssktvGEZT08LJ06bQX7l/JGZdHqp6Ytx6ETkcl4AttO5HAY+IyDE4i/4Ab/dRuM5jw8jD3Di1445njwLOzlqMhicNH76I7IIrGbszTif/WFWr79SKoVRY5nnA+cBBEZ/97sCvKz2pqj4B7OOdZw0wUVVfE5G7gFtE5Hpcp20n8FCl5zIal1HTe1IbdTtyycZcLh6wfDxpEy1ibi4zR0qdtluBD6rqZhEZBCwVkZ+r6oOlDiyXUhb+LcDPgX8DrvDWv6Wqr6ctDICqPhl0Fj+Fi/W/wCJ0jChm3RuZo6TSIauqCmwOFgcFU02cRaV8+G8AbwCfBBCRfYBdgCEiMiStAiiqOjayPA+Yl0bbRuORC8Gs0jqK5s8fuWSjWZ01YPCiZUUzyrYzaXXaikgHbozTwcB/quqyEodURKmBV6Ewfykiq4HngV8Ca3CWv2EYRvuiCSfYS0SWe9PMvGZUe1X1SFyf5TEi0lULcZN22l4NvA+4T1WPEpEPEFj9hpEFcdWxzLpPFz/rpdGfMgugvKaqE0vtpKqbROR+YCrQU7FwBUhk4QPbVHUjMEBEBqjq/wJHpi2MYRhG05B00FUJP7+I7C0iw4L5wcCJwNO1EDmphb9JRIbgQoduFpFXaJHkaY2In+kRLMGbYTQqKUXp7Ad8L/DjDwBuU9W7U2k5QlKFfxougdosXHDuHsA/10Igw0iCuW9qy42dB3PP+scBOPmkMzOWpnFJo9NWVVcAR1XfUmmSFjF/21v8Xo1kMQzDaB4UaKUShyLyFvHxoIILHx1aE6naHHPhGI3Ah/c/AoCOmHgRe8MKaC59XzIOf/di2w3DMCBMxeA4fdyjGUpSX1oqeZphGIafVsEs+wgNlPo4CabwDcNIjNUK9tDmK4BiCt8wjH4s2XF7nnIPc+EXop3cOCFu4JVZ+IZhGO2BWfiGYRjtgVn4hmG0BKGP3vz2BWjhileGYRhGHoq00sCrWiMinwP+AZeX56eqelmw/krgXKAXuFBV78lOSsNob8yqL4K5dJIRpFg+DZigqluD4iqIyGHAWcB7cCUO7xORQ6zqlWEYDUUThmUmTY9cC84DrlHVrQCq+kqw/jTgVlXdqqrPA88Ax2Qko2EYRmFSSI9cT7JU+IcAfyEiy0TklyLy3mD9SOD33n5rg3X9EJGZYQWZV199tcbiGoZhREhe8aohqKlLR0TuA0bEbLoqOPe7cJW03gvcJiIH4cYzRIm9ZKraDXQDTJw4sYEuq2EY7YCFZXqo6omFtonIecAdQcX2h0RkB7AXzqI/wNt1FLC+lnIahmGUjQK9pvCT8hPgg8D9InIIsBPwGnAXcIuIXI/rtO0EHspMSsMwCrJjwyGx6weM+F2dJak/gpqFXwbfAb4jIj3AO8A5gbX/pIjcBjyFC9e8wCJ0DMNoSEzhJ0NV3wFmFNg2D5hXX4nSp5D1A7W1gIqdtxFoB+uvlWn0+6uupKDwReQA4Pu4/s4dQLeq3lB1wzHYSFvDMFKh7R7kSlrJ07YDl6jqIyKyO/CwiCxR1adSad3DFL5hGEaFpOHDV9WXgJeC+bdEZCUuFN0Ufpak+SpbTlu+5dQKr9NJP0PbWYwNTCvcd+mjsCPdobYiMhY4CliWasMBpvANw6iKtn0wK+X48PcSkeXecncwjiiHiAwBFgIXq+qb6QiZT+so/G09LWuFtOrnKkXc525b5VJnyrnndmw4pH2/l+QG/muqOrHQRhEZhFP2N6vqHSlIFkvrKHzDMIw6k4YPX0QE+DawUlWvr7rBIpjCN5qKrEJdjeI0kpVf1zfidOLwjwP+GnhCRB4L1n1RVX+WRuM+pvANwzAqQRV6q++0VdWlxOcQSx1T+IZhGJViI20Nw0hKte6H0I3Sbh37DfN5TeEbhtFuNIr/vq4oYDVtDaP9yMribBhLt8Y05udU0OaqcWgK32ha2tKqNBoLc+kYRuvQmJZlc9Jy11JJJUqnnpjCNwzDqJQms/AzK2IuIkeKyIMi8lhQiPwYb9uVIvKMiKwSkQ9nJaNhGEZh1Cn8JFODkKWF/xVgrqr+XERODpZPEJHDgLOA9+BKHN4nIodY1SsjSsu5CJqYtvwulNSzZdaazCx83OUaGszvQV+h8tOAW1V1q6o+DzwDHBNzvGEYRraYhZ+Yi4F7RORa3IPn2GD9SOBBb7+1wbp+iMhMYCbA6JHWHWEYRp1pIGWehJpqSRG5D1enMcpVwIeAWaq6UETOxGWLO5H4nBKxVzXIJ90NMPGIXZrryhuG0dyoor3N5WmuqcJX1RMLbROR7wMXBYu3A98K5tcCB3i7jqLP3WMYhtE4NNlI2yx9+OuB44P5DwKrg/m7gLNEZGcRORDoBB7KQD7DMIzimA8/MZ8FbhCRgcCfCHzxqvqkiNyGK+C7HbjAInQMw2g4NP2atrUmM4Uf5ID+swLb5gHz6iuRYRhGmTSQ9Z4EC20xDMOoCOu0NQzDaA8sPbJhGEYb0WTpkbOM0jEMw2haFNAdmmgqhYh8R0ReEZGeWspsCt8wDKMSNCiAkmQqzU3A1NoKbC4dwzCMiklivSdqR/VXIjI2lcaKINpkYUWFEJFXgbeB17KWJcJemExJaUS5TKZkNJtMY1R172oaF5HFwTmSsAtuvFFId5Aaxm9vLHC3qnZVI1cxWsbCV9W9RWS5qk7MWhYfkyk5jSiXyZSMdpRJVWvugkkb8+EbhmG0CabwDcMw2oRWU/jdpXepOyZTchpRLpMpGSZTFYjID4HfAONFZK2InFuT87RKp61hGIZRnFaz8A3DMIwCmMI3DMNoE1pC4YvIVBFZJSLPiMgVGcuyRkSeEJHHRGR5sG5PEVkiIquD/++qsQz9hmkXk0FErgyu3SoR+XAdZfqSiKwLrtVjInJynWU6QET+V0RWisiTInJRsD6za1VEpsyulYjsIiIPicjjgUxzg/VZ31OF5Mr0vmpoVLWpJ6ADeBY4CNgJeBw4LEN51gB7RdZ9BbgimL8C+HKNZXg/cDTQU0oG4LDgmu0MHBhcy446yfQl4NKYfesl037A0cH87sDvgnNndq2KyJTZtcLVmR4SzA8ClgHva4B7qpBcmd5XjTy1goV/DPCMqj6nqu8AtwKnZSxTlNOA7wXz3wM+VsuTqeqvgNcTynAacKuqblXV54FncNe0HjIVol4yvaSqjwTzbwErgZFkeK2KyFSIesikqro5WBwUTEr291QhuQpRF7kamVZQ+COB33vLayn+A6k1CtwrIg+LyMxg3b6q+hK4HzSwTwZyFZIh6+v3DyKyInD5hC6BussUDGs/CmclNsS1isgEGV4rEekQkceAV4AlqtoQ16mAXNAg91Wj0QoKX2LWZRlrepyqHg18BLhARN6foSxJyPL63QiMA44EXgKuy0ImERkCLAQuVtU3i+0as64mcsXIlOm1UtVeVT0SGAUcIyLF8r3U7ToVkKsh7qtGpBUU/lrgAG95FLA+I1lQ1fXB/1eARbhXxpdFZD+A4P8rGYhWSIbMrp+qvhz8YHcA/03f63XdZBKRQTjFerOq3hGszvRaxcnUCNcqkGMTcD8ulW/D3FO+XI1yrRqRVlD4vwU6ReRAEdkJOAu4KwtBRGQ3Edk9nAdOAnoCec4JdjsHuDMD8QrJcBdwlojsLCIHAp3AQ/UQKFQWAdNw16puMomIAN8GVqrq9d6mzK5VIZmyvFYisreIDAvmBwMnAk+T8T1VSK6s76uGJute4zQm4GRcNMOzwFUZynEQLgrgceDJUBZgOPALYHXwf88ay/FD3KvsNpxVc24xGYCrgmu3CvhIHWX6AfAEsAL3Y9yvzjJNxr3SrwAeC6aTs7xWRWTK7FoBE4BHg3P3AP9U6r6u0/dXSK5M76tGniy1gmEYRpvQCi4dwzAMIwGm8A3DMNoEU/iGYRhtgil8wzCMNsEUvmEYRptgCt8wDKNNMIXfRojI5tJ7ld3mqRKkpBaRj4nIYRW0cb+ITCxz/1UicmrMtrHipWBudUTkUyKyv7d8s4i8LiIfz1IuozExhW9UhareparXBIsfw6WgrQdnq2pNR1SLSEct20+JTwE5ha+qZ5PRSHOj8TGF34aI46si0iOuWMsngvUnBNbzj0Xk6cBalGDbycG6pSLydRG5O1j/KRH5DxE5FjgV+GpQdGKcb7mLyF4isiaYHywitwbZDH8EDPZkO0lEfiMij4jI7UESsVKf58/EFcH4DXCBt74j+Jy/Dc71d8H6ASLyDXFFM+4WkZ+FFrG4Ajb/JCJLgTOCz7FYXPbT/xORdwf77S0iC4O2fysixwXrj5e+whuPhqk2Csj9BU+2ud76nwTne1KCjKvBZ7nJ+85mBTJPBG4Ozje40LkMA2iN1Ao2JZuAzcH/6cASXPGYfYEXcYU3TgDewCWVGgD8BjfUfxdcWtkDg+N/CNwdzH8K+I9g/ibg49757gcmBvN7AWuC+c8D3wnmJwDbcYprL+BXwG7BtssJhstHPkeu3WB5BXB8MP9VgiIrwExgdjC/M7AcV/ji48DPgs84AvhDKDeugM1lXtu/ADqD+UnA/wTztwCTg/nRuNw3AP8PlzEVYAgwsMB3cRLQjcvgOAC4G3h/sG3P4P9gXMqA4cCf4dL/hscPi7sWcd+DTTaF00CMdmQy8ENV7cVlPPwl8F7gTeAhVV0LIC7P+FhgM/CcuqIR4BT+zH6tJuf9wNcBVHWFiKwI1r8P5xL6dfBisRPuoVMQEdkDp/x+Gaz6AS41NTilOsHzZ++BS5g1GbhdXTbFDSLyv5FmfxS0PQQ4Frg9kAfcgwNcoq7DvPVDA2v+18D1InIzcEd4LWM4KZgeDZaHBLL9CrhQRKYF6w8I1q8CDhKRfwd+Ctxb7LoYRhym8NuTuLzgIVu9+V7cPVJs/2Jsp89tuEtkW1wSJ8FZsZ8s4xxSoK1w2+dU9Z68lSKnlGjz7eD/AGCTunzrUQYAf66qWyLrrxGRn+ISnj0oIieq6tMFZPs3Vf1mRLYTcA+TP1fVP4rI/cAuqvoHETkC+DDObXUm8JkSn8Mw8jAffnvyK+ATgV94b5zFXSxN7NM463JssPyJAvu9havDGrIG54oA50bxz382gLiCFROC9Q8Cx4nIwcG2XUXkkGIfRF0e9DdEZHKw6mxv8z3AeeLyyyMih4hLW70UmB748vfFubLi2n4TeF5EzgiOl0DpgrOw/yHcV0SODP6PU9UnVPXLOBfSuwuIfg/wmbCPQkRGisg+uLeQPwTK/t24tx5EZC9ggKouBP4RVx8Y+l9zwyiIKfz2ZBHO7/048D84n/WGQjsHVuz5wOKgM/NlnK8/yq3AF4LOynHAtTiF+wDOPx9yIzAkcOVcRvCwUdVXcX0CPwy2PUhhhenzaeA/g05b3+L+FvAU8Ii4UM1v4t5YFuJSNIfrlhX4POAeIOeKSJjyOqyXfCEwMehwfQr4+2D9xUHH6uOBLD+Pa1RV78X1A/xGRJ4AfoxT3IuBgcHn/5fgGoArxXd/4Ga7CbgyWH8T8F/WaWskwdIjG4kQkSGqulmc0/o/gdWqOj8jWe4HLlXV5VW0EX6e4bgHznHFHnrNhIjchOtU/3HWshiNhVn4RlI+G1iXT+LcDt8ssX8teR24SWIGXpXB3cHn+T/gX1pI2d8MHA/8KWtZjMbDLHzDqCEicjgucshnq6pOykIeo70xhW8YhtEmmEvHMAyjTTCFbxiG0SaYwjcMw2gTTOEbhmG0Cf8fRGtNsr6doHUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ds_soil.soiltype.where(ds_soiltype_o > 0.0).plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {},
   "outputs": [],
   "source": [
    "ds_soil.to_netcdf('soiltype_grid_spec.nc')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'soiltype' ()>\n",
       "array(7204)"
      ]
     },
     "execution_count": 61,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.sum(ds_soiltype_o>0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'AREA_LND' ()>\n",
       "array(7204)"
      ]
     },
     "execution_count": 62,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.sum(ds_grid_spec.AREA_LND>0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'AREA_LND' ()>\n",
       "array(10844)"
      ]
     },
     "execution_count": 63,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.sum(ds_grid_spec.AREA_LND==0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'AREA_LND' ()>\n",
       "array(0)"
      ]
     },
     "execution_count": 64,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.sum(ds_grid_spec.AREA_LND<0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<xarray.DataArray 'soiltype' ()>\n",
       "array(10844)"
      ]
     },
     "execution_count": 65,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.sum(ds_soiltype_o==0.0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}